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1. — Introduction. 


In order to display the relativistic form of Maxwell’s equations, one usually 
writes them either in the familiar tensor form or the spinor form of LAPORTE 
and UHLENBECK (!). 

OPPENHEIMER (2) was able to cast Maxwell’s equations in a form very similar 
to Dirac’s equations. He analyzed the transformation properties by appealing 
to the usual two-component spinor formalism. 

In the present work we shall obtain a form for Maxwell’s equations which 
are identical to Dirae’s equations for particles of mass zero and which have 


(*) Prepared under Contract DA-30-069-O0 RD-1258 with the Office of Ordnance 


Research, United States Army. 
(1) O. Laporte and G. E. UHLENBECK: Phys. Rev., 37, 1380, 1552 (1931). 


(2) J. R. OPPENHEIMER: Phys. Rev., 38) 720 (19st). 
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sources (34). Rather than using the powerful but somewhat complicated spinor 
formalism to analyze the transformation properties, we shall use a simpler 
method which is analogous to that used by DIRAC (5). We shall show that 
the requirement that one component of the spinor be an invariant leads to 
transformations for the electromagnetic field different from those for the par- 
ticles of spin 3. 

The notation which we use is motivated by that used by MOLIÈRE, who 
uses a truncated set of Maxwell’s equations for another purpose (°). 

The fact that the Dirac operator is the square root of the wave operator 
will enable us to obtain in a simple manner several results usually obtained 
through the use of vector analysis. Furthermore, it will be easy to introduce 
vector potentials with the aid of spinors and to give some of tehir properties. 
We shall also give formulas by which one can obtain a vector potential from 
a knowledge of the field strengths and shall show, in the case where there 
are no sources, how the anti-potential (7) (which is the potential one can intro- 
duce by interchanging the roles of the electric and magnetic fields) is related 
to the conventional vector potential. 


2. — The Dirac form for Maxwell’s equations. 


We introduce the four component vectors 


Wo ‘Dy 
(00) 
(2.1) VE ra o li 
Wo @, 
Ws Dg 


(3) The present paper is a rewrite of the report by H. E. Moses: The Spinor Trans- 
formation of Marwell’s Equations, Report No. IMM-NYU 238, January, 1957, New 
York University, Institute of Mathematical Sciences, New York. The author’s at- 
tention has been drawn to reference (4) which has appeared in the meantime, where 
substantially the same form of Maxwell’s equations are given as in the present paper. 
However, the transformation properties are not discussed. 

(4) T. OnmuRA: Prog. Theor. Phys., 16, 604 (1956). 

(?) P. A. M. Dirac: The Principles of Quantum Mechanics, third edition (Oxford, 
1947), Chapter XI. 

(*) G. MoLIÈRE: Ann. der Phys., 6, 146 (1949). 

(7) This word is used by M. Born as quoted in B. B. BAKER and E. T. Corson: 
The Mathematical Theory of Huygens Principle (Oxford, 1939), p. 105. 
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where 
Wo = 0) (278 = OQ 
(2.14) y, = H,—th, D, = je 
py. = H, —iH, Di = 4, 
ys = H, — iB, ®, = ji. - 


The Maxwell’s equations can be written 


(2.2) SE i 


Or o A 


In (2.2) « (i=0,1,2,3) are matrices which act on yw and are given by 


uf 0 0 0 0 —1 0 0 
x 0 iF 0 0 : — 1 0 0 0 
“eo = T= 5 UA : ; 
0 0 1 0 0 0 0 — i 
0 0 0 I 0 0 a 0 
(2.3) 
0 0 —1 0 fio) 0 0 —1 
, 0 0 0 a 0 0 — è, 0 
x° = a3 == 
n a AO 1 = WIE ARL OO 
0 —èi 0 0 —1 0 0 0 
Furthermore, in (2.2) and later we write 
cu ay, agg 


It is clear that the matrices «* are Hermitian and that 


ate? = ia = — ao} 
(2.5) 
dle 
Pia = Wadi ae 
CO A eS 


Hence the matrices «’ satisfy the following anti-commutation rules 


(2.6) ata + gigi = 2041 , (4, 7 = 1; 2, 8) 
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where 6% is the usual Kronecker 6: 6%=0 for i#j, 6#=1. Thus the mat- 
rices « are a representation of the Dirac operators. 
In particular, when there are no sources, Maxwell’s equations can be 


written as 
Lee 
Se hy ask 
(2.7) i E 
where H, is given by 
13 , Cc 
(2.8) laly = ala pur 5 


Hence, Maxwell’s equations without sources are identical to Dirac’s equations 
for a particle of mass zero. If one solved Dirac’s equation for a particle of 
mass zero with initial conditions such that the top component of the spinor 
were identically zero for all time, one could interpret such a solution as a 
solution of Maxwell’s equations. 

Equations (2.6) lead to the usual result that the Dirac operator is the 
square root of the wave operator: 


RSI ie eG 
sa ESSE =| ea a = 


1 3 0 1 5 3 a 
= Lee = Tale 
| a Be sat) ( a 22 so) =e [v: 2 i 


On applying the operator — (1/i) > «‘'(0/0x;) to both sides of Maxwell’s equa- 


I=0 


wo 


tions (2.2), one obtains 


2.10 UNA 
( ) Dey fe peo 


By identifying the components on both sides of (2.10) and using w= 0, one 
immediately obtains the equation of continuity 


0, 20 a 
Aah aay LV ey 20) 


(2.11) x 


as a necessary condition that Maxwell’s equations have a solution.’ One also 


sees that if there are no sources, the electric and magnetic field components 
satisfy the wave equation separately. 
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3. — The transformation properties of the fields and sources under a Lorentz 
transformation. Spin. 

We shall now discuss how py and © transform under a proper Lorentz 

transformation and why the Maxwell spinor transforms differently than the 


usual Dirac spinor. 
Let us introduce the Lorentz transformation. 


3 
(3.1) a= Vaia, (i =0, 1, 2, 3). 
j=0 


Then Maxwell’s equations (2.2) take the form 


il 3 
(3.2) —-Ya-y=—-479, 
ren cl 
where 
3 
(3.2a) i we Shee 


In order that Maxwell’s equations be invariant we must look for a y’ and 
a ®' which are related to y and @® such that (3.2) leads to 


RENE 
(3.3) — ida y=®. 


V j=0 


It is important to note that we must have 


LÀ 


(3.34) v= 0. 
Let us make the Ansatz that py transforms as 


(3.4) y= Sw, 
where S is a matrix, 


(3.40) 8=(8); (0) 201200) 
3 
Since there is no a priori reason to expect ® (or equivalently >Y «4(0/0x7)y) 


j=0 
to transform in the same way as y, we assume a different transformation law 
3 
for D (or > «/(0/dx)y), namely 


j=0 


(3.5) D'= TO, 
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where T is a matrix 


(3.54) SEA, (1, 7 = 0, 1, 2, 3). 


On substituting (3.5) and (3.4) into (3.3) and comparing the result with 
(3.2), we obtain the result 


3 
(3.6) S = o'Ta! = Yd'jaiTai, (6 ==0) tenor 


j=0 


as equations for S and 7. It is to be noted that (3.6) constitutes four equa- 
tions, all of which must be satisfied by the same matrices S and 7. 

Equations (3.6) are not sufficient to guarantee the invariance of Maxwell’s 
equations, for we must also use equation (3.3a) which states that the top com- 
ponent of the Maxwell spinor must be an invariant. This equation can be 
interpreted as a condition on S, namely 


(3.6a) So Sas Ses 0 . 


We are now in a position to see how our transformation differs from the 
usual Dirac transformation. Equation (3.6) also holds for the usual Dirac 
equation. However, instead of condition (3.64), one has another condition 
involving the 6 matrix: 


(3.6D) S = pTe, 


In the Dirac case (3.65) leads to a solution for 7 in terms of S 


(3.66) TESI 


where the carat means Hermitian adjoint. On substituting (3.6c) into (3.6) 
one gets the usual S for the Dirac equation. 

For the Maxwell case, however, we shall solve (3.6) for S and 7 under 
condition (3.64). We obtain the following theorem: 


Theorem: Let the Lorentz transformation (3.1) be a proper one. Then equations 


(3.6) have only one solution for S and T under the condition (3.6a), the solution 
for T being given by 


(357) D5 = Aa’; , 
where 2 is an arbitrary scalar, and the same solution for S being given by any 


of the equations (3.6) where T is given by (3.7). Furthermore, the matrices S/2 
will form an orthogonal representation of the proper Lorentz growp such that 


(3.8) Soo/A => il , Sie — Sea ZI Sse = 0 . 


ni 
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Before sketching the proof, we shall make some comments on the theorem, 
taking A=1. First of all, we note that the sources ® transform like a four 
vector as in the usual theory. Secondly, since (3.6a) and (3.8) hold, we can 
see that four-by-four matrix S can be reduced to a one-dimensional matrix 
whose element is unity, and a three-by-three matrix which is orthogonal. 
Hence the matrices S form essentially a three-dimensional orthogonal 
representation of the Lorentz group. It should also be mentioned that on 
writing w;= H;j— iE;, (j=1, 2,3), it can be shown that the field strengths 
transform in the usual relativistic way. In Appendix I the matrix S is given 
expliticly in terms of a',. 

In particular, it should be noted that when the Lorentz transformation 
involves only a rotation in space; i.e. when 


t3.9) dg, a7 0, (7 0), 
(hen on using (3.6) with += 0, we have 
(3.10) S = T = (a';) 


and consequently the field strengths transform like an ordinary three-dimen- 
sional vector as they should. 

The theorem can be proved in several ways with varying degrees of ele- 
gance. We shall sketch a simple proof. 

Our proof will be based on two lemmas. 


Lemma 1. Any proper Lorentz transformation (a transformation such that 
a, >0, det ai, =1) can be decomposed into a product of no more than six ro- 
tations, each rotation involving only two axes. 

A proof of Lemma 1 is given in Reference (2). 


Lemma 2. Consider a Lorentz transformation L which can be written as a 
product of other Lorentz transformations 


3 n 
CEL 


Wie 


il 
(3.11) fight 


1 2) n 
such that equations (3.6) and (3.6a) have solutions for L, L,..., L, the corres- 


1 2 
ponding solutions for S being orthogonal matrices denoted by S, S,..., S, and 


those for T by î, tT, iis q. Then equations (3.6) and (3.64) have solutions for 


(8) H. C. Len: Quart. Journ. Math., 15, 7 (1944). 
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the transformation L, those solutions being given by 
(3.12) 


Proof: We shall prove the theorem for n= 2, since the proof can be 
extended in an obvious way for n > 2. Let a‘, and a‘, be the transformation 


coefficients for L and i respectively. Then (3.6) and (3.6a) have the form 


| 
(3.13) G0, 4, 258): 
| 


Now, 
SS = Si aja! Vai S 
1. 2 
satana) 
k 
L. = 1 2 2 
+ aia'Ta (> azar Tat) 
(3.14) = 


i at, cct l'a? (Os 02, 02 Tx) 
k 


+ igo Poet es ara? Tar) 


k 


In (3.14) we have used the fact that («‘)?}=J. If we denote the coeffi- 


cients of the transformation L by a‘; we have from L= ia 


(3.15) ai;— > Oa; , 
k=) 
and hence 
3 
(3.16) Se Vaia TT. 
j=0 


But equations (3.16) are just equations (3.6) for S and 7 for the transfor- 
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mation L which we see is satisfied by 


Le 
pears is 

(3.17) 
T=Th 


for all è. Furthermore, since S and $ are orthogonal and satisfy (3.64), S will 
also be orthogonal and satisfy (3.6a). Hence we have proved our lemma. 

Now using Lemmas 1 and 2 it is clear how one may prove that for proper 
Lorentz transformations equations (3.6) and (3.a) have the solutions given by 
the theorem. One can show by direct calculation that for rotations involving 
only two axes, equations (3.6) and (3.6a) have a solution for 7° given by (3.7) 
and the same solution for S is given by any of the equations (3.6). Further- 
more, the matrix S is orthogonal (not unitary), if A is taken to be unity. For 
an arbitrary proper Lorentz transformation, one uses lemma 1 and then 2 
to obtain (3.7) for 7 and to show that equations (3.6) have as a solution for S 
an orthogonal matrix within a scalar factor A. 

To show that the matrices S form a representation of a group, when A =1, 
one must show 


a) If the Lorentz transformation is the identity, S is the identity. One 
proves this statement by direct calculation using equations (3.6) and (3.6a). 


b) If the Lorentz transformation L leads to S, then L leads to the 
inverse matrix S-!. This statement is proved in much the same way as 
Lemma 2. 


c) A product of Lorentz transformations leads to a product of corres- 
ponding matrices S. Lemma 2 essentially proves this statement. 


| We have not yet shown that the solutions S and T of (3.6) and (3.6a) are 
unique within a factor A and are given by (3.12). The uniqueness can be proved 
by considering the operators 


n n n-1 1 


S- gs. 8, P= TIT 1, 


By methods similar to those used in the proof of Lemma 2, one can show 
(3.18) S=aiTa', (i=, 1, 2, 3). 


Furthermore 


(3.180) Soi Ba Sos 0 . 


10 H. E. MOSES 


Then by direct calculation one can show from (3.18) and (3.184) that 


~ 


(3.19) LAL SÙ 


where 7 is an arbitrary scalar. Hence the uniqueness is proved. 

Having obtained the transformation properties of y, we can now describe 
the spin of the electromagnetic field. 

Let us consider an infinitesimal rotation èg around the «# space axis. Then 
the component of angular momentum M* corresponding to this axis is de- 
fined by (°) 


1 
(3.20) Miydg = 7 oy , 


where Sy is the change in y. In general one can write 


10 
a: 


(3.21) w= 
Op 


+ 0°. 


The matrix o° is defined to be the component of the spin about the axis x*. 
In our particular case, we know that for any space-like rotation (using 4 =1) 


(3.22) S=T= (at). 


On considering small rotations about the space axes, and using (3.22), 
(3.21), (3.20), we obtain as the three spin components 


— 


0 0 0 0 0 0 0 0 
0 0 0 0 0 0 0 a 
ot == 9 (ola == 
0 0 i) ap 0 0 0 o] 
0 0 a 0 0 —i 0 0 
(3.23) } 
0 0 0 0 
0 0 —4 0 
ue 
0 a 0 0 
| 0 0 0 0 
It is easily verified that 
(3.24) ox o=— ic 


(9) P. A. M. Drrac: loc. cit., Section 35. 
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and 


(3.25) (or)*=P (a*)* (0°) 82 = tot), (i= 1). 


Furthermore, it is clear that each component of the spin has the eigen- 
values + 1, 0. Hence the electromagnetic field is a field of spin 1 as required. 


3 
Incidentally, the operator o- P= (1/i) > 0'(0/0x') commutes with the Hamil- 
j=l 


tonian H, as in the usual Dirac theory. 


‘4, — Improper Lorentz transformations. 


Up to now we have considered only proper Lorentz transformations. Since 
improper Lorentz transformations can be considered as a product of a proper 
transformation by an elementary improper one, we shall restrict our attention 
to the elementary improper transformations. We shall merely summarize the 
results. 


There are four elementary improper Lorentz transformations. We shall 
consider them in order. 


4°1. Space and time reversal. — In such a transformation we have 


(4.1) ot = — a, (t= 0, 1, 2, 3). 


We again make the Ansatz 
(4,2) 
The equations for S and 7 analogous to (3.6) and (3.64) are 


sh tes li (G0) Leo) 
(4.3) 


So Soe Sos 0 9 
for which the solution is 


(4.4) T=iI =—S, 


where 4 is an arbitrary constant. If we choose A=1, then @ transforms 


like a four vector. 
In the remaining improper transformations, the Ansatz (4.2) will not work, 
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for it leads to inconsistent equations for S and 7. Instead we must assume 
an anti-linear relationship y’ between y and we write 


y' = Sy*, (asterisk means complex conjugate), 
(4.5) d= FD, 


I 
Po 0, 
for the remaining improper transformations which we now consider. 


42. Time reversal. — Here we have 


go = — x 
(4.6) 
the me (45253) 
The equations for S and 7 are 
SUL 
(4.7) Saia (Gates 


for which the solution is 


(4.8) SETA 


2 > > | 
> > °° & 
= P-_ ©S 
= > S&S & 


If we pick 4= 1, then the sources transform like a four-vector, the electric 
field E reverses sign, while the magnetic field H remains unchanged under 
the transformation. 


43. Space inversion. — In this case 


(4.9) 
VIRA (1= 1, 2, 3), 
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and on using (4.5) our equations for S and 7 are 


oe Sui 


Se Sak (4 =1, 2, 3), 


So == Soe = Sos == 0; 


for which the solution is 


(4.11) RF] 


SM 


0 
0 
1 
0 


dr I 


On picking Z= 1, © transforms like a four-vector, E transforms like a vector, 
while H transforms like a pseudo-vector. 


4°4. Space reflections. — In this case we invert a single space axis, the a! 
axis, for example, while the other axes remain unchanged, 
The equations for S and 7 are 
S=—T 
Ni = ER 
(4.12) INR 
S =— a’ Sa'* 
So tai Soe a Sos = 0 


for which the solutions are 


COAT) 
ie eS IO 
(4.13) gee tee! 
OMO 1 eet 
(“hse OST 


If 2—1, then © transforms like a four-vector, E again transforms as a vector, 
and H as a pseudo-vector. 


5. — Vector potentials. 


It is quite easy to introduce the vector potential in terms of the spinor 
notation. We shall summarize some results. Let us denote, as usual, the 
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electromagnetic field by y where 


(5.1) i ; y= 0. 


(5.2) A 


is defined to be that vector which satisfies 
Delo 
(5.3) (7) HS Di dal , 


where y' is the same as y except for the top component: 


Yo 
Wy, 
(5.4) y =| |, 
Po 
\%; 
| oi, (v= 1, 253)s 


I 


y, an arbitrary function. 


(5.4a) | 


If A is real, then one can show that the field strengths as given by (5.3), 
and (5.4) are related to the four vector A in the usual Way, le. 


(5.5) pu_ OP __ Oy’ 


Gi 


where 


(5.6) pi= A; ; 


A SPINOR REPRESENTATION OF MAXWELL’S EQUATIONS 15 


0 aa E, sr E, uni E, 
Di 0 Dr H, 

(5.64) Ala A : È (i= 0,1122300 
E, a Ho He 0 


Tf, on the other hand A is purely imaginary, 
(5.7) PA Tee 0.5 B real: 


then the components of B are the components of the anti-potential which one 
can introduce if there are no external sources. We shall discuss the anti- 
potential in more detail shortly. 

The most general gauge transformations correspond to changes in A which 
leave the last three components of y’ unaltered, but which may in general 
alter the arbitrary function y,. 

An example of a gauge transformation which will not alter any component 
of y’ is given by 


(5.8) EA 


where A’ is the new potential and where A is a four component entity which 
obeys the wave equation: 


(5.9) Cid = 0. 


This statement is proved on using A’ in (5.3) and then using (2.9). This gauge 
transformation is the one most frequently used. 

The vector potential A can be solved for in terms of y' by using any 
Green’s function G(x, t|x’, #) which satisfies 


(5.10) =a 


Vjzo 0%; 


G(x, t|x’, 1) = d(x — x')d(t—?t), 


for then we may write 


(5.11) AG) = | [G(x t]x', 0) ye Uda at. 


In (5.11) and elsewhere x refers to the space variables only. 

That an infinity of Green’s functions satisfying (5.10) exist is proved in 
the usual treatments of Dirac’s equations. Each Green’s function will lead 
to a different vector potential A, all of which, however, are equivalent in that 
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they lead to the same field strengths and, in fact, even to the same arbitrary 
function y,. The corresponding gauge transformation is usually more ge- 


neral than (5.8). 
If one requires the vector potential A to satisfy the Lorentz condition 


dA 047 
(5.12) + div A = » a = 9) 
then y,=0, y'=y and 
ip 2 d 
= — ÎT A. 

(5.13) y Da on 

3 
On multiplying both sides of (5.13) by (—1/i) Y «(0/0x’), using Maxwell’s 

420) 


equations (2.2), and (2.9), we obtain 
(5.14) O?°A=—-47nD, 


These are the usual equations for the vector potential which satisfy the Lorentz 
condition. Our potential, however, may be complex. 


51. Relation between the potential and anti-potential. — If there are no 
sources, one can introduce instead of the usual electromagnetic tensor F' 
and potential g° (equations (5.5), (5.6), (5.6a)), another electromagnetic tensor 
fi where the electric and magnetic fields are interchanged with a suitable 
change in sign, 


Of HRR DO 
(5.15) MESE : HG data Pea ec ae 
i ’ (è, pe 0, 1, 2, 3). 
Hees pi 0 De 
Wak Thy seal Bl. 0 


(5.16) di RCA 


We can now give a relation between the potential and the anti-potential 
which gives rise to the same electromagnetic field. We shall assume both po- 
tentials to satisfy the Lorentz condition. Denoting by 9 a 4-component real 
column vector whose components are the components of the potential, and 
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by » the real column vector whose components are components of the anti- 
potential, the relation is 


: 2 e) | 1 


0 
20 dai 


j=1 


(5:17) p(x, t) = — a VR gt) Cae 
|e — x" | 


where o° are the spin matrices (3.23) and x, x’ refer to Space co-ordinates only. 


92. A real vector potential corresponding to an electromagnetic field. — We 
Shall give a real vector potential from which an electromagnetic field can be 
derived in the manner of (5.3) or equivalently (5.5). Let y be the spinor 
which gives the electromagnetic field in the sense (2.1), (2.1a). Then a real 
vector potential A(=%), which will reproduce this field is 
: ieee eae | dì 
(5.18) A(x, E) 33 x 2 pa Li (x, Di 7 
where Re means real part and « are the matrices given in (2.3). The vector 
potential given in (5.18) does not satisfy the Lorentz condition. Instead its 
gauge is such that the 0-th component is the electrostatic potential, i.e., 


(3.19) 
divA4A=0, A CAT AS As) 


The writer wishes to express his thanks to Professors K. O. FRIEDRICHS 
and W. MAGNUS and Dr. R. KRAICHNAN of New York University, and Pro- 
fessor M. MARcUVITZ of Brooklyn Polytechnic Institute for helpful discussions. 


APPENDIX 


Elements of matrix S for a proper Lorentz transformation (/ = 1). 


Sty = al . 


q Qy 
Sri Sas ‘Sas Sha lato, S 30 0 . 
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fi lg 


Re ee 


2 


0 ZO) 
o 404 


(aa 


I 


a 


Y 
IST: 


0 
2 


lin i(ala 


TG 0 
13 = (4300 — 40943 


S 


ye 


DE 


DS 

CRE a SÈ = 
3 3 8 
ee QUE I Ce COD 
S (Si Sì S Ss 
3 3 3 3 3 
SE TOSI, Toe 
3 3 3 S 3 
x_— x_— x_ — ~~ 
TS 3S 8S ‘© cs 
| I | 
ie eee 
o o SI o N o 
3 8 3 Ss 3 
aS E° al aS AS 
[sj iS 3 Ss DS 
DE of Sl SE eS 
S 3 S S S 
GS al pi eS E 
8 Sì 3 S S 
x_ x_— x_ x x_ 
ui Del = a Del 
an El = 0D, 00 GD) 
Zp) (ay Cel si (22) 
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Statistical Weight of a System of Particles with Arbitrary Spins. 


V. S. BARASENKov and B. M. BARBASEV 


Joint Institute for Nuclear Research 
Laboratory of Theoretical Physics - Dubna, USSR 


(ricevuto il 29 Marzo 1957) 


In some cases it is important to know the number of possible states of a 
system of particles with a given angular or isotopic spin. An instance of 
this is the multiple production theory (4?) the analysis on the basis of the 
charges of separate particles (*4) being often impossible because of considerable 
complexity (*). The statistical weights of the systems of particles of spin 
S =}; S =1 have been studied in details in papers (*) and (*). 

The formulae which have been obtained in these works are complicated 
and the drawing up of tables of statistical weights for different numbers of 
particles proves rather tedious. 

The difficulties considerably increase if the particles of spin S = 3 are 
also taken into account, which appears to be indispensable to bring the 
mul- tiple production theory in agreement with experiment (**). However, 


(1) E. FERMI: Progr. Theor. Phys., 5, 570 (1951). S.Z. BeLENKIJ and L. D. LANDAU: 
Usp. Fiz. Nauk, 56, 309 (1955). 

(2) S. Z. BELENKIJ and A. L. NrxiSov: Zu. Éksper. Teor. Fiz., 28, 744 (1955). 

(3) V. M. MaxsimeNnKo and A. I. Nikisov: Zu. Éksper. Teor. Fiz., 31, 727 (1956). 
(4) V. S. BarAaSENKOYV, B. M. BarBASEV, E. G. BuBELEY and V. M. MAKSIMENKO: 
to be published in Zu. Eksper Teor. Piz. 

(*) Sometimes a partial analysis of the system charge states turns out to be pos- 
sible by separating a part of the particles into groups (e.g. pions and nucleons from 
strange particles) which are characterized by their total isotopic spin. No charge ana- 
lysis is carried out within each group. The statistical weight of such a group is deter- 
mined by the number and kind of its particles and the value of the total isotopic 
spin. The isotopic spin of the group may be considered. 

(9) G. Racan: Rend. Lincei, 17, 386 (1933). 

(8) Y. YEIVIN and A. pm SHALIT: Nuovo Cimento, 1, 1147 (1955). 
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the drawing up of tables of statistical weights for systems of particles 
with different spins can be substantially simplified by using recurrent 
relations. 

The statistical weight W: _£,_.e,(8) of the system (Ej ... Gen) Lobo ione 
sists of €, particles of spin S,, & particles of spin S,, etc., for the given total 
spin of the whole system may be defined as a sum of statistical weights 
Wi «,-1..¢,(8:) of the system (€, ...€—1...&,), from which we may pass 


to one of the states of the system (é,...&,...é,), by adding a particle of 
spin S,: 


(1) Wi eee (S) = vi Wi, albp ab ries (S;) . 
ie 
It follows from the definition of the spin addition law that S; =|S — Sk: 


(S—(S,—1)|... (8 + 8,).. Evidently W, . 31..<(S,)=0 if the system 
(€,...€,...€,), cannot have the spin value S = 8;. 

Further, let us assume W,,0..0(0) =1. Ina specific case when (€,...é ...€,) = 
== (69. 45.0) = (7, 0), and (51; €y), = (0)-6,_,) = (0; 2),, relation (1) 1s identieal 
to the corresponding relations in (°). 

Below we give the tables of values Mine eS. = W,..m(S) which have 
been obtained with the help of the recurrent relation (1) for the values of 
v, m, n and which are necessary to compute the multiple production of par- 
ticles at energies E < 10—15 GeV. 

The tables do not contain the values of W,,,(S). These are presented 
in (6) (*). The calculations have shown that the use of the recurrent relation (1) 
is more convenient than that of closed formulae for W. (8) (e.g. of for- 
mulae (1), (7) and (9) in (*)) and considerably simplifies the drawing up of 
tables of statistical weights (‘). In a number of cases the calculations for 
particles of large spin are simplified by preliminary transformation of the 
statistical weight expressions on the right-hand side of relation (1) with the 
help of the recurrent relation (1). 
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In conclusion the authors take pleasure in thanking Prof. BLOHINCEY 
for his interest in the work. 

Thanks are also due to A. A. ANDREANOVA and R. A. MEDVEDEVA for 
their assistance in carry ng out the calculations. 


* ‘ NS È Tariz (3 7 = i . 
(*) The tables of the work (3) were computed anew with the help of recurrent relations 
and the necessity was found to make corrections: 


Weol2) = 3; Wia(5) = 7; Waao(1)= 5798; Wyo) = 49360;  W;,, = 21726. 


We take an opportunity to thank E. G. BuBeLEv, L. G. ZAstAvENKO and K. Hstex 
for their help in correction of the tables (*) and in discussion. 
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1. — The classical part of the theory of nuclear reactors (*) is founded on 
the Boltzmann-equation, well-known from the kinetic theory of gases. Es- 
sentially all practical problems can be reduced to the search for approximative 
solutions of that equation fulfilling various initial and boundary conditions. 

A considerable number of papers (') have been publised on the various 
methods of approximation and actual calculations performed with these me- 
thods. A few monographs (?) have also appeared which summarize systema- 
tically the results so far published. However, up to the present the processes 
occurring in nuclear reactors have not been considered in terms of the proba- 
bility theory. Certain initial attempts were made in this field as well (8) no 


(*) The non-classical part of the theory of nuclear reactors, namely, the theory 
of cross-sections etc., belongs to the field of theoretical nuclear physics. 

(1) R. E. MARSHAK: Rev. Mod. Phys., 19, 185 (1947); A. M. WEINBERG and L. C. 
NODERER: Theory of Neutron Chain Reactions, AECD-3471 (1951); O. OuLsson: Ark. 
for Fys., 10, 129 (1955). 

(2) S. Grasstone and M. E. EpLuND: The Elements of Nuclear Reactor Theory, 
(New York, 1952); H. SoopaK and E. C. CampPBELL: Hlementary Pile Theory (New 
York, 1950); Selected Reference Material United States Atomic Energy Program, Vol. 2. 
Reactor Handbook: Physics (1955). È ; 

(3) V. V. Gavranipze and 0. D. CrsSvinr: Zu. Éxper. Teor. Fiz. 28, 3 (1955). 

(4) L. PAL: Publiactions de VInstitut Mathématique de Academie des Sciences de 
Hongrie, 1, 41 (1956). 

(9) L. TAKAGS: Publications de VInstitut Mathématique de V Academie des Science 
de Hongrie, 1, 55 (1956). 
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results, however, have been obtained, except when the simplest processes are 
involved. 

One of the most general formulations of Boltzmann’s transport equation 
describing the variation in time and space of the average neutron density was 
given by L. N. Ussacnorr (5) who generalized the transport equation in such 
a way that at the same time he considered the fact of neutron multiplication 
and the role of delayed neutrons as well. Obviously the solution of this gene- 
ralized transport equation can be given so much the more only under certain 
strongly simplifying assumptions. 

In the course of developing solution methods for the transport equations 
numerous authors (7), have made use of the adjoint of the transport equation. 
According to UssaAcHorr, it was B. B. KApomrsEv who for the first time 
(in 1955) pointed out the physical significance of the adjoint equation. Ussa- 
CHOFF’s paper quoted above gave an extremely clear formulation of the phy- 
sical contents of the adjoint transport equation. He introduced the concept 
of «neutron importance » and proved that the adjoint transport equation is 
equivalent to the integro-differential equation for the neutron importance. 
The definition of neutron importance may be given in the most simple way 
as follows. Suppose that at a time t =0 the power of the reactor is zero, 
i.e. no neutrons are present in the reactor. Now, if a neutron of energy Hy 
and direction of velocity 92, at a time t = 0 is injected in the volume ele- 
ment dVr, at r,, the importance of the neutron at a time ¢ is determined by 
the average number of neutrons produced by this neutron. 

We note that the definition of neutron importance used here somewhat 
differs from that given by USssAcHOFF (5). It is worth mentioning that both 
the equations fer neutron density and those for neutron importance are in 
close relation with certain generalized forms (1°) of the Kolmogoroff-Feller 
equations (*) well-known from the theory of stochastic processes. 

The purpose of this paper is to formulate in the most general manner the 
fundamental equations describing the stochastic processes in reactors, and to 
derive from these the generalized equations for the average neutron density 
and neutron importance, as well as for the fluctuation of neutron density. 
It would prove especially difficult to derive the equation describing the fluc- 


(©) L. N. Ussacnorr: Proceedings of the International Conference on the Peaceful 
Uses of Alomic Energy, Vol. 5, p. 503. 

(*) A. M. WEINBERG: Am. Journ. Phys., 20, 401 (1952); R. EHRLICH and H. Hur- 
witz: Nucleonics, 12, No. 2, 23 (1954). 

(8) E. D. PENDLEBURY: Proc. Phys. Soc., A 68, 474 (1955). 

(9°) A. KoLmoGoRoFF: Math. Ann., 104, 415 (1931); W. FELLER: Math. Ann., 113, 
113 (1936). 


(1°) L. Pit: Zu. Exper. Teor. Fiz., 80, 362, (1956); Soriet Phys., 3, 222, (1956). 
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tuation of neutron density without the results of the general theory of sto- 
chastic processes in reactors. 


2. — First of all, let us start from the assumption that the medium in which 
neutrons move is a mixture of 7 types of atomic nuclei. Denote by N, M., 
.., N, the average number of the different types of nuclei per unit volume. 
If the spatial distribution of the various nuclei in the medium is not homo- 
geneous, the density of the nuclei in a point defined by the vector p, should 
be denoted respectively by M(o); No(Po), ---» N(po). Denote further by of, 
O;,---, 0; and oj, 03; ..., 0; the microscopic cross-section for the scattering 
and absorption of the different types of nuclei. Let w’/(H,, Q,; H, 2)dHdQ 
be the probability that the energy of a neutron which suffered elastic or 
inelastic collisions at the j-th nucleus, is inside the interval (#, E+d£) and 
the direction of its velocity is inside the elementary cone dQ about Q provided 
that before the collision the energy of the neutron was 4, and its direction of 
velocity Q,. 

It is expedient to introduce the following density function: 


> N;ojwi(E, Lo; E, 2) 
si 20s(Po; Ky, o; E, Q) = i! 


fi 
pati 
> N;0 j 


Gia 


For the sake of simplicity we introduce for the notation of the microscopic 
cross-sections the quantities 


(2.2) OS MON 
DICA î) 


Although for the general treatment several types of fissionable nuclei 
ought to be taken into account, in practice, however, usually only one single 
type of fission material is present. Thus in our treatment we consider only 
one type of fissionable nucleus. Denote by N, = N,(po) (*) the density of 
nuclei of the fissionable substance at the point defined by the vector e,, and 
by Q,= N,o' the macroscopic fission cross-section. The total cross-section is 
evidently 


(2.3) Qo a Qs "a Qa + Q; . 
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cated manner, and in case the spatial distribution of nuclei is inhomogeneous 
it depends on the space co-ordinates as well. 

Regarding the number, the energy distribution, ete., of neutrons arising 
in the fission process most general assumptions will be made. Be p(y, ko, ..., ki) 
the probability that in a fission caused by a neutron of energy #4, k, prompt 
neutrons and in succession k,, k,, ..., k, fission products in group / with decay 
constants A,, A,,..., 4: will be produced. The delayed neutrons belonging to 
the i-th group are emitted by the fission products of the groups 7. Evidently 
we have 


(2.4) darei 
Toga) 
and 
1 
(2.5) D (kot + hy) p(B, Bos +++) a) = kp. 
Hig niesss Key î=0 
where <k)> = ” is the average number of prompt neutrons per one fission 


process, and <4k,) the average number of fission products emitting delayed 
neutrons of the i-th group, also related to one act of fission. It can be seen 
directly that 


(2.6) =o, (=e we, 


v; being the average number of delayed neutrons of the 7-th group. Denote 
Dyas (i= 0,1,.,-;.0) the mean squaresof We (=O) leona). 

Suppose that both the prompt neutrons and the delayed neutrons are 
independent of each other. Let us denote by w}(H,, 2); H, 2)dHdQ the pro- 
bability that the energy of some of the prompt neutrons lies inside the energy 
interval (£, H+ dE), and its direction of velocity inside the elementary cone dQ 
about the vector £ provided the energy of the initial neutron was E, its 
direction of velocity Q,. 

The delayed neutrons belonging to the i-th group are emitted by fission 
products of type i. We denote by w‘(H,, 2); H,Q2)dHdQ the probability 
that the energy of the delayed neutron emitted by the fission product of type è 
is inside the energy interval (2, E+d£), its direction of velocity inside the 
elementary cone dQ about the vector Q again under the assumption that the 
prime neutron had the energy 4, and direction of velocity Q,. The question 
of a suitable choice for the actual form of the probability density functions 


w(E,, 2); H,Q) (i=0,1,..., 1) will be dealt with when discussing the appli- 
vations. 


2 


è. — There exists in many respects a considerable resemblance between 
stochastic processes taking place in reactors and the cascade processes of cosmic 
‘adiation. Thus it seems to be practicable to make use of the relevant results 
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of cosmic ray cascade processes as well in the theory of stochastic processes 
occurring in reactors. 

Suppose that the space and energy distribution, as well as the distribution 
of the directions of velocity of the neutrons are known at a time t = 0. 
Since each neutron initiates independently of the others a chain reaction with 
a certain probability, it is sufficient to consider a process initiated by one 
Single neutron only. Denote the energy of this neutron by H,, the direction 
of its velocity by Q,, and let it be located at the time t = 0 in a point de- 
fined by the vector r,. For the treatment of this chain reaction initiated 
by one neutron the cascade theory of Janossy (1) is extremely well suited. 
This theory makes use of that characteristic of certain stochastic processes 
(Markoff-processes) that by a suitable choice of probability variables the process 
can be determined by means of the first discontinuity of these variables (the 
so-called first collision method). 

This procedure is in fact not a new one since among the fundamental 
equations of the theory of stochastic processes the so-called first Kolmogoroft- 
Feller equation concerns exactly the first sudden change of probability variables. 
In case of processes which do not involve a change in the number of particles, 
the advantage of this method is not so evident; this advantage, however, is 
almost invaluable in the theory of the so-called cascade processes implying 
a change in the number of particles. It was JANossy (4!) who has called 
attention to this circumstance for the first time. 

In the following we also make use essentially of the results of Janossy’s 
cascade theory (!!). However, our method slightly differs from that developed 
by JANossy. One of the differences is that, while JANOSSY chooses some of 
the space co-ordinates (approximately the direction parallel to the axis of the 
cascade) as parameter for describing the process, we choose, of course, the 
time as parameter. The application of one of the space co-ordinates as a para- 
meter involves also in the cascade theory certain difficulties, as during the 
development of the process not only an increase but a decrease of the space 
co-ordinate may occur, although in cosmic ray cascades the probability of a 
backward scattering of particles along the axis of the cascade is insignificant. 
Thus the possible backward scattering may, as « catastrophic collisions », be 
excluded from the process. 

However, in the theory of the slowing down and diffusion, as well as of 
the multiplication of neutrons this method of exclusion can by no means be 
utilized since, unlike cosmic ray cascade processes of high energies, in this 
case an accurate investigation of backward scatterings is of decisive importance. 

We note that the possibility of the way followed by us has already been 


(11) L. JANOSSY: Proc. Phys. Soc., A 68, 241 (1950); L: JANOSSY and H. MESSEL: 
Proc. Phys. Soc., A 63, 1101 (1950). 
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pointed out by JANossy (22); introducing certain generalizations (7) he even 
eliminated the difficulties referred to in the case of processes implying no change 
in the number of particles, and applies the results of this calculations for 
determining the energy distribution of particles slowing down in nuclear emul- 
sions. The further generalization of the equation given by JANossy was dealt 
with by the autor (7°) as well. 

The other difference is due to the fact that, on the one hand, the cosmic 
ray cascade process is usually treated as one taking place in an infinite space, 
on the other, the processes in atomic reactors take place naturally in a medium 
of finite dimensions, and thus it should be given what happens to neutrons 
arriving at the boundary of the medium. Suppose that the reactor defines 
a convex region and the space surrounding this region is empty. The reflector 
and shield surrounding the active core of the reactor are of course considered 
part of the reactor.) If a neutron in such a system has arrived at the boundary, 
it will get out into the empty space from which it cannot return to the reactor 
any more. The boundary condition has to express this circumstance. In the 
course of the following it will be shown that the treatment applied here espe- 
cially simplifies the formulation of the boundary condition. 

Let us turn now to the determination of the stochastic state of the reactor. 
We choose as the origin of our co-ordinate system the geometrical centre of 
the region V, which is filled with a multiplying and moderating medium. Denote 
by V, the convex region around the end-point of the radius vector r. Although 
no restrictions need be made on the dimensions of the region V,, and this 
region may also include a part of the empty space surrounding the reactor, 
it is still advisable to consider the region V, as part of the region V,, occupied 
by the reactor. The neutrons should be classified at a time ¢ according to whether 
they are inside or outside the region V,. (For sake of accuracy we note that 
the neutrons being at the time ¢ just on the covering surface of the region V, 
are also considered to be inside the region.) Denote the number of neutrons 
inside and outside the region V, by m and m’ respectively. The neutrons 
both inside and outside the region may be classified according to the following 
definitions (1): 


Type 7, denotes neutrons of energy lower than £, with directions of velo- 
city lying in the cone 2, about Q. 


Type 7, denotes neutrons of energies lower than H, with directions of 
velocity lying outside the cone Q, about Q. 


(12) L. JAnossy: Zu. Heper. Teor. 26, 518, (1954); Ac'a Phys. Hung., 2, 289 (1953). 
(13) L,. JAnossy: Zu. Éxper. Teor., 30, 351 (1956). 
(4) V, A. SEVAST'JANOV: Usp. Mat. Nauk, 6, 47 (1951). 
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Type 7; denotes neutrons of energies higher than H with directions of 
velocity lying in the cone Q9 about 2; and finally 


Type 7, denotes neutrons of energies higher than £ with directions of ve- 
locity lying outside the cone Qg about &. 


Suppose that inside the region V, the number of neutrons of type 7, is m,, 
the number of neutrons of type 7, is m,, of 73 is ms, and of 7, is m, respectively. 
Denote further by the non-negative integers m,, m,, m, and m, the number 
of neutrons of the respective types outside the volume J,. 

In the following we describe the state of the system at a time t by intro- 
ducing two four-dimensional vectors with non-negative integer components, 
i.e. the vectors 


(3.1) 10 = (Gthe Wns Wl Vln) 
and 
(3.2) m'= (m,,;m,, m,, ™,) 


which indicate that at the moment t the states 


(3.3) 6,= ml, mT ,4- mt, ml, 
and 

« ~ I i I I 
(3.4) S__=MMT,+mMT,t+mMT,+ m,T, 


exist. The state S,, actually expresses that inside the volume V, at the time ¢ 
the number of neutrons of type 7, is m,, that of type T, is m,, etc. The 
meaning of G,,, is that outside the volume V, at the moment ¢ the number 
of neutrons of type 7, is m', that of type 7, is m,, ete. 

Denote by A(t, r,, 2), E; m, m') the probability that at the moment ¢ 
the neutron distribution defined by the vectors m and m' is fulfilled, provided 
that at a time ¢=0 a neutron of energy E, and direction of velocity 2, was 
present in the point defined by the radius vector rj: Our task consists really 
in the determination of the distribution function H(t, ry, ®, 4); m,m’). It 
is evident that, if we know this distribution function, we may obtain in prin- 
ciple exact answers to any question of neutronics. Unfortunately, however, 
the determination of the distribution function H(t, ro, 2), Lo; m, m') has 
proved impossible in the majority of cases. Nevertheless, it is not so much 
the distribution function itself, as especially its first and possibly second mo- 
mentum, that are needed. The momenta, however, can be most symply de- 
rived from the generating function G(t, ro, 2), Ho; g, g') of the distribution 
function H(t, ry, 2), Hy; m, m'). 


DI 
bo 


ip, led AB 


4. — In the following we discuss the determination of the distribution 
function H(t, r,,2,, By; m, m'). The first collision of a neutron of energy È, 
(denoted 6-neutron in the following) which starts at the moment ¢ = 0 from 
the end-point of the radius vector r, in direction Q,, results in either absorption 
or scattering (elastic or inelastic) or fission. A further possibility is that the 
neutron does not suffer collision in the time interval (0,7) at all. Taking into 
account these four possible and exclusive starting events of the process, we 
construct the distribution function H(t, r,, 2), E; m, m'). Thus 


(4.1) Hi ah ee Ee oe es 


where the distribution functions on the right hand side represent the following 
in succession: H, is the probability that the neutron of energy H, starting 
at the moment t= 0 from the end point of the vector r, suffers no collision 
in the interval (0, #), and at the moment ¢ the state of the system is defined 
by the vectors m and m'. H,, H, and H, are in succession the probabilities 
that at the first collision in the interval (0, ¢t) the ¢-neutron is captured, scat- 
tered, and causes fission respectively, and at the time ¢t the state of the system 
is defined by the vectors m and m'. 

We calculate now the distribution functions on the right hand side of the 
expression (4.1). Let A(V,, P,) =1, if the end-point P, of the vector pe 
is inside the region V,, and zero in any other case. It is convenient to gene- 
ralize the expressions (2.1), (2.2) and (2.3) as follows: 


(4.2) erp py Ste de, ol ERA PI) 
and 
(4.2') COMICI 155) "i Qi(Po, E)A(Ry, Le) (i= 5, a, f, 0) ’ 


We define now the probability that the neutron of energy H, starting at the 
moment #= 0 from the end-point of the vector r, in the direction Q, suffers 
no collision in the interval (0, #). Denote this probability by F(t, ro, 8), Po). 
Obviously we have 


F(t + At, Fo; Soy Lo) = Fit, ro, Qo, E){1 == oo (To + voto; E,) At} + 0(At) 


from this we get 
t 


(4.3) F(t, ry, 9, Hy) = exp |- Vo | extent, E] ar ’ 


0 
where 


pot) = r04 DLL 
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41. — Let us proceed now to determine the distribution function 
H(t, ro, 2, Hy; m, m'). We introduce the following functions: Let 
A(V,, P,,) =1, if the end-point P,, of the vector p, is inside the region V,, 
and let is be equal zero in any other case. Further, let A(Q,,2,) =1, if 
the vector Q, lies in the cone Q about the vector £, and equal zero if the 
vector £, is outside this cone. Finally we put A(#, 2)) =1 for E > B, and 
A(#, E)= 0 for E< E,. Denote the quantities A-(V,, P,), AL, 2,) and 
A-!(E, E,) as the inverse of the above quantities. Let us denote by y, the pro- 
bability that a y-neutron is of type 7,. It can be shown easily that 


| y, = A (Qg, 2,) A(E, E), ve AQ, 2,) ACE, E), 


| Y3 = A(Qg, LI AME, E), y, = AQ, LAME, B,). 


After introducing the generalized Kronecker symbols dom =i, and 
oe Bl Ò,m. We may easily construct the strongly singular SRO 


del 
function H,. We find that 


(4.6) H(t, ro, 82), Hy; m, m') = 
4 4 . 
= Oar Oar: {4 (¥., TAs a Vi ee Pte) > vi coni] Fit, To; o, Ey) b) 
î=1 Om; t=1 Dom; 


where 
Po(t) = ro + vot + 


4°2. In the case when at the first collision the 6-neutron is captured, we 
can write the following expression: 


(4.7) Hero, Sto, Lo; Mt, mM ) = Pam don bo| QL eul MER Sting Lig) bee 


43. — If the neutron is scattered at the first collision then, after the scat- 
tering, the neutron if it changed energy may initiate in the remaining period a 
similar process to that caused by the ò-neutron. Suppose that in the period 
(0, #') no collision occurs, further in the time interval (t’, t’+-dt’) a scattering 
takes place, as a result of which the energy of the d-neutron is reduced from E, 
to a value inside the interval (H', H'+ dH’), and the initial direction Q, of 
its velocity is changed into the direction inside the cone dQ’ about £'. During 
the time t— #' still available the scattered neutron initiates a process which 
at the moment t results in the neutron distribution characterized by the vectors 
m and m’. Taking into account all possible cases of the moment of collision, 


3 - Supplemento al Nuovo Cimento, 
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the changes in energy and in direction of velocity, we find 


(4.8) H,(t, To, Qd, Hy; mm = nf fe » To; S205 E) ) QLeo(t’ ), Ho]: 


‘ws [Po(t’), Hy, o; E, 2'] H[t—t’, po(t’), Q', E'; m, m']dt' dH’ dQ’. 


Considering the thermal agitation of the atoms in the moderating medium, 
it is to be expected that, in particular for small values of the energy £,, at the 
collision the energy of the ò-neutron does not decrease, but increases. 

The probability density function w* in the expression (4.8) was accordingly 
chosen. 


44. — Finally, in the case when the first collision leads to a fission, the 
further development of the process is determined by the forming of new 
« branches » initiated by prompt and delayed neutrons which were produced 
in the fission. Suppose that in the time interval (0, #) no collision occurs, 
and in the subsequent time interval (t’, ¢’-+ dt’) a fission takes place resulting 
in ky, prompt neutrons and k,,..., k, fission products. The prompt neutrons 
and the delayed neutrons emitted by the fission products initiate further pro- 
cesses and in the period ¢—t’ left they produce the neutron distribution de- 
fined by the vectors m and m’. We put K[t—#, p,(t'), 2,, E; kb; m,, mj] 
for the probability that in the period t — the k prompt aoe arising 
from the fission which occured at a place determined by the vector pol’), 
produce the neutron distribution defined by the vectors m, and m,, provided 
the energy and the direction of velocity of the neutron initiating the fission 
was respectively E, and Q,. Similarly we denote by K[t— #, e,(t'), Lo, Eo, kis 
m,, m,| the probability that in the period # — #' the delayed neutrons emitted 
by the k, fission products of type è which were produced in the fission taking 
place at the end point of the vector e,(t’), give rise to the neutron distribution 
defined by the vectors m, and m;, provided again that the energy and di- 
rection of velocity of the neutron initiating the fission was 4 and Q, respectively. 
With the help of the distribution function the distribution function H, can 
be easily expressed, namely 

t 
(4.9) H,{t, ro, 82), Hy; m, m') = nf FW, Fo, Bo; E)Q*[polt'), El: 
0 


42 PB) toy ski) > >» TL Kt, ult ), Qo, Ho, ki; mi, mi ]} dt’ . 


m=m,+.. ~+my, my Sioa! .+m,= =m’ i=0 


It can be easily shown that the distribution function K[t — t’, e,(t)’, Lo, Ey, 
lio3 To; my] may be reduced to the distribution function H(t, ro, Lo, Bo; m, m’) 
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describing the process initiated by the 6-neutron. Namely as each prompt 
neutron initiates a process independently of the others, we can write 


(4.10) K[t—t, po(t'), Lo, Hy, ko; mo, mo] = 


— SOI ih [ferma e Q'). 


my=mj+.. 5 + mp me Da j=1 


‘H[t—t', p,(t’), Q', F';m, , mj]dE'dQ'. 


We can determine the distribution function K[t — #, e,(t), 2,, E, k,; m,, mi] 
(î=1,2,...,2) as follows. In general not all of the k, fission products suffer 
decay during the period t — t. Suppose that during the time t—1?' a part %! 
of the fission products suffers decay, whereas % = 4, — ki does not. At a 
small distance from the place of the fission and with certain delay as compared 
to the moment of fission the neutrons arising from the decay initiate a process 
similar to that initiated by the first é-neutron injected into the system at 
the time #= 0. We denote by K,[t—t’, e,(t),, ®, Ho; m,, m'] the probability 
that one of the k, fission products, say the j-th, suffers decay in the time in- 
terval (#,t), and out of a point of the region around the end-point of the 
vector po(t), which can practically be regarded infinitely small, the neutron 
brought about by the decay initiates a process as a result of which a neutron 
distribution defined by the vectors m, and m, appears during the time # — #, 
assuming again that the initial neutron has the energy £, and direction of 
velocity £,. With the help of the distribution function A; we can write now 


ARA 
a 


e Kilt, pt), 2,, 25, hy: mami) = (7 es | (kE) ¢—1)]- 


k.=0 
a 


ky; 
; 5 DI i | Meer Pol t), Q,, H,; m,m,]. 


mj=my +... bmg, JOH IE COS UL =m; i=1 


Since the fission products recede into a very small distance only from the 
source of the fission, it can be assumed that the delayed neutrons start the new 
cycle (of course, with a time lag) practically from the spot of the fission act. 
Thus 


t—t' © 
(4.12) K,[t—t’, pot), Lo, Ho; my, mi] | || A; exp [— "A; |wy( By, Qo; BH’, Q’)- 
o 0 
Pipe ee Opt Qta (db de dol, 


With the help of a somewhat more careful procedure the movements of the 
fission products may also be considered. This is, however, not specially needed. 
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In this way the distribution function H(t, r,, 2), H); m, m') has been 
determined which satisfies the non-linear integral equation (4.1). It seems to 
be practically not feasible to find the general solution for equation (4.1). 
However, as has been referred to, for practical purposes it is not so much the 
distribution function itself, as rather its momenta which are required. The 
following task should thus consist in determining the generating function of 
the distribution function. 


5. — We choose two four-dimensional vectors, denoting the first by 
q= (4,444) and the second by G = (45 21 Isr a)3. WE Construct the 
scalar products 


(m-q) = MQ, + MQ: + Msd3 + Mss » 
(5.1) 


(m'-q') = mig, + mg, + ma, + mg, 
The generating function of the distribution function is given by the infinite 


series 


(5.2) G(t, ro, 2, Ly; q, q’) = DI exp [(m-q)+ (m'-q’)] H(t, ro, &), Hy); m,m'‘]. 


m, m' 


First of all we have evidently 

(5.3) G(t, ro, 2), Hy; O, O)=1. 

We are interested in particular in one of the special cases of the generating 
function; this can be derived by putting did 0 (i=2,3,4) and q;,=0 
(1=1, 2, 3,4). Introducing the notation 


(5.4) CASA ai) 


the first derivative with respect to q of the generating function g(t, ro, 9%), Po; 9) 
gives for g=0 a special, but for our purpose important first momentum of 
the distribution function. The equation 


RM: To) Lo; Ei La, E), 


(5.5) mò) = (ae Fo, Los Bo; a 


eq »; 
gives obviously the average number of those neutrons which at the moment + 
are inside the region V, around the end-point of the vector r, and which have 
vectors of velocity inside the cone 2, around the direction Q with energies 
smaller than £, provided that at the time t = 0 one single neutron of energy 
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E, and direction of velocity ®, was present at the end-point of the vector ry. 

The average neutron density ;(t, ro, Lo, Do; r, 2, E) which from a prac- 
tical point of view is of more importance, can be easily determined from the 
function M,(t, r,, 2,, Ei Va; 2g, E), as 


E 
or) =| | fuor Sa Bir 2 2) AV, da". 


Ves 29 0 


ni(t, To, Lo, Hy; r, 2, HE) dV, dQd# is evidently equal to the average number 
of neutrons being at the moment # in the volume element dV, aroung the 
end-point of the vector r, with vectors of velocity lying inside the cone dQ 
around the direction 2, and with energies inside the interval (E, E+d£), 
provided all neutrons originate from that single neutron of energy H, and 
direction of velocity Q, which at the time t = 0 from the point defined by 
the vector r, initiated the whole process. 

In the course of the following it will be shown that in the integro-diffe- 
rential equation for determining M(t, r,,2,, Bi V., La #) no operation 
refers to the quantities V,, 2,, E. Thus for the sake of convenience, in the 
following the notation of these quantities will be omitted, their role showing 
itself namely in the initial condition 


(5.6) lim My(t ry, Loy Bai Vis Las D) = AV, P.)A(0q; AB, Bo) 


only. 
So as to determine the fluctuation 


(5.7) Dit, Lato) Ve E) = M(t, Ty; LI, E.) me Mit, To) DI E) 
the second momentum which can be derived on the basis of the relation 


(5.8) mp = (POTITO 9) E oO ona 
q=0 


E 
-| | [radio Qo, Bose O, E')dVr,dQ'dE', 


V, 290 


is required. After constructing the generating function we write the equations 
for the functions M, and M,, by which both the average neutron density #1, 
and the mean square of the neutron density n, may be readily determined 
with the help of (5.5) and (5.8). The fluctuation of the neutron density is 
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given as in (5.7) by the expression 


(5.7’) eur, QI) n,(t) 7,5 Qe, B,) — Why Fo, Ae B) 


0) 


Before determining the momenta M, and M, we have to make following 
two remarks. 

Denote by Z(ro, 2), E) AV. dQ,AE, the probability that in the volume 
dV, around the end-point of the vector r, there is a neutron with energy 
lying inside the interval (E, H,+d#,) and direction of velocity in the cone 
dQ, about Q2,. If Z(r,, 2%), E) is known, the neutron density 


(5.9) n,(t, r, 2, E) =| | Jz. 92), Hy) ni (t,ro, Lo, Hy; r, Q, E)dV,, dQ, dk, ’ 
Vp 
which is valid for arbitrary initial distribution, can be easily calculated. 


The neutron importance, on the other hand, introduced by USSACHOFF (°) 
is given by 


(5.10) INC Te ad =| J fne ro ob By sey aie 


Vp 0 


According to the instruction under (5.2) the equation satisfied by the ge- 
nerating function G(t, ro, 2), Ly); q, g') can easily be obtained. After simple 
calculation we find 


(5:11) G=G,+G,+ G,+ Go, 
where 
(5.12) Got, Fo, Lo, Hy; q; q)= 


= Fit, r,, 2, #,) {A(V, r) La, Sy, eee AY, Foa) e ets} 


t=1 


and 


t 


(5.13) Ga(t, ro, Lo, Ho; q, q’) = Uo Qi pot), EJF, ro, 2, E) dt’; 
0 
further 


to 
(IO) vol {fee Fo Stes ta) Opal )ietele 
00 


‘ws [polt'), Ho, 2; E, Q'E[t—t’, polt'), L', E'; q, qg']at db’ dQ’ 
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and finally 
t 


(5.15) Gls Qo, Berg) = Po] PU Faas BOF Tol Bolt Lee 0) ko, ka) 


ky 
0 


-Gilt_—t, P(t’), QQ, Ho, ko; qa ile st, Poll’), 925, Eo, ki; gi q']}at'. 


The functions G? and G° can be reduced in one or two steps to the function G, 
namely, 


© ve 
6.169) a7 ={| | Bo, Do; PDAL, pull), DB q, a’ 1AE AL 
v0 
and 


(5.17) = {G;( (t—V, pol), 82, Hy; g, g’] + exp [— (¢—-# Wace 


where 


t-t o 

(5.18) Grlt—1, poll), 2, Hos gg |= i) [exe [— t"A,]w; (Ho, 8293 E’, 22')- 
Sans) 
-G[t—t’—t’, po(t’), 2, E'; gq, q'] dt’ dH#' dQ’. 
According to the relations (5.4), (5.5) and (5.11) the equation for the first 
momentum of the distribution function can be written as follows: 
t o 

(5.19) M,(t, To; NPA H,)= M(t, ri, Lo, E) + af | Q*[polt'), EF, To, o, E,): 


-w* [eo(t'), Ho, Qo; HB’, QW M[t—t’, colt), Q', H'] dt! dH’ dQ" 


t 
1 
+05] OFLeu( E,\F(t, ro; Q,, Ey) dI vAlt—#, Po(t’), o, E,.Jdt', 
i=0 


where 


(5.19) M(t, ro, Lo, Ho) = Apre )A(Qo, Q,)A(E, Eo) F(t, ro, Lo Zo) 


Po{t) 


and 
(5.20) Alt-#, pot’), Lo; Bo] = i | w8( Ey, Qq3 PD Mld pol’), 2, BAH dQ", 
0 


whereas for 1=1, 2,..., J 


t—t' © 


(5.20') Aift—t’, polt), Lo, Lo] ike exp [— t"4.)w;(Eo, Qo; BH’, &2')- 


-M[t-t-#, p(t’), 2', E]diaz'aQ' . 
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For the second momentum of the distribution function the equation 


(5.21) M,(t, ro, Q, Hy) = Mot, Ty) 829, Hy) + m J{ [OFLeu, E,JF(t, ro; Lo, E): 
00 


-w*[eo(t’), Ho, Lo; E, Q' | M[t—t’, eo(t’), Qo, E'] dt’ d#' dQ’ + 


U 


t 
at to| OF Lou EF, To) Do, Ey) DI v.B.{t — l, Poll), do, E] dt' | 


i=0 


t 


+ 05/0? [eo(t’), EF, ro, Q,, E), Poll), Lo, E,\ dt’ , 


0 


can be derived, where B, (î=0,1,...,) is identical to the expression A; 
(¢=0,1,...,1) of Eq. (5.20) only M, has to be replaced by M,, while 
l 


(5.22) CE (AO ATE 


i=0 (VIE 


Eqs. (5.21) and (5.19) are nearly of the same construction. The only difference 
is due to the last term representing inhomogeneity. Eqs. (5.19) and (5.21) 
are often favourable in determining M, and M,, being suitable for iteration. 

Finally, let us derive the transport equations from Eqs. (5.19) and (5.21). 
Making use of the fact that F(t, ro, 2), Ho) = F(t’, ro, 20, Ey) F(t" pol), Lo, Lo) 
and carrying out some simple transformations, we get from (5.19) 


OM, (t, ro, Lo; Bo) 
ot 


(5.23) iy Qo (To, E)voMi(t, ro, o, E) = 


t 
= vd grad M (t, ro, Lo, Lo) + Die (has Em} [od ro fy, So; E', 92’) S 
0 
l 


~My (t, ro, 2, EB") dH dQ! + OO Eo) » V:Ai;(t, To, Lo, Hy) - 


i=0 
With a similar procedure we obtain from (5.21) the integro-differential equation 


OM, (t (940% 
(sas) rt rar Ser Te) OF (ry, Da), Malt, ro, Lo, Hy) = 


= vi: grad M,(t, ro, do, E) + OO or 5 | [we (ro, Ey, Lo; B',Q') 
0 


‘Ms(t, ro, Q', E')dE'dQ'+ VQ? (To, E) 
=F 07 (To, E) C(t, To; Lo, Ey). 


v;Bi(t,r0, Lo, E) + 


a 


iM~ 
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With the help of Eqs. (5.19) and (5.21) the initial condition of (5.23) and 
(5.24) can be easily written: 


(5.25) lim My(t, ro, Lo, Ey) = lim M(t, ro, 2, Hy) = lim Molt, ro, 2, Ho) = 


t—+0 t—0 


at AV, P,)A(2q; 2) A(#, E) < 


So as to formulate the boundary conditions we may proceed as follows. 
If P, is inside the region V, the Hqs. (5.23) and (5.24) are valid, only Q; 
(i=a,s,f,0) has to be replaced by Q; (iî=a,s,f, 0). If PF, is outside the 
region V,, on the other hand, according to the definition of the cross-section 
Q* (i=a,s,f,0) the equation 


0g(t, To) o, E) 
ot 


= vio grad p(t, To, Lo, Lo}, 


holds both for M, and M,, as long as the neutron does not reach the boundary 
surface of the region V,. The initial condition for 
Eq. (5.26) is equivalent to that for Eq. (5.25). 
Concerning the solution of Eq. (5.26) the following 
should be noted. Suppose that from the end-point 
P,, of the vector ro the region V, is viewed, 
inside the cone 2, (Fig. 1). If Q, is outside 
this cone, the solution of (5.26) is given by 


(5.27) pt, ro, Lor Boi Vi» La, E) = Mt, I o, Qa; H) = 


= M,(t, r,, Boi Ip; Qo, E) = A(V,; P.. 0) Aa; Q,) ACH, E) - 


If Q, is inside the cone Q, defined above, the following may be said. Let ty 

be the period during which the neutron 
cosd=(20 25) COS Hy = (2p Qs) of energy 4, starting from the point P,, 
Ao 7 in the direction Q, arrives at a point 
defined by the vector r, on the bound- 
ary of the region V, (Fig. 2). It is evi- 
dent that in the time interval 0 <t< to 
the solution of Eq. (5.26) is identical 
with the expression (5.27), but for a 
t>t, already (5.23) and (5.24) hold. 
Denote by @, the normal belonging 
to the end-point of the vector r, oriented outwards. Finally the boundary 
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conditions may be written as follows 


A(V,, Pr) A(Qq, 2) ACF, 25) 
for) Oc (Qi 


lim M,(t,r,,2,, Ey; V., 29, E) = 
rors M(t, MU, Quo Va La; E) 

for 1 (Qi 2920 
(5.28) | and 


AVERE 


Pol8) 


)A(Qp, 2) A(E, Ey) 
for. 0h 


lim M,(t, To) Lo; Ko; Me, Qa, E) Fi 


| rors 


This general theoretical treatment of stochastic processes taking place in 
reactors will be followed in a subsequent paper by some examples showing 
the application of the theory. 


M,(t, r,, Lo, Boi Ven La: E) 
for —1<(Q,:82,) =O. 
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1. — Introduction. 


In recent years the problem of construction of powerful electronic syn- 
chrotrons with energies of the order 1GeV and more has attracted great 
attention. 

Machines of this type are destined for the production of and experiments 
with hard Bremsstrahlung y-rays and for experimental work with high energy 
electron beams. The use of y-rays, an electromagnetic radiation which has 
been extensively studied, greatly facilitates the performance and interpretation 
of experiments on elementary particles. 

Great electronic synchrotrons are essentially characterized by the presence 
of intensive relativistic electromagnetic radiation of electrons in the magnetic 
field of the accelerator. The power of this radiation rapidly increases with the 
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increase of the electron energy / 


626 
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where r, is radius of curvature of the orbit, and 


Wolk ee 
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Consequently, already at not very high electron energies the major part of 
the power of the accelerating system is expended on compensating for radiation 
vosses. Thus at 109 eV the losses amount to ~ 20keV/r. An increase of 
energy merely by ‘a factor of 10 (to ~10GeV) results in losses up to 
20 +30 MeV/r. 

Another essential feature inherent to the radiation of relativistic particles 
in general is its sharp directivity: practically all the intensity is concentrated 
within a narrow cone of angle y+ about the direction of the instantaneous 
velocity of the electron. Calculations showed that within y~? all radiation 
may be assumed to be tangential to the trajectories of the electrons. 

The radiation spectrum is very rich in higher harmonics of the funda- 
mental frequency ¢c/R. In contrast to the non-relativistic case the maximum 
of the radiation spectrum has the frequency 


(1.2) Da = = È y, 
which at high energies (> 1 GeV) corresponds to the region of soft X-rays. 
The foregoing considerations as well as formulas (1.1) and (1.2) are de- 
rived from the classical theory of radiation. To determine their limits of 
applicability the following simple criterion may be used ('): the quantum cor- 
rections for the fundamental characteristics of the radiation are small if its 
wave-length considerably exceeds the de Broglie wave of the electron, or, 
which is the same, if the energy of the electron is much greater than the energy 
quantum £ >hiw,. Formally these conditions may be resumed as follows: 


me] 

(1.3) E<E,=me al 
As seen from (1.3) the movement of electrons up to the highest electron 
energies possible in a synchrotron may be regarded as classical, while in ge- 


(1) V. VLADIMIRSKIJ: Zu. Éksper. Teor. Fiz., 18, 392 (1948). 
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neral the effect of quantum (statistical) fluctuation of the radiation itself 
cannot be neglected (2). The fundamental character of this effect was first 
pointed out by SoKkoLov and TERNOV (*). 


2. — Motion of electrons in an axially-symmetrical accelerator in presence of 
radiation. 


Let us first consider the movement of a relativistic electron in an axially- 
symmetrical magnetic field under action of an external accelerating field taking 
into account the radiation. 

In the present paragraph we shall discuss the results obtained within the 
frame-work of the classical theory, basing on the relativistic equation for the 
movement of a point electron (?). 


(2.1) Sie 
ade ds 


2e2 (d2U; a? U,, 
Bq edge oh ase |? 


é 
=— FU; + 
(6) 


where U, are the components of the 4-dimensional velocity and F;, the com- 
ponents of the field tensor. The quantity in brackets is the « self-acting » force, 
due to radiation. The Z-axis is directed along the magnetic field. Since we 
consider the movement about the equilibrium orbit with a constant radius £, 
we must introduce the electric field €, which compensates for energy losses. 
(We shall consider the familiar case of € having only the azimuthal com- 
ponent €). The field components have the following form: 


— Jit 
| Cae Afr) TH.) Tee ; 
(2.2) J x R 
deli ; TR 
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In the case of a synchrotron € may be represented as a travelling wave (see 
for instance ref (°)): 


2) A. KoLomenskiJ and A. LEBEDEV: Zu. Eksper. Teor. Fiz., 30, 205, 1161 
(1956); Symposium CERN, Geneva (1956), p. 477. 
3) M. SANDS: Phys. Rev., 97, 470 (1955). 
(4) A. SogoLov and I. T'ERNOV: Zu. Eksper. Teor. Fiz., 24, 249 (1953); 28, 431 
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(1955); A. SogoLov: Suppl. Nuovo Cimento, 3, 743 (1956). 
5) L. LANDAVU and E. Lirs1c: Theory of Field (Moscow, 1948). 
6) D. Boum and L. Foupy: Phys. Rev., 70, 249 (1946). 
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where V, is the total amplitude of the voltage at the accelerating gaps and 
g the phase of the particle relative to voltage. 

From (2.1)-(2.3) we obtain the linearized equations of motion in the fol- 
lowing form: 


È CA ee Vn sing, Erase pE- E 

(2.5) Ti (H — E.) 300 W.|(1-2n) R + 2 E, : 
iL Gl W, e? 

9 DEIR A I | 2 0 

(2.6) E di [E] 4 E, a+ pa 3 


(where the index s denotes the equilibrium values). In addition, in the right- 
hand part of (2.4)-(2.6) appear quantities of the type (W,/H,)(R?/c?y?)o ete., des- 
cribing the low « proper » radiation damping of the betatron oscillations. Since 
these terms are of the order y~* relative to the terms describing the funda- 
mental radiation «friction », they can be ignored. To equation (2.4)-(2.6) 
must be added the linear relation which bounds the values 7 =qg—g, and 
o=r— È 


a È (6) 
(2.7) ia dro 


where gq is the number of r.f. harmonics. From the system (2.4)-(2.6) we derive 
the equation of radial motion 
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where 2 (within small corrections) is the frequency of the synchrotron oscil- 
lations 


qereV 
= 2nk =) 


V= Vp sing, . 


21. Vertical oscillations. — Consider first the more simple equation of 
z-oscillations (2.6). It differs from the familiar equation of betatron oscil- 
lations (valid in absence of radiation) by the presence of the term (W,/E,)é 
which corresponds to the vertical component of the radiation reaction force. 

In absence of compensation for radiation losses H, = — W, and the oscil- 
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lations amplitude a, remains constant. But if, as is always the case in practice, 
the accelerating system not only compensates for radiation losses, but imparts 
to the electron the energy EH, = eH(R)R, the amplitude a, according to (2.6) 
varies as 


t 
AVI 
DIO Spad = dz). 
( ) a, ~ E; exp| [| 


DI 
4 
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For the evaluation of the corresponding relaxation time we obtain the ex- 
pression 


2-107? 
2.44 Idi 
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where T,, is the time of acceleration of the particles up to the highest maximum 
energies. 

It must be noted that even if the energy is imparted to the particle not 
continuously, but by discrete portions, in the accelerating gaps, radiation 
damping would practically be still described by the same formula (2.10). 


3:93. Radial oscillations. - Now consider the radial oscillations 0, described 
by the equation (2.8). It is easily seen that in the zero approximation, in 
which terms of the order W,/H, #./E,, etc., are neglected, the solution of 
this expression is a superposition of rapid 0; and slow 0, (betatron and 
synchrotron) oscillations. 

The solution of (2.8) which is a first approximation equation, will be sought 


in the form: 
t 


invi n ) t+ Qpelt) EXP i | Qar + c.c. , 


(2212) 0 = 4»(t) exp 


0 


where 4, 4. are slow functions of time, the first deviations of which are to 
be considered as being of the same order of smallness as W,/E;, H,/H., hereby 
W,/E, < Q:(c/R)V1—n. This last assumption means that the radiation damp- 
ing during a period of betatron and synchrotron oscillation is small, which is 
true up to H~10" eV. For the law of variation of the betatron oscillations 
amplitude the following results are obtained 


(2,13) ao ~ E;* exp 
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while for synchrotron oscillations the expression obtained 


(2.14) ann VB st exp|— 5 [ ae, 
; ff ‘ 2 1—n', 
0 

coincides with the result of ref (**). 

As seen from (2.13) and (2.14) the oscillations 0, are subject to a powerful 
radiational friction. The essential difference between them and the above 
mentioned case of ¢-oscillations (see (2.10)) consists in the dependence of the 
sign and magnitude of friction on the value of the field index n. 

If the condition 


(2.15) Oni 


characteristic for weak focusing is maintained, the damping of o, occurs. If 
this condition does not hold the antidamping occurs (*). 
The damping of o, takes place in the intervals 


(2.16) ZE po 
The corresponding relaxation time in this cases has the same order of mag- 
nitude as for 2-oscillations. 

The radiation friction for radial betatron oscillation (see (2.13)) is a result 
of the superposition of several effects differing in character: a) the effect of 
the radial component of the radiation reaction force in presence of compen- 
sation of radiation losses leads to damping similar to the damping of vertical 
oscillations; 6) the dependence of the power of radiation losses W (see (1.1)) 
on the value of H in a given point, W~ H?(r), leads to a perturbation of É 
equal to AH = 2nW,(oR/)(see (2.5)); c) the effect linked with the fact that 
for a particle deviating by the amount o the length of the orbit is increased 
as compared with the equilibrium orbit, while the power of compensation 
of the losses is not directly dependent on 0, leads to AH = — W,(0/R). 

Owing to effects 5) and c) the radiation losses and their compensation de- 
pend on the phase of betatron oscillations and are subject to rapid variation 
occurring at the frequency of this oscillations. 

The simultaneous effect of b) and c) may be qualitatively explained if o, 
is considered as the oscillation of a damped pendulum under action of a quasi- 
elastic force f, ”—(1-— n)(0/R). According to the above, the equilibrium 
point of the pendulum, the position of which is characterized by AE= E — E,, 
itself oscillates and is shifted with a velocity AE ~ (2n —1)W,(0/R). This 


(*) Antidamping of this type and methods of suppression are discussed in Sect. 4. 
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may lead either to damping if A# and /, have opposite signs, or to antidamp- 
ing if this condition does not hold. 
In conjunction with the effect a), it leads to damping conditions (2.15). 


8. — Consideration of radiation fluctuations during the motion of electrons in 
an axially-symmetric accelerator. 


As noted in Sect. 1, to take into account quantum effects in the pre- 
sent case means that on the basis of the classical theory we must additionally 
take into account the corrections due to statistical fluctuation of radiation. 
We make now the following assumptions, valid for H< E,. 


a) The emission of individual photons on the equilibrium orbit is sta- 
tistically independent, i.e. obeys the Poisson law. The correlation correction 
is not essential as it can be easily seen from (1.1), (1.2) according to which 
hw, E<1. 


b) The emission of photons occurs instantaneously. The uncertainty of 
the moment of emission is of the order of dt ~ 1/a, ~ R/cy?, i.e. is negligibly 
small as compared with any of the characteristic time scales of our problem. 


c) Spectrum and intensity are derived on the basis of the classical theory 
and are valid up to terms of the order [E,£,]?. 


The validity of this treatment in the region E < È,, i.e. at hon< E follows 
from the general assumptions of the radiation theory of quantum electro- 
dynamics (see for instance ref. (?)). 

As an illustration it can be noted that in a 5 GeV accelerator the electron 
moving at the highest energy emits — 100 photons/r with a wave-length 
in 21À (radius of the orbit ~ 15 m), the ratio E/E, being ~ 1078. 

It follows from the foregoing that in order to account for the statistical 
character of radiation losses, the equation (2.5) for È must be completed by 
introducing 
(Sal) w(t) = > e,d(t—t,), 


a 


where ¢,—fiw, is the energy of the 1-th photon, and 6 the delta-fvnction. 
By taking into account w(i) the system of equations (2.4), (2.6) is turned into 
an inhomogeneous system and accordingly (2.8) becomes an inhomogeneous 
equation with the rigth-hand part (c?/R)(w/H,). Further changes in (2.8) con- 
nected with the consideration of w(t) lead to corrections of higher orders of 
smallness since the fluctuations of losses may be regarded as perturbations. 


(7) A. AntesER and V. BERESTETSKIJ: Quantum Blectrodynamies (Moscow, 1955). 
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The solution of an inhomogeneous system corresponding to zero initial 
conditions has the following form: 


t a t 
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where 2D is the Wronskian of the homogeneous equation (2.8), and D,, D, 
are the corresponding algebraic supplements. 

Let us determine by means of (3.1) the value of the mean square ampli- 
tude (<a> = 20>. In order to do this we average over the photon spectrum 
and over the phase of betatron and synchrotron oscillations. Since the oscil- 
lation frequencies are widely differing, the averaging over each is performed 
independently, using the expression for the spectral density of radiation of 
the relativistic electron moving along a circumference (see ref. (?:*)). 


[co] 


fa 2 a 
(3.3) U(y) dy = SE = vty dy Jr da . 


From (3.1), (3.2) (3.3) we obtain 
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where A is the Compton wave-length, and 7, = e? me? the « classical radius » 
of the electron. 

Now consider briefly the qualitative physical picture of the discussed effect. 
At the moment of the emission of the photon e;, the radius of the closed orbit 
is decreased by the amount do. = — (Re;/(1 — n)E,) with the following change 
of 0, by the amount 50,=— do.. It is easily seen that during one half 
period these jumps lead to a decrease of the oscillation amplitude and during 
the other half period, to its increase. It is well known however that, owing 
to the statistical character of the process, the amplitude in absence of damp- 


(8) D. IvANENKO and A. SogoLov: Classical Theory of Fields (Moscow, 1948). 
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ing increases, in the mean, as the square root of time (this is true for 
radial betatron as well as for radial synchrotron oscillations). 

The process outlined above might be prolonged at infinitum. The taking 
into account of the familiar adiabatie damping (~ E;?* for 0, and ~ E> 193 
for 0.) does not substantially alter the picture, and the characteristic stochastic 
rising of oscillations ~ Vt is virtually maintained. This picture however is 
fundamentally changed if we take into account the radiation friction (2.13) 
and (2.14) mentioned above, which limits this rising to a damping time equal 
to & ((1— n)/n)(E./W.) for 0, and a ((1—n)/(8 — 4n))(E./W.) for o,. This 
problem will be considered more in detail in the section devoted to arbitrary 
magnetic systems. 


31. Radiation fluctuations and vertical oscillations. — The mechanism of 
excitation of ¢-oscillations by fluctuations differs from that of o-oscillations 
outlined above. It is to be remembered that the photon is emitted in a di- 
rection deviating, in the mean, by a small angle from the direction of the 
electron velocity. At each act of photon emission a recoil pulse is imparted 
to the electron, which leads in the average to the excitation of oscillations 
similarly obeying the law ~ Vt. In the case of g-oscillations this effect leads 
to only small corrections to the expressions (3.4), (3.5), while in the case of 
z-oscillations it resumes the total effect of fluctuations. 

It is to be noted, however, that radiation damping (2.10), limiting the du- 
ration of rising by the time of relaxation, impedes any significant increase 
of the amplitude a. Performing a series of computations according to the 
scheme outlined above and averaging over the photon spectrum, the angles 
of emission and the phases of vertical betatron oscillations, we finally obtain 


the following expression for a. 
t t 


rks ie [ at 
(3.6) <a, > = 924 v3 Rny, exp |— | E. ac] y da . 
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4. — Betatron oscillations in magnetic periodic systems in presence of radiation 
and methods for the damping of oscillations. 


In the present Section we consider some peculiar features of radiation in 
magnetic periodic systems. 

In the general case, when the curvature x of the equilibrium orbit and 
the field index n are functions of the azimuth, the linearized equations of 
motion for 0, 2, n are easily obtained (ref. (°)). As to their form these equations 


(9) A. KOLOMENSKIJ: Symposium CERN, Geneva (1956); Zu. Telm. Piz.: 26, 1969 
(1956). 
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coincide with (2.4)-(2.6), but differ from the latter in that these parameters 
as well as W,(x) are functions of the generalized azimuth. 

As usual the solution for 0 will be sought in the form of superposition of 
a periodic solution corresponding to a closed orbit and betatron oscillations 
near this orbit 


(4.1) 0 =o. + ¢ exp [ium]f(9) + c* exp [— ium] f*(9) , 


where f(#) is the function of the azimuth depending on the magnetic structure 
and m the number of the successive elements of periodicity. 

We shall limit our discussion to the more characteristic case of a system 
consisting of different sectors, with an axially-symmetric field in each of them 
(strong-focusing synchrotron, race-track). 

From the results of Sect. 2, obtained for an axially-symmetric field, it follows 
that in each sector of a strong-focusing accelerator with |n|>1 we will not 
have radiation damping as in (2.13), but radiation antidamping effect with 
approximately the same decrement. The physical meaning of this phenomenon 
has been elucidated in Sect. 2. A simple consideration shows that the same 
may be extended to accelerators in general. In the discussed periodic systems 
the amplitude of betatron oscillations varies according to the law 


da 


AE, 
Soli (aces 1) ion ; 


where « is the momentum compaction factor. 

According to evaluations in a strong-focusing accelerator the increase of 
amplitude by a factor e must take approximately 2-10-?/H (GeV) s, which 
gives to the effect of radiation excitation great practical importance. 

Owing to this effect electron acceleration in strong-focusing accelerators 
up to energies of 5--10 GeV and even less becomes very difficult and even im- 
possible unless some special provisions are taken. 


41. Methods of overcoming radiation antidamping. — We have devised and 
studied several methods of overcoming the radiation antidamping. 


A) One method of overcoming radiation antidamping of radial oscil- 
lations consists in the use of artificial strong coupling between radial and ver- 
tical oscillations. The latter in all systems are always subject to radiation 
damping according to the law (2.10), independently of n. 

The working point must be located on the line of differential resonacen 


(ref. (‘)) 
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while the coupling is performed by rotating the sectors around the azimuthal 
direction or by introducing special magnetic lenses into the straight sections. 
Computation gives for radiation decrements of coupled radial and vertical 
oscillations the following expression 


(4.4) Cros = — 4% @ (2— a), 


where x is the value characterizing the oscillations « bounding » which is equal 
to 1 in absence of coupling or at great enough distance from the resonance. 

In the opposite case x +0 and we obtain for the decrements the ex- 
pression 


We 
Li 


È 1 
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It will be seen that practically neither 0, nor 2-oscillations experience any 
damping or antidamping, while at the same time the amplitude of vertical 
oscillations equals that of radial oscillations which rapidly rise in absence 
of damping owing to fluctuations. 


B) A damping of radial oscillations might be achieved also by creating 
the accelerating voltage amplitude dependence on the radius Y= Vo(7). This 
method is not directly related to radiation, but may be used with good results 
in powerful electronic accelerators in which the accelerating voltage has to 
be very high owing to radiation. 

It is physically obvious that this dependance will affect both betatron 
and synchrotron oscillations. Successive computations give the following 
results: 


a) Synchrotron oscillations are subject to additional damping with the 


decrement 
ra e {yn,) HEV, cos po _ 1 oie n) 34 IR @Va 
(4.6) hd Ae ED, DI lls dee ; ae 


where for definition of y see (5.1) and L, is the equilibrium orbit’s length. 
Since in strong-focusing accelerators the value of « is small («~5/|n|) the 
effect of damping becomes felt only at great n,. 


b) The damping of betatron oscillations, owing to the dependence V, (7) 


is characterized by the decrement 


1 xeVo COS PC de; 1 xW, 
CIR e IE, 
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As seen from (4.6)-(4.7) the damping of synchrotron oscillations is followed 
by equivalent antidamping of betatron oscillations and vice versa. 

Nevertheless this last method may prove useful for the damping of syn- 
chrotron oscillations at the expense of betatron oscillations, since the former 
may prove more dangerous and vice versa. 

A better result may be achieved in a strong-focusing accelerator by using 
a peculiar damping system by means of a variable azimuth voltage gradient, 
for instance according to following scheme 


OV OV | 
3 SI : oi) 
cr cr 


| | t 
\n<0|n>0 | i S| ae 0) | n<0/n>0|. 


As shown by computations, this may lead to the damping of betatron oscil- 
lations, while the synchrotron oscillations remain practically unaffected. An 
alternative, less convenient, method consists in making the sectors n > 0 
and n<0 of different length. 


5. — Effect of radiation fluctuations on betatron and synchrotron oscillations 
in magnetic periodic structures. 


The effect of radiation fluctuations for this case can be taken into account in 
the same manner as for the case of an axial-symmetrical accelerator (see Sect. 3) 
by solving the corresponding Hill’s inhomogeneous equations and averaging. 
In order to make this picture physically more clear, we give here another 
deduction of the expressions for the mean square values < don) and Ka), 
<a> separately. 


51. Betatron oscillations. — Suppose that on a certain azimuth the electron 
emits a photon with energy e = haw. This gives rise to a jump-like change 
of the localization of the closed orbit and its slope in the said point, with a 
simultaneous change of the analogous characteristics of the betatron oscil- 
lations 


| so. = —3e, =—— y(0), 
(5.1) ic 
Sole Ope = y'() 


where y(2) is the known periodic function depending on the type of the mag- 
netic system and having the same period (see ref. (ye 
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Using (5.1) and the familiar normalization of the Wronskian W (ff*)= — 24 


we obtain 


(5.2) Kan) = |f 
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The total change of the value {aw consists of two parts: the systematic changes 
of (4.2) and of the stochastic part (5.2). 

Taking this into account and using (5.2) we obtain, according to a scheme 
similar to that used in Sect. 3, an expression for the determination of <a?,» 


oie 5D <nWw ,/(1 — n)> E, ne 
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where F is a factor depending on the parameters and configuration of the 
system and equal to 


— 
DI 
ra 
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where x is the curvature of the equilibrium orbit. The averaging in (5.4) is 
performed over the element of periodicity, including the sections free of field, 
to which corresponds x = 0. 

In the case of an axially-symmetric field “= 1/R(1— n)? and we obtain 
again the formulas of Sect. 2, 2. 

By integrating (5.3) we obtain the expression: 


t 


dx 
(5.5) Cag»? = CM) EXP | — ((a—1) W.+ El + 
n s 
t t 
3 Wee rae 
ce 24/3 nAel | exp | — fee —1)W.+ 5))) zl da; 
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which is a generalization of (3.4). 
If, owing to the small value of energy or of the acceleration time, the radiat- 
ion friction has no means of manifesting itself, (5.5) leads to: 


t 


5 LAI 5B rAceF 
(5.6) Qob)in = (Ceo 0 Fg, 24.4/3 y yo da 


The first quantity in (5.6) corresponds to the « classical » part of betatron 
oscillations, while the second expresses its « fluctuation » part. Their behaviour, 
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in general, essentially dependent on the sign and magnitude of the para- 
meter 


If £,> 0 both parts of (5.3) are subject to the antidamping discussed above, 
the nature of which, as was demonstrat- 
ed, is purely classical. 

If a damping system is used we 
obtain from (5.5) an established value 
of the amplitude. For example 

PIRA 
(5.7) E Calda = se LE va 
We give below a table showing the va- 
lues of F as a function of the selection 
of the working point in a strong-focusing 
magnet. 


For the sake of simplicity consider 

the case when the gradient in the sectors 

SY i of different types has the same abso- 
3 lute value, and the sectors and straight 


DE sections are of the same size, 


la = 1 15 Ne eh ee enti 
0.4 x 
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In case of a weak-focusing accelerator 
with straight section (race-track) the 
factor F differs from 1/R(1— n)? by 
a few percents only. Fig. 1 shows a 
typical dependence of the amplitude 
of radial betatron oscillations on the 
magnitude of £ = t/T,, where T,, is 
the time of acceleration; the approxim- 
ating curves I and II correspond to 
(5.6) and (5.7). 

Note that for vertical oscillation, as a result of a similar computation, we 
obtain: 


0.2 


0.1 
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The established values of their amplitude are independent of energy and for 
all practically interesting cases are negligibly small. 


TABLE OF THE VALUES OF FACTOR I, 


uu t/Rv=0 Ry == 0,3 1/Ry = 0,5 


0.78 850 430 | 330 
1.30 530 220 180 
1.65 510 210 150 
1.94 620 310 | 110 


52. Synchrotron oscillations. — Consider separately the slow synchrotron 
oscillation of radius 0,, phase n. and energy AH,. For simplicity we shall 
discuss only the case when the curvature of the orbit in all magnetic sectors 
has the same value x = 1/R. According to phase determination 


(5.8) Ne = G(Ors dp Or) 3 


where ©, is the frequency of revolution, we obtain the following useful re- 
lation: 


5 quae AE, 
(39) Ya = R20? SYP? x i? 
where <p>, denotes an averaging over the sector with magnetic field. The 
radiation power of the particle W(r.) moving with energy EF, = E; + AH, 
in a magnetic sector along a closed orbit of radius 7, = R+o0,(0), according to 
(1.1), has the form: 
ap eee 


(5.10) Wir.) = W. {a +2 | a) eel 


where W, denotes the power of the radiation losses during the motion along the 
circular section of the equilibrium orbit with equilibrium energy W,=W(f, £.). 
Now using the law of conservation of energy for the non-equilibrium particle 


and taking into account the presence of radiation fluctuations, the phase equa- 
tion takes the form 


d 


mc 
(211) di 


d ce il 
NES a Az.) — 3 AH.) — = 2 €,(@; — @;) . 


ae 


The quantities ((d/dt) AE.),, and ((d/dt) AE), denote the variation of the value 
of AE, per unit of time due to the accelerating field, and the radiation losses 
respectively. The last quantity in the right hand side describes the radiation 
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fluctuations characterized by the following values: e; the energy of the emitted 
photons of the j-th harmonic, and a, the number of photons emitted per unit 
of time by an electron moving in a magnetic field; factor o = 1--(l/Ry) takes 
into account the presence of sections free of field and radiation. 

Using the expressions derived above and (5.11) we obtain the linearized 
phase equation for slow synchrotron oscillations 


Sta cy) = 
(oily) | (= | 3 a y+ 27 = ata, 


where Q is the frequency of phase oscillations 


Bali: a _ AVS 
o) “ein 


and % denotes the parameter 
(5.14) Q=2+ 


From this general equation may be derived the result obtained by SANDS for 
the case of a race-track (ref. (*)) by correcting an error in (*), i.e. by mul- 
tiplying the right-hand part of the equation by a factor 1/o which takes into 
account the decrease of fluctuation effects due to the presence of straight 
sections. 


Solving (5.12) we obtain the mean square value Kai > 


t t 
E: 5b geAccp>y,V—1 05 ( OW, 
5.15 UE SEA, sul xa “dé 505 da 
(5.15) <a; 1273 Riot Jesp lease de ye VY tala. 
0 x 
t 
Hence at | (QW,./oF,) dé >1 the established value may be easily obtained 
0 
st VA Lay 
(5.16) ida = O 
326 IAT DERE 


In case of non-established oscillations, when radiation damping is still weak, 
t 
which corresponds to | (QW./oH#,)d&<1, we obtain from (5.15) the following 
0 
expression 
t 
(6.17) Cai.) = 55 4 roAe? pnd cts ps ae [yesves da 
ee 244/83 Roty Von , 


0 


For the case of weak-focusing 


(5.18) Ge eee 
il = 4) 
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For strong-focusing, according to (5.14), we have for the centre of the stableregion 


(A) er AE 


cl 


The total behaviour of the function <a?,)*(É) is shown in. Fig. 2 for both weak 
and strong-focusing. 


It must be noted that the amplitude of fluctuation excited phase oscil- 
lations is rather large and slowly decreases in the region of high energies 
(asy-*). Therefore the azimu- 


. 2 
thal size of the accelerated part- Canc) 
rad 


% 


icle beam is great enough to 
practically inhibit all coherent 
radiation (31°). As seen from 
(5.15), (5.17) (see also Fig. 2), in 


2 \t 
ne/ 


increases. The maximum value, 
however, is not necessarily rea- 
ched at the end of acceleration 
but might be reached long before 
that, depending on the accelera- 
tion time and other parameters. 
In the latter case (a?,>* pas- 
ses through a maximum with a 
further decrease according to 
(5.16). If the maximum value 
of <a? >? is comparable with the 
size of the separatrix, it will 
(unless some special measures 
are taken) lead to the loss of a 
fraction of accelerated particles. 
We give below the formulae 
for the mean square of radial 
synchrotron oscillations. Note, 
that the maximum deviation 
o. which plays the role of the 
amplitude a,,, is equal to 


the beginning the value <a 


nel 


) 0 0.2 04 0.6 0.8 1.0 
ax . = . - . - . - . = . = 0 
(5.20) a : a (AE je ) Emax=56eV ; /n/=100 ; p=, 5 R=15m; 0=13:; G=100: Ps= 60 


oe E, Ì 
Inita D2. 
where y,,, is the maximum va- 
lue of (9). Using the relation between 7, and g, (see (5.9)) and the expression 


(10) J. Nopvick and D. Saxon: Phys. Rev., 96, 180 (1954). 
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(5.13), from (5.16), (5.17) and taking into account adiabatic damping, we obtain 


ee 


(5 21) al i DV3 Ato Wax? 
O.41 c/est — j / 
£ a 48 Qo 
t 

~_— 2 T 5 fe 

On valenze VO IRE 
(5 22) EI a o Panna eel = 0:5 da Ss 
O. dA oe/ in 924 v3 R*o%y15 


0 


The absolute value of the excited radial synchrotron oscillations is very small. 
Thus from (5.21) which gives the highest possible value at a given fixed E, 
and we obtain the approximate evaluation 


E? (GeV 
(5.23) (pi — = om, 


which at |n|= 100, even for H = 10 GeV, gives only a few mm. 

In resuming it must be noted that the outlined scheme may be extended 
with small modifications to the computation of the influence of a variety of 
perturbing effects, involving statistical and fluctuation phenomena (brems- 
strahlung on the residual gaz nuclei, Coulomb interactions) as well as to the 
computation of corrections due to discrete feeding of energy into the acce- 
lerating gaps (« slit» oscillations). 
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Contents. — 1. Introduction. — 2. Rigorous solution of equations of dif- 
fusion type. — 3. On the Feynman’s « integral-over-all-paths ». — 4. Tran- 
sition to continuous #-variable. — 5. The usual perturbation calculus. 


1. — Introduction. 


We begin with the rigorous solution of some equations of diffusion type. 
Further we investigate Feynman’s «integral over-all-paths » and use our pre- 
vious results to modify the computation technique of this method. One may 
say roughly that the relation between our method of summation and that 
used by FEYNMAN (+) is similar to the relation between the Lebesgueian and 
Riemannian integral. Some considerations, especially of combinatorial and 
probabilistic character, in the second part of the note will be only outlined. 


More exhaustive treatment will be published soon (?). 


2. — Rigorous solution of equations of diffusion type. 


21. — Let us call an equation of the form 


(1) op), do; t) A 


ot At->0 


(1) R. P. FEYNMAN: Revs. Mod. Phys., 20, 367 (1948). 
(2) J. Lopuszznanski: Acta Phys. Pol. (in press). 


= lim | [20 IVA fost 29) li st 
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an equation of the diffusion type. The unknown function is (i, jo, t) and 
and 2(¢ 
a time interval At; the initial condition for t = 0 is 


k) are given functions; they are elements of the transition matrix for 


me (j — jo) if 7 is continuous 
(2) pI; Jos $= 9) 
ae? if j is discrete. 


To this type of equations belong e.g. the Schrédinger equation, equations 
of stochastic processes as well as equation describing the sum-over-states 
(in these examples, except the last, t has the meaning of time). If the j spectrum 
of states is not continuous we replace the integral on the right-hand side of (1) 
by a sum (or Stieltjes integral). 

For the sake of simplicity we confine ourselves in our considerations to 
the case when the j-spectrum is discrete and the number of levels is bounded, 
say 1. We assume also that the variable t is discrete. Finally we will discuss 
the case of continuous t. Thus we consider the equation 


1 


(1a) WI, To; n) SZ w(", Îo; REL), 


whesekk=smAt and 9; elia 1.2 


22. — To solve (la) we define the generating function for y(j, jo; #) 
(3) Dj, lo; A) = DI Aj, Îo; N) 


n=0 


convergent uniformly in the neighbourhood of 2 = 0 (*). It follows from (1a) 
that @ satisfies the equation 


(4) DG) 07, Sei (j|y) Dy, jo) - 


Equation (4) yields the solution 


DG Aes 4) = Tie 


(*) If we want tcuse a continuous t-parameter, we must replace the power series on 
the right-hand side of (3) by a Laplace transform of y(t). 
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where 
—Azi|1) +1, —A2e(2/1), (400) 
i — A2(1|2), = A2(2\2) 1, AAA 
JA = (2|2) (2|2) 
— Az(1]|1), — A8(2|1), we.) —A2(0/0) +1 


is the determinant of the «stochastic » matrix (we write inverted commas 
because 2(i|k) must not be probabilities) and J,, is the minor of J(/). We 
put this solution into the form 


(5a) PESO 
lode C2(Jo| 4) fi 

with 

(50) O ve n J(A) . 


Tt is obvious from (3) and (5a) that, if we develop Q in a power series in res- 
pect to A 


(6a) Q(a) = > A a(n) 
we get i 

sE co(n) 
6D J, Jo pars 
dI YA 103) Behalf) 


The solution (5) or (6) has an interesting feature: each solution of a problem 
concerning a process of diffusion (Markoff processes, Schrodinger equation, etc). 
may be written as a partial derivative in respect to a transition matrix element, 
say 2(% |x), of a« potential function » w(t). This function is symmetrical in 
2(i|k). For example we are able to write all elements of the S matrix as « gra- 
dient components » of only one such « potential function » in the 2(t|k) - space. 


2°3. — We prove now the relation 
1 i l l 
(7) PE a TI II (i | ky 
J 2 \ i le II I bi! i 
where the summation is submitted to following restrictions 


Dus= +1. 


l l 
=1j=1 


l l 
(7a) DI Mij = > bi = Mi 
j=l cli 


J 
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Let us write the right-hand side of (7) in the form 


0 ipa gibt) 


ML Hi 


Cy C1 
where we sum over u;, for which 
l 
(8a) Fu=n+1 
î=1 


holds, and the integrals are taken along closed paths €; in the $, complex- 
planes, enclosing the origins of co-ordinates resp. Now we are able to per- 
form in (8) the summation over n. We have 


l i HralC;, 
(9) eas > Bla) P Pp U2 GE = 


We introduce new variables 


then we have 


Li ipa 
TL do, = ap dd, 5 
consequently 
Tales CEE dh, Tae A i 
10 (Lp E = © 
ey 2ri ¢p Gp Ie AD a WANT IAT > D It dr Jr 
ci ci ci 0’ 


where the paths C; are similar to C,, but in d,-planes. 
24. — The function 


(11) ~—1=AK(A) 
tends to zero together with 7. Consequently 


1 
QA) = Fn (1 + AKA), 
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may be developed in the neighbourhood of 2 = 0 in a power series in respect 
to A. In this way we get 


n+1 1 m 
(12) CLC (EA ee, Seve too Ml Bae 
m=1 m its jn uy By i=] 
m U 
Vui=n+1 
s=1 j=1 
where 
l un! 
(12a) d ini II a(t | ky“ aR) > oe > Gy. wy 3 
pu k=1t=1 Hix: Hi Hi 


the summation on the left-hand side is submitted to condition (7a). 


3. — On the Feynman’s « integral-over-all-paths ». 


31. — FEYNMAN (reference (!)) proposed a new method to solve the Schrò- 
dinger equation by means of so called « integral-over-all-paths ». 

In this section we will present a modification of this method; we have 
not, however, in view to modify the main idea of this method but only the 
mathematical formalism. 

Let us illustrate firstly the method of Feynman and further our approach 
on the example of a Schrédinger equation in one space dimension (*). 

We ask for the solution 


0 
Y(x, ®; t) = (a exp e n); 
of equation 
NL) he Oy | dn 
do) di dt 2m da | e: 
with the initial condition 
(14) WL, Lo; t=0) = de — x). 


This solution may be written by means of the Feynman’s method 


(15a) p(n, Lo} t) = lim | aes Jan 200 ane, (Tr (spa) exp I (2, ve) ; 


At>0 . 
nAt=t =_—_ 
nel 


(*) We perform our computations in one space-dimension, the generalization to 
more than one dimension seems to be easily done and is outlined in the last but one 


section of this note. 


5 - Supplemento al Nuovo Cimento, 
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where 


(150) L(#@;, 0; 1) = 


and #,=«#. In (15a) 


CETRA: 


exp i 


m \} i | 
(16) (Fan exp 7 L(0;, Wi) At = % fe , 


is the transition probability amplitude from state x,_, to state x, in thetime 
interval At. If we denote, as in Sect. 2, the transition amplitude from state k 
to j in the time interval At by 2(j|k), the solution (15) may be, of course, 
written 


(15c) (x, &a; t) — lim | a fae, os Xp [Ts] 23-2) . 
AES ‘= 
me ln) i 
n=1 


We conclude from (15c) and Sect. 2, that Feynman’s method may be 
applied to each equation of diffusion type; the solution of such equation has 
the form (15c) and may be immediately written down, if we only know the 
transition matrix 2(j|k). 

So far Feynman’s idea. 

In formulae (15a) the (n —1)-fold integration appears; this summation 
procedure was already modified by some authors (see e.g. (*)). We propose 
next another way of computation: we replace the (n —1)-fold integration 
by sums over groups of « similar » paths. What will be understood by « similar 
paths » will be clarified further. 


3'2. Let us demonstrate our attempt on the simplest example, when «& 
has discrete levels and the number of levels is bounded, as in Sect. 2. Although 
this case is meaningless for physical practice (except perhaps the case of spin 
waves in a crystal), nevertheless this example is suitable to present clearly 
the main idea of our speculations, without troublesome computations. This 
case was already investigated and rigorously solved in Sect. 2. As in Sect. 2, 
we assume also the time to be a discrete variable, the shortest interval being At. 

Thus the solution (15c) becomes 


U l n 


(15d) Veda [| 


@=1 @,=lt=1 


(3) B. Davison: Proc. Roy. Soc., A 225, 252 (1954). 
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This solution consists of a sum over « chains » of length n, which all start in 4° 
and end in x. Such a «chain» may be also expressed by 


n 


TI (a: 


î=1 


(17) 


ti) = [[T [] 2G|*)"*, 


j=1lk=1 


where m,, denotes how many times the term 2(j |k) appears in the concerned 
«chain ». We have, of course, 


(18a) 


U l 
I 
J=1k=1 


n 
Each « chain » [[ 2(@;|<;_,) may be written in a unique manner in the form (17) 
i=1 


on the right-hand side, but not vice versa; the number of « chains » which may 
be characterized by the same ensemble of exponents m,, and which begin in 2, 
and end in x we denote by 


(19) <a 


Ma), 


we quote by M the dependence of (19) from all matrix elements m,, 


Ma 5 Miz 4 Mi, 

Mo 5 Mo 5 Mor 
(19a) M = 

Mau 4 Miz 4 Mar 


We mention that if M is given the length n may be immediately fixed by 
means of (18a). 
The necessary condition to form at least one «chain », if M is fixed, is 


(180) > Mn = > Ma; fork! 
and 

(18€) > Mg "i L= d Mp5 b) 

(18d) = 


ie al 
9 


Relations (18b-d) mean simply that in 


> (ea 
4) 


a «chain » the number of transitions 


from a fixed state, say k, must be equal to the number of transitions towards 


this fixed state, excepted « and %. 


If we exclude from considerations the 
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case when M is only diagonal, then (18) isalso a sufficient condition to build a 
«chain » of length n. It follows from (18) at once that 


(20) > > (j — k)my, = © — %, 


consequently the final state x is uniquely determined by M and 4. Thus we 
call some «chains » or paths «similar », if they are characterized by the same 
M and «. 

Using notation (17) and (19) the solution (15d) may be written 


l l 
(15e) pl, co n) = > IT ITG|M)"*<|M]o), 
M j=lkb=1 
the summation is extended only to such M’s for which (18a) holds. Of course, 
the <0|M]|x)?s, for which (18b-d) do not hold, vanish. 

The expression on the right-hand side of (15e) has the form of a multiple 
power series of the variables 2(j7|k); thus the <«|M|a,>’s play the ròle of 
power coefficients. 

It should be emphasized that the nature of the physical problem is entirely 
involved in 2(j|k)’s (e.g. the mode of interaction of the particles of the system 
is described by their potential energy ap- 
pearing only in the exponent of the tran- 
sition probability amplitude, in a similar 
manner as already shown in this section 
for the one dimensional Schròdinger equa- 
tion in «-representation). On the other 
hand the expressions <«a|M]|x,?s do not 
depend at all on the physical nature of 


< nat > the problem; to find them we must only 
Pigi: employ geometrical or combinatorial con- 
siderations. 


The problem of evaluation of <#|Mj|a,> may be called «the problem of 
construction of fences » (see Fig. 1). Indeed: let us imagine, we 
have got a transport of planks of At cm width. Each of the planks 
is painted on the one side blue and on the other red. From among | 
these planks m,, is cut-off slantwise so that, looking at the red side, 
the left border has the h-ight k cm and the right — j cm (see ‘ 
Fig. 2). The number of all planks is n. Imagine further, that we | | 
have two poles in the distance n At cm one from another. One has 
the height x cm, the other — « cm. 

We ask in how many ways one may build a fence between 
ends so that the upper sides of the planks form a continuous, broken line with 
these poles, of the height of the pole tops resp. (see Fig. 1). Each side of the 


Fig. 2. 
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fence is uniform in colour; if we look of the red side we should have the pole 
of x, height on our left-hand side. 


3'3. — We show now how we can evaluate approximately the expression 
<a|M]|a,). We use combinatorial considerations as well as make extensive 
use of the law of large numbers. We omit here, however, these considerations, 
because they are rather long, although they are interesting and instructive. 
These heuristic considerations will be published soon in extenso in a longer 
article (see reference (?)). We quote now only the final result, giving a good 
approximation for large n, 


ITs! 

(21) Cx Maya i 
Pr II TT es! 

k i 

for 4, sqtisfying (7a), and the generating function 
= lo = 
22a asd = = i 
(22a) Bla os 9) = ay LA» 
with 
~ ib fi OU bya 

22 = -—-1). 
(220) COSE i ; | = 


We show, that the functions È satisfies approximately the equation (4); 
deviation appear only for small n. Since 


eo) pie a 
ard: AJ A \2ai J een CAD ati] Re, CAD elio) 


we have 


i dA (DIO! e degen: 
83) AY aWwly)d SR a 


L 53 dA J dita, dé, 
CIT gala (oni) Pe Pd, > easing 


[9 1 | RT | 
È le 71> ali) " G-AD Ail EW] 


The second term in the bracket of (23) may be written 


di i dico , 1 È: 
ea [ea p » PI BAS Ai RX Ala ~ 

1 dé, 5 i Mu — Dl AO 
ae ia Pe PI Ci AD alii LAZ ali lanes RI: 


IR 
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because only the terms belonging to large n are of importance and thus we 
may replace 


da An 
[Zee TO 
by 
1 An 
z {ardor lan 


The first term in the bracket of (23) vanishes for x 44%, because 


"i dl d0.t, 9 i 
(25) pps a oe lea A a — 


i 


"PP Uri) 2° 
7 bo ae, \ LT i AYA) 


È il 
lesa 


as 


is an one-valued function of ¢, in the neighbourhood of ¢,= 0. 

Thus, taking into account (24) and (25), we see from (23) that equation (4) 
is satisfied for x # x. 

For x=% the first term in bracket of (23) yields 


aad 
Ane 


Thus the right-hand side of (23) may be written, taking into account (7), 


if dZ lo ae di 1 
A} A 0(@|a) J Asda 


foo) An | Tu 
nl or } | ky) Hin — 932 One ———______) — In) al 
Za iat Ci HT #6] TITTI wx! Ale 
$i k 
(co) n 
“Lx k) Uj — 95x OK Pr = ; s k ve Mr» F 
+ Bley Pee Detlev et o 0 


Clup= n 0 COMI LX Upg = n 
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It is evident that for large n and sufficiently large yw, the last term may be 
neglected in comparison to the first. Consequently equation (4) is approxi- 
mately fulfilled for x = x, too. 

If we substitute in (22b) for J~! the expression (7) and than perform the 
integration over 2 we get an expansion of @ in powers of 4; the n-th power 
coefficient reads 


Laat LT a 
SATO Tat 


j=l % 


(26) (n) 


n = Il 
where summation is restricted by (7a). Since 


d0(n)- 
P(e, Lo 5 n) = O2(x, (a |2) 
relation (21) follows immediately from (26) by comparing it with (15e). We 
see also from (12) that our @(n) is the first term in the alternating series of 
(n) divided by the number of terms in this series. 


4. — Transition to continuous t-variable. 


41. — We have shown that the exact solution of an equation of diffusion 
type is given by (12) and an approximate solution, proper for large n, has 
the form 


is 1 9 IO Li 
(27) (a, Lo; N) cena TT ’ 


where summation is restricted by (7a). This all was done for the case when 
t-variable takes discrete values. Let us investigate now the case when At 
tends to zero. 

We define the function p(s), s= 0, 1; p(0)= 6, p(1)= a, where a and D 
are arbitrary complex numbers. Further we introduce by linear transformation 
a new variable 
s(a+ b) — 

V ab ab 


E =- 


We write p(s) =q(&). It may be shown, as usual, with the help of characte- 
ristic functions that for 


ste 2 aE) DIG) 
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where the summation is taken over &, satisfying 


ee 


aye 


a central limit theorem is true; indeed, for large n, 
(a + by (a +b) (2 a_\\° 

(28 P(é) = exp | — ———— .— ; 

Cr v ‘Qanab i 2 a+b 


in other words, 


n! (a5) exp|— (a + 6)? (o na J] 


29 ~ ine E SSE, arbr-*~ 
(29) al(n_ a)! x. V27nab 2nab 


We apply this result to perform the computations in (27) for At + 0. We 
confine ourselves to the case when the number of levels is bounded or at least 
the concerned spectrum is discrete, the investigations of the case of continuous 
spectrum meet some obstacles. 


The most suitable case for computations is that of stochastic processes, 
I 


because all 2(j|k) >0 and (a+b) = Y 2(j|k) =1. If we consider the sum- 
dei 
over-states some complications already arise, although the condition #(j|k) > 0 
holds too. The most troublesome case is, however, that of Quantum Mecha- 
nics, because 2(j7|k) are complex numbers (*). Let us pay some attention to 
this case. 
On account of (29) we are able to replace for large uy, 


Un! D È 
(30) IT an II sie, 


by 
4 
(31a) (2 2 (j|k)) ) He ( (Qrru,) =D! 2 o-1- 


ak. pe 2(1]k) Hee | ‘| 
ang 1 Vail: Ae Selo |? 


‘exp|— 


(*) In virtue of (A) and (B) as well as of the continuity of the solution in respect 


U 

to the t-variable we have for each representation >¥ 2(j |b) =1+0(Ad); hence normal- 
j= 
ization as in case of probabilities holds for At +0. 
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where Q7 is the reciprocal quadratic form to 


Q(1,, «5 Ci.) 5 SO Sa, dip ow Ris Tye on YY 


=07Tg=0 1T,-,=0 


and o? is the determinant of the matrix of «dispersion » o, .,, belonging 
Ps IO Uh 


to the «distribution function » M(S,, 823, 8i-1)5, S19 Say -+-5 Sia = 9, 1, 


p(0, 0, ...,0) =e (116) (¥ 207 |) 


p(0, 0, .... 1) = 2(—1|k) (de (9 \\k))- 


j=1 


otherwise =O. 


However, (29) enables us to put (30) into another asymptotic form; indeed, 
we have 


DI (Nix — 5-1 Di)” 


LA 1-1 
)|h)) & dan; +-Aan0 se —t ex 
ae) (> #01) I CCE ENTI fj nied sx 


j=l j 
where 

Ny, = Ln — SI [ork j=0,4,..: (= 1; k= 1, 25 
and i. 

1 
3 a(r|k 
ae ci ik) ia Si 1 i 
ZA, 4 Pegli 


the additional condition (7a) reads now 
L L 
DI 7 T > Nin = Noi - 
k=1 k=1 


The essential feature of both forms (31) is that they have the form of a mani- 
fold Gaussian distribution and that the dispersion is proportional to yy, or 
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n;+, resp. in the exponent. If we go over to continuous # variable we must 
carry out the transformation 


| Mix IN ESF a 
(32) = Ie ba 11526 


| (22 At = ti, 3 
or 
Nix Nt tee 


Nok Ate tie 


resp. Thus we get a new Gaussian distribution in both cases, but with dis- 
persion equal to t,o? At or t;_1x%;,x%;x At resp. proportional to At. It is obvious 
that, if At tends to zero, these distributions will tend to Dirac delta. We shall 
not, however, replace these expressions straightforward by Dirac delta’s; we 
see that the main contribution in (27) comes from the region 


t.2(7 | k ; 

(33a) in & dEi in (31a), 
> e(r|k) 
r= 

(33) Podi aio) 


Since Yu, = > py, it follows immediately from (33a) and also, after short 
c 


o) 


algebraical calculation, from (330) that these t, must fulfil the equations 


(34) Skt x ESE 


from these equations we can find t;. Thus, of interest are the contributions 
due to groups of paths, which show a deviation in regard to the stochastically 
typical group of paths, leading to stationary state. Computations of these 
deviations will be the subject of a separate paper; the starting point for these 
computations will be the equations (12) and (27) as well as (31). 


ON THE MODIFICATION OF FEYNMAN’S « INTEGRAL-OVER-ALL-PATHS ) METHOD 75 


42. — We give now some crude estimations for the case of transition pro- 
bability amplitudes of usual form in x representation, i.e. given by (16) and (155). 
What belongs to the factor 


(Det) 


in (31) one gets 


x(a! | ee att oy nll pee mi DA I |" 
Oe | a) #2 (i) | dw’ exp oF Al (e' — 0)? == V(x yar oe 
= exp |— LV (0) My Rx FT wrt 


Here we have employed the method of stationary phase as well as that 


(A) | exp bi 


What belongs the solution of (34) for t,, ie. the stationary state, so 


e 


(B) 2(0' |a) = dla — w')[1 + O(At)] + (OAT) , 


consequently all ¢, are approximately equal, hence e.g. 


(35) Il exp 


av’ =1 


i ; dI noe 
AAC Ju. |= exp ia fre Ja" 3 


These results change from one representation to another; e.g. for k repre- 
sentation we get 


1 


(Da |) ( | dk! C 


k'= 


exp | Uh iH i exp | 5 t, (È i ro) : 


i 2m 


4°3. — In connection with representations we mention that we are able 
to construct 2(j|k) in any quantum representation, if only the transition pro- 
bability amplitudes are known in solely one representation 


m 3 Telly It I mi ‘pal m\2 
(carnal) [fe da <a' |a') exp OF AI (x’ — 2) 


44. — Further we mention that the idea, presented in this note, may be 
immediately generalized from one to several dimensions. To this aim the 


a 


-V (a) At| <a |a> = 2z(a'|a). 


> 
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indices j and k in ¢(j|k) as well as in m, must be replaced by vectors, e.g. 


J Si (Î1, oor 44) 5 k= (ky, sony ky) . 


5. — The usual perturbation calculus. 


If the transition probability amplitudes have the shape of (16) and (15d) 
and if VAth-!<1 or tends to infinity for considered set of states, as e.g. in 
the case of helium, it may be shown that the application of the method of 
slow varying factor in the u,,-space leads to results which are obtained by 
means of the usual perturbation method; we must, of course, work in the 
Schrédinger picture, where interaction energy is independent of time. This 
new approach to the known solution of perturbation method may call some 
attention; perhaps other mathematical conditions as usual are required for 
the approximation to be valid. 

Finally we mention that from the crude result (35) one may conclude that 
the first approximation step of the perturbation method yields the solution 
in form of a product of the undisturbed solution and the unitary factor given 
by the right hand side of (35), which is of course physically meaningless. We 
get rid of this factor by normalizing the potential V(#) in such a way that 


fre da = 0% 
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1. — Introduction. 


When an elementary particle, such as a proton, meson, etc., strikes another 
particle or group of particles, one or more new particles may be created in 
the process. For a given incident particle and a given target particle, particles 
chosen at random may or may not be produced. For one thing, very small 
cross sections are connected with the direct production of neutrinos, electrons, 
and muons as well as with the production of photons in many instances, so we will 
be solely concerned here with the production of particles more massive than 
muons. In addition, if the total rest mass of the product particles is greater 
than the total rest mass of the initiating particles, conservation of energy 
requires that the incident particle have a certain minimum energy, the thre- 
ghold energy, for the reaction to proceed. It is further presumed that electric 
charge must be conserved in the reaction. Furthermore, two other quantum 
numbers (1), called here the « attribute » A and the «heaviness» B, which 
refer to intrinsic properties of the particles and which add algebraically for 
a state of more than one particle apparently have to be conserved (2) if the 
production reaction is to occur with any appreciable cross section. Table I 
gives the A and B numbers for the known and suspected particles to be consi- 
dered here. 

B, the heaviness quantum number, seems to be rigorously conserved in 


(1) M. GELL-MANN: Phys. Rev., 92, 833 (1953); T. NAKANO and K. NISHIJIMA: 
Progr. Theor. Phys., 10, 581 (1953); R. G. SACHS: Phys. Rev., 99, 1573 (1955). 

(2) See Proceedings of the Sixth Annual Rochester Conference on High-Energy Physics 
(New York, 1956), for a discussion of the experiments relating to conservation 


of these numbers. 


78 C. O. BEASLEY jr. and w. G. HOLLADAY 


TABLE I. — The «attribute» A and the « heaviness » B for the class of strongly interacting 
particles. The bar over a particle indicates an anti-particle. The numbers in the squares 
are the masses in GeV of the particles used in the threshold calculations. 


A 
1 0 == 1 
A° 
1 PN 1.115 == 
939 Sexo 1.321 
1.189 
Bie so K*K° mtn KK? 
493 .139 493 
AO 
=] E = 1.115. NP 
1.321 Dei .939 
1.189 


all reactions. A, which appears in the expression (*) q= 4/2 + 73, where ¢ 
is the charge of the particle and 7; its third component of isotopic spin, seems 
to be conserved in all « fast » reactions such as the collision processes contem- 
plated here, but is not conserved in the «slow » processes among the part- 
icles in Table I. 

For a given particle being bombarded by a given incident particle, many 
possible reactions (called here valid reactions) exist consistent with the conser- 
vation of the above-mentioned quantum numbers, and, of course, there are 
many others (called here invalid reactions) which do not conserve these numbers. 

From what has been said, it is easy to answer the question as to whether 
a given set of final particles will be produced from a given set of initial part- 
icles. If the total charge, attribute, and heaviness are not conserved, the cross. 
section for the reaction will presumably be so small that the reaction will not 
be observed, but on the other hand if these quantities are conserved, it is 
expected that the reaction will occur, provided, of course, that energy-mo- 
momentum conservation holds (*). 


(3) A number of people have emphasized the efficacy of this relation. R. G. Sacus: 
Phys. Rev., 99, 1573 (1955); B. D’EspAaGNAT and I. PRENTKI: Nuclear Phys., 1, 33 (1956); 
J. Tiomno: Phys. Rev., 103, 1589 (1956). 

(4) Since our initial and final states are in the continuum, no account need be 
taken of parity and angular momentum conservation. Furthermore, for the particles 
listed in Table I, conservation of B automatically takes care of the selection rules based 
on conservation of angular momentum which would forbid, say, an integral and a 
half-integral spin particle from making two half-integral spin particles. 
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Not as easy to answer is this question: for a given pair of initial particles, 
what are all the final sets of particles which may be produced? Even if atten- 
tion is restricted to those final states involving no more than three or four 
particles, the possibilities are so numerous that a haphazard attempt to write 
down all the final state products could very well leave some unsuspected but 
possible reaction unlisted. In the next section of this paper we provide a 
systematic procedure for writing down all the possible valid reactions for a given 
initial pair of particles. The procedure automatically eliminates the invalid ones. 

This procedure has been applied to determine the possible reaction pro- 
ducts of a nucleon, a pion, photon, electron, kaon, anti-kaon and anti-nucleon 
incident upon a nucleon. These reaction products are listed in Table III for 
reactions including as many as four product particles In addition, the threshold 
energy has been calculated for each of these processes under three different condi- 
tions: 1) when the target is a free nucleon; 2) when the target nucleon is 
in the nucleus, in which the nucleon is considered to have kinetic energy of 
25 MeV with momentum toward the incident particle (9); 3) by means of 
a double stage reaction first considered dy FELDMAN (6), in the first stage of 
which a nucleon strikes another nucleon (at rest or in motion in a nucleus) 
and creates a single high energy pion, which then becomes the incident par- 
ticle for the second stage of the process. Since in process 3) there are two 
target nucleons involved, the threshold energy for a certain reaction by way 
of the double-stage process is lower than that of the single stage process in 
which an incident nucleon strikes a target nucleon and creates the products 
directly. This lowering of the threshold is somewhat analogous to that which 
occurs when the mass of the target is increased; for such an increase there 
is a decrease in energy lost in center of mass motion. 


2. Algebra of elementary particle reactions. 


In this section we wish to provide a general procedure for writing down 
all valid reactions involving known or suspected particles heavier than muons. 
A valid reaction is defined to be one which conserves the numbers A and Bb. 
We will not distinguish the charge states of a given particle, but if that is 
done, the procedure can still be carried out. 

Each particle is given its designation by its A and B numbers; this design- 
ation for a particle is called a term (7). For example, the term denoting 


(5) W. McMirran and E. TELLER: Phys. Rev., 72, 1 (1947). 

(6) G. FELDMAN: Phys. Rev., 95, 1697 (1953). 

(7) Each term represents a particle which may have different charge states. It 
should also be pointed out that the (01) and (0—1) terms have further non unique 
properties since they represent both A®, & and A°, X respectively. This fact should 
be kept explicitly in mind when the techniques presented here are applied. 
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the nucleon is «11», that of the pion «00», or that of the K particle « 10 ». 
The A quantum number is the first number of the term, and the B quantum 
number is the second. There are nine possible terms. We define a set to be 
any group of these terms. Now each of the nine terms may be expanded into 
a set of two terms, the sum of whose A and B quantum numbers gives that 
term which was expanded. This expansion just gives pairs of particles into 
which a given particle may transform and still conserve the A and B quantum 
numbers. In Table II we see a listing of all possible expansions of the nine 
terms into sets of two terms. 


TABLE II. — All possible expansions into sets of two terms. 
11 +11 + 00 10> 114+ 0—1 1-1- 10+0-1 
10 + O1 10 + 00 00 + 1—1 
1—1+401 
01 +00 + 01 00 > 00 + 00 0—1+0—1+4 00 
11 + — 10 O—1+ 01 —10+1—1 
10/211 — 10210 | CIT 2410 
—1-1+11 | 
—ll+1—1 
11-01 10 —10--11+4 0-1 —l-1- 10+0-1 
00 + — 11 — 10 + 00 00 + —1—1 
—1—1+401 


Now suppose we start with a given initial pair of particles. First, all valid 
sets of two product particles are easily written down by methods of exhaustion. 
Now we wish to find all valid reactions in which the final state consists of 
three particles. To do this, we take a term of one of our sets of two product 
terms and expand it by using Table II. This expansion of a single term, plus 
our unexpanded term, gives us a set of three terms. To assure that all valid 
reactions are obtained and that duplications will not arise, the following rules 
have been formed for the expansion of a group of n-term sets into a gronp 
of n-+1-term sets. These rules are stated here without proof, but a proof of 
them may be found in the Appendix. These are the rules: 


1) All valid reactions in which n-+1 particles are produced from a given 
pair of initial particles can be formed by expanding all the different terms 
of the valid sets of n terms. 


2) Im the expansion of the terms of a set, the expansion of the first 
term being expanded in that set may contain a 00 term in its expansion; there- 
after, to avoid duplication, no expansion of any of the remaining terms of 
that set may contain a 00 term. 
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3) A set of n+1 terms from the expansion of a set of n terms has al- 
ready appeared if n — 1 terms of the set of n-+1 terms are present in one of 
the sets of n terms which has been previously expanded. 

We shall use as an illustration of the procedure a nucleon-nucleon (N-N) 
reaction. In this case, there is no valid reaction which leads to a final state 
of two particles different from the initial state. To obtain the three-term sets, 
we note that there is only one term to expand, since both of our only set of 
two terms are 11°s. Therefore, all valid reactions with three product particles 
are 


(1) 11 + 11 + 00 11 +10 + 01 


These represent the product particles of the reactions 


(2) N+N-+2N+4x and NE Ne NK Se 
SN eC Ln 


Now we have all valid reactions with three product particles. From the sets 
listed in (1), we may obtain all valid reactions for four product particles by 
two of the three terms fixed and expanding the third according to 


holding 
In our example, we see from our two sets of three 


the same expansion rules. 
terms that there are five different terms to be expanded in order to get all 


possible sets of four terms. In the expansion of the 11 term, in the first set 


of (1), we get 
(©) 11 +11 + 00 + 00 11 = 00. 10 =F OL, 


and in the expansion of the 00 term, we get 


reer SONE eee LL Tese 
(4) i 
Ben se ae 0) 4 10 11+11 11 +1—1. 


Notice that one term of the expansion of 00 was left out because of the appli- 
cation of rule number two. Now in the expansion of the 11 term of the second 


set in (1), we get only one term 


(5) ee Ole 001 


the other having dropped out under rule number three. Likewise, the expan- 
sions of the 10 and 01 terms have only one term each; all other terms are locked 
out by rule number three. The two unique terms are 


ME 1-01, (LL FAO 11+ 10, 


(6) 
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Thus the nine sets in (3), (4), (5) and (6) are all of the valid reactions with four 
product particles. The procedure may be extended to include the 5-particle 
product reactions and so on. 


3. — Threshold calculations. 


In this section we record the formulas used in calculating the threshold 
kinetic energy 7, in the laboratory system for the various processes. If E, 
is the total laboratory energy of the target nucleon with mass M, which is 
assumed to be moving directly toward the projectile which has mass M,, and 
if M, is the mass of the i-th one of n product particles, then straightforward 
relativistic mechanics with conservation of energy and momentum yields 


_ AE, —[(4? — MiMs)(Ei — Mi)} 


(7) T, mi 


SUI ; 


where A = 3[(> M,)?— Mi — Mj]. Here, c=1. If the target nucleon is 
t=1 


at rest, then #,=—WM,, and 


(8) Li 


1 T \2 
2M, (2 MI) cò (Mi = M,) | . 


From Eq. (7) and (8) all threshold energies for a single stage process can 
be calculeted. However, as mentioned in the Introduction, there is a double 
stage process which yields somewhat lower thresholds than those for the single 
stage nucleon-nucleon process. In this reaction, a nucleon strikes another 
nucleon creating a pion which in turn creates the new particles. The threshold 
energy of the pion-induced reaction can be calculeted by Eq. (7) and (8). Then, 
the problem remains to find what minimum energy a nucleon must have to 
give a pion of this energy. 

If (E,, P.), (E., P.) and (E_, P_) are the total energy and momentum for 
the target nucleon, the projectile nucleon, and the produced pion, this minimum 
energy is given by 


(9) BH, =~ (1+ 38) +3[1+991+9}F, 


where 9 = E_— E, R=P_+ P,, S = M?/R? —Q?, and M = nucleon mass. 

From Eqs. (7), (8) and (9) thresholds for all sets of two, three, and four 
product particles have been calculated for the types of projectiles and targets 
‘mentioned in the Introduction. These thresholds are listed in Table ITI in 
order of increasing energy within each set. 
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TABLE III. — Reaction products and threshold kinetic energies (GeV) for various initial 
pairs of particles. The column labeled N refers to a target nucleon at rest. The column 
labeled N refers to target nucleon with 25 MeV of energy moving toward projectile. 


v+N 


nen NN =2N+r]iN4N>2N4 7 
Reaction e-EN (*) | | ct N a+ N 
products Ee ware a fo 
N N N N N N NN 
N+r .15 SEO 29 DOME. 14 
K+A° 91 To .76 Se ell Sal 0.04. 0,967 1 0276 
K+ 1.04 JD .89 GS | 1% 1.02 1.08 | 0.87 
N+2n .32 .26 Ml; SR 160 39 032 025 
K+A°+r 1.15 .91 1.01 ‘7700 1:39 06:13 81:20 0.96 a) 
K+Z-+r 1.29 INOS outa. Ls | ey | ey A | 
N+K+K GOs 19 oo ele aa CRA deg DILSG 
2K+E 2.36 eee | ETO 60 ORI yee aie, 
2N+N 3.75 | 2.98 | 3.61 | 2.88 | 4.10 |-3.30-] 3.9 | 3.0 
N+A°+ A° ASS 0 3 ATL 05-3 land s3 5.1 4.0 
N+A°+% 5.13 4.07 | 4.98 3.92 5.5 4.5 5.3 AS 
NEX+A° RO re RO ALOT MO 4.5 5.3 4,2 
NELZA-X RR O TORTE e oe) 5.6 Ado 
2A°+ 6.24 | 4.95 | 6:10 | 4.81 6.6 | 5.3 6.4 5.1 
N+E+E 6.35 «| 5.04 | 6.21 -| 4.89 | 6.7 5.4 | 6.5 5.2 
MA+Z+È 6.53 | 5.17 6.38 | 5.02 | 6.9 5.5 6.7 5.3 
22+E 6.82 5.40 G:C Tun 22 5.7 7.0 5.6 
N+3n 51 .40 36 .26 | 0.68 | 0.57 .52 Al 
K+ A°+ 27 1.41 Lal Oa 27 98 | 1.66 1.36 1:48 WS 
Dp eae ee are ey i A od I lB E VASTI? 1.30 | 
MEN re Ko 1:79) 1:42 7) 61.65.) 1:28." 2.07 1.68 1.86 | 1.49 | 
2K4+r+E Dre | etto 2:56 2:01 Tent 24 23 2.200 
2K+A°+K I | 2.47 192.97) 2.880) 3.4 || 2.8) 08.2 2.6 | 
2k eK = agli 263°) 02:17 12.49 0003.0000 72.0) ast: ZI 
2N+r+N 4:18.) 3.31 | 4.03 | 3.16 | 45 3:64: 3A 
N+r+A°+A° | 5.36 | 4.25 (eS Ome 090s. Ta O o Nik Saal 
Neer Sa A | 5.61 4.45 | 5.46 | 4.31 6.0 4.8 | 5.8 | 46 | 
N+r+A°+% 5.61 | 4.45 | 5.46 | 431 6.0 AS SO 
Nee dae Sea | 5.88 4.66 5.73 | 4.53 6.2 | 5.0 | 6.0 | 48 
pal ICG 5.09 E 4 76 |. 5.85.) 4:62 64° 5.1 oe | O 
2N+A°+K | 5.99 4.76 5.85 4.62 6A. eb. 1 6.2 | 4.9 | 
NIPZ-ENK 627 | 498 | 6.13 | 483%) 6.6 | 5.3 6.4 BL | 
2N+E+K | 6.27 | 4.98 | 6.13 | 483.) 66 | 5.3 | 64 Lat | 
2A9+E+r | 6.77 | 537 | G:62N O oe) I al 5.7 6.9 | 5.5 
N+E+E+r Giggoi ento does 1a | 5.8) | 71 5.6 
NOES eee | 7.07 5.60 6.92 | 5.46 7.4 COM arias 5.8 
a+ 7.36 nee i oto 6.69 Te wl Lie | 7.5 6.0 
|} 2A°+A°+K | 7.36 ee) a Se OTT TI 6.0 
N--E-4+A°+K | 7.49 | 5.94 | 7.35 | 5.80 | 7.9 6.3 TI 6.1 


| (*) An, electron should be included as a-reaction product for this reaction. 
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TapLe III (continued). 


Reaction 
products 


N+A°+E+K 
NESSO 
NOSED TEI 

NES ESE 
N+Z+E+K 
DSL ACL K 

TNS Ep ya K 
ps SNK 

AQ BK 
Sn BALK 


Reaction 
products 


NE ET 
N+A°+K 
N+X+K 

2N +27 
N+A°+K+7x7 
N+24+K+472 
2N+K+K 
2A04+2K 


N+2K+E 
2K+A0+Y 
2K +25 


3N-+N 

INL AP AO 
IN + APS 
JIN DHE IN 
INS SESS 
N+2A°+38 
N+E+E 


Ci O. 


+N ee IN+N+2N+ xiN+N-2N+r 
e+N (*) | ) t+N =+N 
N N N N N N N N 
oa Pe a 
740" 6.94. 0.55 We 80 a7. 6.3 rove 0 
7.68 | 6.09 | 7.53 | 5.93 | 8.0 6.4 7.95 az 
7.680 \0 6:00 O 6.4 79 | 62 
es 6.2 7.66 | 6.05 | 8.1 6.5 he 8.0! a 64 
7.8 6.2 7.66.5 6. 0G mes. 0 6.5 8.0 ts. 6.4 
8.0 6.3 78 6.2 8.3 6.7 8.1 6.5 
8.0 6.3 | 7.8 6.2 8.3 6.7 8.1 6.5 
8.3 Mie) SRE 6.4 8.7 6.9 8.5 6.7 
9.1 Fes eee 0.) Te emilee O56 7.6 9.3 7.4 
9.5 hess 9.3 TAY ie SOS -Ealey.0 9.6 TAI 
| | 
NE Reaction = Ke 3 
x | N products N da n 
— ES cel 
0.29 0.14 N+K+r 0.22 + 0.12 
1.58 1.10 A°+2K 1.26 0.91 
175 1.26 x2+2K 1.43 1.05 
0.59 0.35 2N+A° 3.68 | 2.82 
1.96 1.40 IN +S 2.01 as Ok 
2.17 1.57 N+A°+8 4.97 3.85 
2.49 1.81 N-EZ$E 5.24 4.06 
3.63 2.70 N+K+2n 0.46 0.29 
A°9+2K+xr 1.58 nen 
Z+2K+r 1.76 Jesi 
a: 278 N 22K kK 2.03 1.52 
3.88 2.91 E+3K 3.09 2.36 
4.15 3.11 PANE Ar 4.12 3.18 
5.63 4.29 2N+Z+xr 4.37 3.37 
7.10 5.45 2N+N+K 4.950.) 3.67 
7.43 5.72 N+A°+E+xr CATO NO, 
7.43 5.72 N+E+Z+xr 5-16 0447 
7.76 5.98 Nap Ata AC Kee |, (6.050, 1 47d 
3.97 A 26.8 2N+E4K MER VA MOZZO 
9.0 7.0 NASSP K le 6735) se 
9.2 TOA N+4Z+A°+K | 6.35 4.94 
9.6 7.4 N4Z4 XK | 6.65 5.17 
2A/9+E+K MRE 5.95 
N+E+E+K 105 6.05 
AD SEEK 7.94 6.21 
22+E+K 8.27 6.46 
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(*) An electron should be included as a reaction product for this reaction. 


TABLE III (continwed). 


Reaction 
products 


N+A°-+N 
S+N+N 
2A°-+ A°0 
N+E+A° 
2A°9+% 
AOL AOL 
NESS 
TWEEVDS, 
A°+X+Y 


254% 


KLK4A°+r 
KERESLa 

2K+K+N 
9kK+K+e 


Reaction 
products 


A°+ AP 
D+A° 
A9+% 
uh 
Bs 


N+N+7 


ELEMENTARY 


K-+N 

N N 

Exothermic 

Exothermic 
66 | 45 

Exothermic 

Exotbermic 
DD: 12 
.93 .66 
1.26 .92 
1.43 1.05 
3.68 2.83 
3.92 3.02 
487 SOT 
4.97 3.86 
Hale 3.98 
5.13 3.98 
5.24 | 4.07 
5.40 4.21 
5.40 4.21 
5.68 4.42 
6.42 .| 5.01 
6.72 5.25 
sis .07 
Za 165 
46 .29 
1723 .80 
IST 1.16 
Teo: Rod 
2.02 162 
3.08 2.36 

N+N 
N N 
| 

ET | 0.48 
0.95 0.62 
0.95 0.62 
Telo 0.76 
1.84 1.30 
0.29 0.14 


PARTICLE REACTIONS 


Reaction 
products 


N+A°4+N+7 
N+Z+N4+x 
2N+N+K 
294 AP Lr 
NEPI 
DINO AL DY ET 
AX+ APD +4 
NELLO La 
py EIN Ee 
NOMS EAL as 
N+ A°+ A°+K 
2A°4N+K 
N4E+N+K 
DEE 
Ke Nee Ace 
Ke NEA 
KEN aoe 
K+N+3423 
KEN 22s 
AE+E4+r 
S+8+847 
K+2A°+8 
K AO N+E 


Tec NOS Ss Se 
K+A°+2+8 
K+A/°+3+E 


Reaction 
products 


A+A°+K+K 
N+E+2K 
E+N+2K 
Z+A°+K+K 
A+Z+K+K 
DICK 


29 
.61 
.61 


J +] -J +I 


a 
CS 
re 


8.26 
8.26 
9.48 


K-+N 


N+N 


iN 
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TABLE III (continued). 


Reaction N+N 
products 
N N 
A9+ M47 1.11 0.74 
M+ A°+z 1.30 0.88 
Aad +7 1.30 0.88 
pelea, 1.49 1.03 
ASL N-+K | 1.58 1.10 
N+4A°+K | 1.58 1.10 
X4+N+K 1.78 1.24 
NESSO 78 1.24 
B+84+7 2.23 1.61 
BAe Kk 2.69 | 1.97 
AO+E+K 2.69 1.97 
Bee ee Kk 2.93 2.16 
SL EIR 2.93 2.16 
2K 9 K 0.19 O77 
N4N+2r 0.59 . | 0.35 
A°+A°%+2r 1.46 1.01 
A°+Z+2r 1.66 1.17 
MH4A°+2r 1.66 (oly! 
pe Ea eee. 1.33 
A°+N+K+x | 1.96 1.40 
N+A°+K4+r 196 | 1.40 
Z+N+K4+r Dr 1 yi 
NAS Keer 27 = a melas? 
N-N+K+K B® |) Th AT 
E+E+2r 2.65 1.94 
Hel AOE Kx Sls. || BERD 
MO+8B4K4+7 3.13 2.32 
E+Z+K+xr 3:37 2:50 
ZS+E+K4+xr 3.37 2.50 


Elaol ano Tec 


HOLLADAY 


Reaction 
products 


E+E+K+K 
2N +2N 
N+A°+N-+ A? 
N+A°4+N4% 
INCESE SNC, 
NE SPO 
20420 
2A°4N+E 
N+E+2A0 

| N+E+N+E 
2A9+A°+% 
ALS +20 
N+E+4A49+4% 
N° Z+N+E 
DAP Os; 
A®°+ 30-4 A945 
22+2A° 
N+E42Y 
22+4+N+È 
AGED oy 
20-+A°+ > 
22) 
A°+ 5+ A+ 5 
AETISTTE 
T+E+A°+E 
DLEpSEE 
28+42E 


| 


| 


Or 
i 
ow 


I 
i Db 
= 
OQ 


FS 
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I 


lor) 


Sey Sy cape 
i 
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9.8 


10.2 
10.8 
TELI 
i 
11.5 
13.0 


dI 


SS 
(sy fo ice) 


SI 
SD 


Ses I 


Cap ea 
WO O & BR w dd wu 


UU 


cS 00 00 00 
Sa 


10.1 


(*) The exothermic reactions N + N + up to 13 pions (inclusive) and IN SL INI nike IR up 


t06 pions (inclusive) are not listed. 


Fig. 1, 2 and 3 give an easy means of finding approximate thresholds if 
the total mass of the product particles is known. Each curve is labeled by 
the type of initiating process. Representative reaction products are given 
on each curve. These graphs would be useful in determining thresholds for 
reactions involving more than four product particles which are not listed in 


Table IIT. 
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Figg. 1, 2, 3. - Threshold kinetic energy as 
> a function of the total mass of the reaction 
* 

tar Sa Ao products. Each curve is labeled by the 
È := | nrojectile and the target. The symbol N 

na Ss = pro) 5 Wi 
Lo | | SS represents a target nucleon with 25 MeV 

= . 

È 5 | == Ch of energy with momentum toward the 
= | = ee a projectile. Representative reaction pro- 
= | Sa ducts are given on each curve. A. bar 
= = over a particle represents an anti-particle. 
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Table IV gives a brief listing of the lowest threshold for the production 
of any elementary particle for a given production process. For each particle, 
other product particles (with the exception of an arbitrary number of nucleons, 
a number which depends on the number of nucleons in the initial state of the 
process) produced with that particle are listed. 


TaBLE IV. — Lowest thresholds for a given elementary particle for a given process. Bars 
over letters indicate anti-particles. An arrow over an N means that the target nucleon is 
in a nucleus and has KE of 25 MeV with momentum toward projectile. 


ET Reaction | Di y+N | y+N | N+N | N+N n N+N + N+N = 
(°° | products [N+N|N+N|e+N|e+N|->2N+x > 2N+n| >2N+n| >2N+7 
ico NIE | CAO ea N Ne 

MR 158) 10 091 07 I 00 0.96 0.76 
K | KK | 24900181 | 1.90 171.199) ares Wie 1-46 1.56 1.26 | 
N | N+N | 5.65 | 4.29 | 3.75 | 2.98 is! 3.3 3.9 3a 

MEO SAS pd 539 1080-91 | 0.72 1.13) |“ 0.94 0.96 ON ON 

TIES A SASSI 5.45 4.88 3.87) 5.3 4.3 Bu 4.0) 

| = | K+% 1.78 | 1.26.| 1.04 |:0.82| 1.27 £02" )y 31.08 0.87 

| x | AS-D | 7.4 DT | 5.13 | 4.07 | 5.0 | 4.5 Dao 4.2 
© | Wise) | Gass |) Maric | este TS eae} 2.15 2.43 1.97 
He | o2A0-S 6.0 16.8 | 6.2 14.05 e 6.6255) 15.32 es6.4 eet 

(*) Nucleons as a reaction product are not listed 
(*) An electron is considered to be a reaction product of this reaction. | 


4. — Conclusion. 


On the assumption that the quantum numbers A and B are conserved, 
a procedure has been devised for writing down all the reaction products com- 
posed of particles in Table I. Table III lists these reaction products up to 
four particles together with the production threshold when the projectile is 
a photon, electron, pion, kaon, anti-kaon, nucleon or anti-nucleon and the 
target is a stationary nucleon or a 25 MeV nucleon with momentum toward 
the projectile. Threshold energies have also been calculated for the Feldman 
two stage process. 

As expected, lowest thresholds for the production of any given elementary 
particles are associated with the e+N or y+N process, although the difference 
of thresholds of that process and the two stage process becomes less marked 
at high energies. The two stage process lowers the threshold over the analogous 
one-stage by roughly 20% to 25%. Furthermore, in the single stage process, 
the threshold for a 25 MeV target nucleon is some 25% lower than that for a 
target nucleon at rest. 
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It seems worthwhile to point out that the most massive anti-hyperon, 
the =, which can be made by the process, energy +N > 2A°+ &, has a thre- 
shold by two the stage process of about 5.3 GeV, and therefore presumably 
can be made at Berkeley. 


APPENDIX 


The rules stated in Sect. 2 for the expansion of a set of n terms into a set 
of n-+ 1 terms will now be proved. Rule 1 essentially states that any valid 
set of n +1 terms (which represent, of course, n + 1 particles) can be obtained 
by the expansion of some one term of one of the sets of n terms. 

Analytically, what is meant by a valid set of n-+1 terms is that 


nHt 


sb 2 
xii = DI oe, Bs ; 
d é 


i i=1 


Il 


where «,f; is a pair of the integers (1, 0, — 1), which, of course, represents 
one of the particles in Table I, and «A, and a,b, are the AB designation of 
the two initial particles. The sum indicates that the as are to add algebraically 
as are the f’s. Now rule 1 will be true, if some two terms, say cf + xntiPat1; 
of the first n + 1 set, combine to give one of the nine possible terms, say ab, 
representing a particle in Table I, for then 


n+1 n—1 

> afi = Daf + ab, 

AIN col 
but the expression on the right is a valid set of n terms, such that one of the 
expansion of the ab term yields our set of n +1 terms. Thus, our problem 
is reduced to one of proving that in a valid set of n-+1 terms, there are 
always some two terms which, when added, will give one of the nine terms 
representing a particle in Table I. 

The proof will be subdivided into two cases. In the first case, we suppose 
that one of the numbers of some term of the n+1 set is a zero, say OP. 
Now if the second number of that term is also zero, obviously the term will 
add to any other of the n terms to give a valid particle, since 00,0 aca: 
If the second number happens to be a 1, (the proof proceeds in the same way 
if this number is a — 1), then all the remaining n terms cannot have PSs Toe 
their second numbers, as the sum of the n terms would give a term with an 
n in it (8). Of the initial-state particle set, the one with the greatest sum of 
A and B quantum numbers, the N-N reaction can only be 22. Hence, one 
of the remaining n terms must have either a 0 or a—1 for its second number; 
in either case, that term added to our initial term (01) will result in a valid term. 


(8) Throughout we are assuming that all two term sets can be written down by 
the method of exhaustion and that n > 2. 
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For the second case, we assume that none of the n + 1 terms has a 0 in it. 
Now let us refer to Table V below, in which are listed all possible terms 
without 0’s. 


TABLE V. 


11 =) 
la = eh 


Note that the sum of any n + 1 terms in the same column or same row gives 
a non valid term with + (n+ 1) in it. Therefore, either the diagonal or anti- 
diagonal sets of terms, or both, must appear. But it is obvious that the sum 
of either of the diagonal or anti-diagonal sets of terms gives a valid term, 00. 
Hence, two of the terms must add in this case to give a valid term. 


n 


With regard to the second rule, consider the expansion of the set SI asp; - 


=a 
One expansion of «6, is. x}, + 00; likewise, an expansion of «6, is x}, + 00, 
and similarly for the rest of the terms. But any one of these expansions added 
n 
to the other terms of the set gives > «8; + 00, so that only one such expansion 
t=1 
need be made. 
The third rule may be proved by the following line of reasoning. For any 
given pair of initial particles the sum of terms representing possible reaction 
products must remain constant. Therefore for a set of » terms 


(1) (af, sony a8.) ’ 


which has been completely expanded into n + 1 terms, and for another n-term 
I rf L; I . . 
set (81, ..-, &Pn) which has an expansion 


. Li 
(2) (2181, O05 5) CRE ae a Andy, 5 On+10n41) ’ 


we have 


n n-1 


(3) > ai; = Di cibi + a db, id +10 nti o 


t=1 (eal 
Now if «6; = «fi, when i goes from 1 to n— 1, then from (3) 


KIDS = 6.6, + @ 1b 4, 
Ted b, + 4:40 1 +, is one of the possible expansions of «,8,, and thus the 
expansion (2) has already appeared ‘as one of the possible expansions of (1). 
This completes the proof of rule 3. 
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1. — Introduction. 


Even before the experimental evidence, the multiple meson production 
had been discussed in connection with the criticism for applying the limits of 
the field theoretical description to a high energy region (1), hoping to find there 
some new conceptions or frames of a future theory (2). For the moment how- 
ever, confining ourselves within the frame of the familiar field theory we 
shall formulate a theory of multiple meson production in a way to clarify 
the field theoretical meaning of the FERMI (°) and LANDAU (+) theories pro- 
posed rather intuitively. 

Recently many authors have treated the multiple meson production in dif- 
ferent ways, but their theories may be classified roughly in two types; 4) theories 


On leave of absence from the Yokohama National University, Japan. 
W. Hetsenpere: Zeits. f. Phys., 101, 533 (1936); 113, 61 (1939). 
G. Waragnin: Nuovo Cimento, 5, 689 (1957). 

E. Fermi: Progr. Theor. Phys., 5, 570 (1950); Phys. Rev., 81, 683 (1951). 
L. D. Lanpau: Izv. Acad. Nauk USSR, 17, 51 (1953). 
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which are formulated in the frame of the field theory, such as those of HEI- 
SENBERG, LEWIS-OPPENHEIMER-WOUTHUYSEN (°) (briefly written hereafter as 
LOW), etc., 6) theories which seem to be outside of the usual field theory, 
such as those of FERMI (*7) and LANDAU (*8). 

Now, before entering into the detailed disscussion of the characters of 
these theories we shall see first which theories remain valid to explain the 
recent experimental data in an extremely high energy region such as the 
multiplicity, the angular distribution or the production branching ratio ete. 


a) Multiplicity. Heisenberg’s non linear theory (°) fails to explain the 
energy dependence of the multiplicity of produced mesons. For the multi- 
plicity n_ = 20 measured at the energy 10+'*eV the theoretical predictions 
are given as (1°) 


n> = 30. (LANDAU) 


<n,» = 11,7 (FERMI) 


b) Angular distribution. The angular distribution of air showers in 
cosmic rays was not explained by the LOW type theory. The values calcu- 
lated of the half angle 0,, inside of which one half of the energy is emitted, 
are as following for the energy of 10+14 eV, 


0 = 6107 (PERMT) 


0, = 1,3-10-4 (LANDAU) 


The experimental value being 1-10-', FERMI’s theory seems to give a rather 
large value. 


c) Lateral component. According to NISHIMURA (!1) the lateral component 
k, of the meson’s momentum in an air shower is almost always of the order of 
several hundred MeV independently of the incident energy. This fact is rather 


(°) H. Lewis, J. OPPENHEIMER and S. WouTHUYSEN: Phys. Rev., 73, 127 (1948). 
(5) H. W. Lewis: Rev. Mod. Phys., 24, 241 (1952); H. Fuxupa and G. TAKEDA: 
Progr. Theor. Phys., 5, 957 (1950); H. UMEZAWA, Y. TAKAHASHI and S. KAMEFUCHI: 
Phys. Rev., 85, 505 (1952); T. D. Lup and E. M. HENnLEY: Phys. Rev., 101, 1536 (1956). 
(7) J. Lepore and R. Sruart: Phys. Rev., 94, 1724 (1954); R. H. MILBURN: 
Rev. Mod. Phys., 27, 1 (1954); J. Kovacs: Phys. Rev., 101, 397 (1956). 
(8) M. HamaGucHI: Nuovo Cimento, 4, 1242 (1956). 
(9°) W. HEISENBERG: Zeits. f. Phys., 188, 65 (1952). 
(°) L. RosEntHAL and D: TCHERNAVSKIJ: Usp: Fiz. Nauk USSR, 52, 185 (1954). 
(1) J. NISHIMURA: private communication. 
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difficult to be interpreted by the FERMI theory due to his isotropical equili- 
brium state. 


d) Branching ratio of production. The branching ratio ny./m_. of 0°- and n°- 
production was measured recently by FUIJMOTO (1?) (~ 0,5). In the FERMI 
theory this ratio is determined only by the statistical weight of 0° and 2°, and 
results to be /n,.—= 2. On the other hand Kosa (!) showed that the 
LANDAU theory could yield the value 0,5. 


As seen above, the theory of LANDAU seems to represent well, at least in 
the present stage, the extremely high energy events. So, it will be rather 
important to analyse the essential characters of the LANDAU theory which 
gives us an ample scope to see the nature of these phenomena. 

First, stimulated by Heisenberg’s turbulence theory (!') FERMI had de- 
veloped the theory of multiple production as an evaporation of a perfect 
gas from the ‘hot spot’ in thermal equilibrium created in a contracted volume 
of meson cloud around the colliding nucleons. Then, taking into account 
the meson-meson interaction in the contracted volume, LANDAU regarded the 
production process as a jet of ideal fluid described by the 4-dimensional velo- 
cities u,(7t) ... u,(rt), energy density e(7t) and pressure p(rt). The behaviour 
of the meson field is determined by the hydrodynamical equation 


(1.1) 


Dy == (e+ pyu wy, + On: p 


and further, the STEFAN-BOLTZMANN relation is supposed a priori, as a state 
equation 


(1.3) Pag 


whe 
isp) 


As such, a theory of FERMI-LANDAU type appears as a foreign element in 
the usual field theory, but nevertheless it describes well the high energy phe- 
nomena. If so, in the ultra high energy region other unknown laws are 
dominant differently from the low energy region? Or, the present field theory 
may be still valid even in the high energy region but its solutions behave 


(12) Y. FUJIMOTO: Private communication. 
(3) Z. Kosa: Progr. Theor. Phys., 15, 461 (1956). 
(4) W. HEISENBERG: Zeits. f. Phys., 124, 628 (1948). 
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approximately as something like a liquid subtance? In order to obtain 
some information on the limits of the usual field theory it will be important 
to examine in which way and in which place the FERMI-LANDAU type descrip- 
tions can be situated in the frame of the present field theory. 

In this paper, the descriptions of FERMI-LANDAU type theories are derived 
in a field theoretical manner. Such interpretation may justify, more or less 
from the field theoretical point of view, the descriptions of FERMI and LANDAU, 
but not necessarily their artificial pictures such as « thermoequilibrium » or 
«fluid ». In fact, one can obtain the same features by a more natural model 
in a statistical treatment of the field theory. In Sect. 2 the general structure 
of the S-matrices for multiple production are discussed from the statistical 
point of view. Then, as a special case one obtains in Sect. 3 the perfect gas 
model (FERMI type theory) where the correspondence with the original FERMI 
theory becomes clear. Sect. 4 is referred briefly to the case of the imperfect 
gas model. The hydrodynamical behaviour of the meson field is derived in 
Sect. 5, in which the starting equations of the LANDAU theory are obtained. 
In Sect. 6 some comments on the LANDAU theory are given concerning his 
initial condition which plays an essential role in determining the characters 
of the meson fluid. Finally in Sect. 7 a simpler model as in Sect. 2, based on 
a Statistical field theory, is shown to give good agreement with the recent 
experiments. 


2. — Statistical treatment of the S-matrix. 


The statistical theory of meson multiple production firstly suggested by 
HEISENBERG and realized by FERMI has received many criticisms and objections 
because of its theoretical inconsistencies (!), but it plays an important role 
in the research of the general features of the high energy events. Then the 
concept of «thermoequilibrium » of FERMI’s theory being rather a foreign 
element in the field theory, some trials were done recently to derive the 
equilibrium state by treating the S-matrix statistically (1°). 

Now, to understand this idea clearly we shall consider an example of black 
body radiation. Taking a large hole surrounded by a perfectly reflecting wall 


(15) S. TAKAGI: Progr. Theor. Phys., 7, 123 (1952); H. Fuxupa: Phys. Rev., 89, 
842 (1951); C. AuLuck and D. KorTHARI: Phys. Rev., 90, 1002 (1953); S. NANDAU: 
Progr. Theor. Phys., 11, 605 (1954); W. KrausHar and L. MARKS: Phys. Rev., 93, 
326 (1954); K. Ono and K. Yoxor: Progr. Theor. Phys., 13, 101 (1955). 

(!9) D. Ivo: The essential idea of the statistical treatment of the field theory for 
meson multiple production was published in Japanese in 1955 and 1956. (See the foot- 
note of the paper of Y. TAKAHASHI) (16). Publication in English has been retarded 
by some trivial reasons. Y. TAKAHASHI: Nuovo Cimento, 4, 531 (1956). 
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we suppose firstly that in this hole region there exist only monochromatic 
photon waves of frequency ©. This state is represented by an occupation 
number of photons as n(@) n(o = 0) = %. Then, at the time? = h if we 
throw into this hole an atom which interacts weakly with photons, tran- 
sitions begin to occur between the different states of photons. So, after a long 
lapse of time the transition probabilities to the state whose occupation number 
distribution is n(0) = f(@) should be estimated by the square of the cor- 
responding S-matrix 


(2.1) Prob [n(m) = f(0)] ~| <y[m(@) = f(0)]; Spoln(o) = CE 


Experimentally however, we know that after a sufficiently long time the 
distribution of photon numbers will be following Planck’s law, 


i 
(2.2) Mo) = Sa (ico kT] xi GT. hoy): 
This means that the probability has a large value only for the state in which 
n(w) is given by Planck’s distribution. In other words, the pratically im- 
portant state is only y[n(w) =Planck’s distribution] (*). This state should be 
called as the most probable final state for radiation. 

In the field theory we could obtain, at least in principle, the transition 
probabilities to any state, apart from the practical difficulties of calculation. 
But in the real cases the created mesons are so many that the observed 
final state may be very often in the most probable state or in nearing ones. So, 
taking only the most probable final state the theory will be remarkably sim- 
plified. With such approximation we shall discuss the structure of the S-matrix 
for multiple production. 

The cross section of n pion production by nucleon-nucleon collision is pro- 
ved to be 

27): [1 
(2.3) On = | | n 


Vy 


(pip? |! | Pips; ko 


dpr toe 
1 WS ae Lan 


" §4(Ap — ky — ky — «.» — ka), 


(*) W[n(w) = 1/(exp [ha,/kT] — 1)] is one of the eigenstates of the Hamiltonian IT, 
for the hole radiation system, and satisfies the equation 


noliog = Hy 
has 


> exp [p(ha,/kT)| — 1 


8 


pe nov . 


ee ee Di AON = 
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where 
P}, P}: momenta of the incident nucleons 
V,: relative velocity of the incident nucleons 
P, =(P,, E;), P.= (P,, E.): energy momentum of final nucleons 


k, = (k;, ¢;):energy momentum of emitted mesons 


Pepe Pie Peeper XP 


(1/vn!) <|Z|>: invariant collision matrix of Moller 


Then, the most probable final state is obtained as a solution of the variational 
problem for the transition probability |<|/Z)>|?. For example in the occupa- 
tion number representation of the transition probability P[7,, n. ...], intro- 
ducing the entropy S by P[#,, no, ...;] = exp[S[nin ...j:d] (conservation), the 
variational problem is formulated as 


S| ORA 
Y Sn, — I RC 
On; 


> ta Ne, > ete = WV 


where 7,,”., ... are the numbers of mesons of momentum k,,k,, ... respectively, 
and N is the total number of emitted mesons and W is the total energy. 

But to see clearly the relation between the S-matrix and the most probable 
State we will transform the S-matrix according to the method of URSEL in the 
statical mechanics of imperfect gases (1°). Putting firstly |{P°P}|I|P,P.)|°= 
= C(P,P,) we define a new function f,(k) by |{P°P?|I|P,P,,k)|° = 
= C(P,P.)f,(k); and further f,(k) is defined by |<(P{P?|I|P,P,, kk? = 
= C(P,P.)[f.(K)}f.(k») + f.(k,k.)]}. Repeating such procedure the S-matrix is 
rewritten as 


=0, 


\(PYP!|1|P,P,> |? = C(P,P,) 
KE Pls k) da (P,P.)f,(k) 


(2.4) ) |<PIP}|I|P,P,, k,k,> |? = C(P,P.)[f.(k)f(k) + f.(k,k)] 


<P.P,|I|P,P,, k,k,k,) |° = O(P,P))[f()f.(Ky)f.(ks) + f.(k,k,)f.(k,) + 
+ fa(Fo)fa(Koks) + fi (k)fs(kxks) + fs(kxkoks)]. 


(1?) See for example HixL: Statistical Mechanics (1956). 
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In this expression the S-matrix is characterized by a series of functions fi(k), 
f.(k, ky), ... ete. Though the general form is rather complicated (*), if f.(K1ks), 
fs(k,kxks), ... are negligible compared to f,(k) the S-matrix is factorized in res- 
pect to the variables of emitted mesons, and the problem becomes much sim- 
pler. This is the case of the LOW theory though treated rather classically. For 
example, in the theory of FUKUDA-TAKEDA (°) with the PS(PS) coupling, the 
S-matrix was factorized and their f,(k) was (see Appendix B) 


(2.5) f,(k) ~(4), 


g, coupling costant, w, meson Mass. 
From such simplified model, one can see easily that fi(k) is a quantity 
intimately related to the form of the meson field (meson cloud) around the 


nucleon. 
In the case where the S-matrix is factorized there is no correlation between 


(*) Assigning the final state by the occupation numbers n; of mesons k; and by 
the «cluster » numbers n6_, of the correlated functions f,(kjh;, ...; k,), the cross section 
becomes 


Qn 1 g? Ni 1 1 g? nd 
dai Dall al’ RUI al ma (4) fk) 2,2,4,4) 


pats: ij May 


1 Maga 
a [i (4) ITA OO A; dl 


e (1) 2 
ijk...1 Vi jk...l U! \w 


-d(W — DI eM; — Di (go) 
i ij 


where A, = Vdk;/(27, 2; = (9?/1)(u/e), V normalized volume, 4 meson mass and 


g coupling constant. The problem is solved by a variational method, giving, for example, 
the «cluster » number distribution in the most probable state, 


n*(Key ki) = (kx bey) 21.Q. + Q exP [— Bley te +--+ 21), 
where the temperature 1/6 of the system is determined by the conservation 
Gi 1 
wW= Za dk... dky@y%-- Qifilley «.. ki) exp [— Ble: + - + 801. 
And the multiplicity is expressed finally as 
1 
<n» = Sri dk, oe dk, 2, tee Qf (ky see k,) exp [ Bley + eee + e1)] . 
Ll . 
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the created mesons, and the mode of production may be something like a jet 
of an ideal gas. In fact, as will be described later, FERMI theory just corre- 
sponds to an evaporation of this ideal meson gas from the « hot spot » of the 
colliding region. On the other hand, where the functions f(k1ks), fs(Aik.k,) ete. 
are not negligible there are correlations between created mesons or mesons 
and nucleons. In such cases the production may resemble a jet of an imper- 
fect gas, and if the meson-meson interaction is very strong the meson cloud. 
may behave like a fluid. 


3. — Perfect gas model (Fermi type theory). 


In this case, the S-matrix being factorized, the total cross section takes a 
simple form 


2 5 d 1 2 1( 1 2) see J1 n 
(3.1) on = 92) feto Pr dps f(x) (ks) -- fa (ken). 


son E, n! 
dk dk, dk, 


SURE, En 


SAD oe ee ee 


and in the occupation number representation it becomes 


(2a0)4 ae je 
a > ¢(pips) IT [fh l'a — ¥ nk) 
(3.2) Ak, = way dk 
8 (27)? s 


Then, the most probable final state is given, as was described above, by a solu- 
tion of the variational problem, 


| San. = a l= 0, 


(3.3) m,! &, 
| D ee = 0. 


When n, is sufficiently large, n, is reduced to 


s 


(3.4) n=) n 


where the unknown constants (Lagrange multiplier) B, are determined by the 
conservation relations 


(3.5) (Ap), = > (,),n 
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Or, in momentum representation 


h(k) 


SO gres 


This is our starting general equation obtained by a statistical approximation 
of the field. The detailed comparison with the experiments by this equation 
will be given in Sect. 7. Here, the special form of fi(k) is discussed con- 
cerning the physical meaning of the different statistical theories of production. 

To simplify the problem, neglecting the momentum conservation the above 
results are expressed as 


k 
| n(k) = fi(k) exp [— Bei] =, 
(3.6) : i 

| N | f,(k) exp [— Bev]ak , 


where AW is a difference of energy between the initial and final states of the 
nucleons. If we, as LOW, take for the meson field the form f,(k) = g?//’, 
the distribution n(k) becomes, 


| Mk) = i — vzexp[— fex]dk, 
(3.7) O 
| = 


to) 


. 


This is just the result of LOW. On the other hand, if the eigen field function 
fi(k) is taken as a contracted volume of FERMI 


(3.8) O 


we get justly the result of FERMI. 


| mk) = 2 /1 — vrexp[— Bex] dk, 
(3.9) | er ila 


Now we will remember the meaning of the Fermi theory from a general 
point. of view. The Fermi theory is based on the intuitive supposition that from 
a hot spot of thermoequilibrium state in the collision region mesons are emit- 
ted as by an evaporation. Owing to the high velocity of the incident 
nucleons, the contracted volume of the hot spot is 


= ACNE ib 
(3.10) Lee = Q VAI = vw V1 = do 3 IT (1) . 
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where v, is the nucleon velocity. With a thermoequilibrium hypothesis, the 
energy distribution becomes Planck’s one 


Qrak 
n exp [Bex] SUE 


(8.11) n(k) dk 


where the temperature 1/8 of the system is determined by the conservation 
of energy. 


(3.12) Jenuyas = Vey 


where 4, is the incident kinetic energy of the nucleons in C.M. system. Then, 
the multiplicity is given by 


12 is Ey 3 Ei i 
(3. ) : (n) = N “i , 


in good accordance with the experiments. 

As was seen above (Eq. (3.9)), FERMI results are obtained also by an ap- 
proximate method in the field theory. It should be noted however, that what 
was derived are the results of FERMI, but not his picture of thermoequilibrium 
state. In a field theoretical way, the problem of meson multiple production 
is to calculate any probability 


P = P(njn, 3) 


for production of n, mesons of wave number k,, n, mesons of wave number 
k, and so on. In practice however, it is too complicated because of the large 
number of produced mesons, and even a calculation of the phase volume is 
almost impossible. But for an exceptional process one can calculate the pro- 
bability. This is the case where the eigenfield of nucleons is shaked off at 
random by the shock of collision, and the probability becomes, as is easily 
seen (1°), of the Poisson type, 


(3.14) Pp~ [teen Hayes 


f(k) being a probability of finding k mesons in the eigenfield. This form of 
f(k) is, of course, dependent on the Lorentz contraction of the eigenfield. In 
general, the process of shaking off k, mesons after shaking off k, mesons 


(18) F. BLocH and A. NORDSIECK: Phys. Rev., 52, 54 (1937). 
(9) W. PAULI and M. FIERZ: Nuovo Cimento, 15, 1 (1938). 
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is different from the shaking off in a inverse way, so the supposition « at random » 
is not exact. Even so, this « at random» may be a rather good approxim- 
ation and forms the standing point of LOW. In such way, the f(k) was 
calculated and used by LOW in a form f(k) ~ (g?/u*)(u/éx)- 

Accordingly noting the similar form and same order of two factors, 
(g?/u)\V1 — v2 in LOW and (4/3)x(1/u*)V1—vì in Fermi (eqs. (3.7) and 
(3.9)), FERMI theory may be considered as an intuitive one for the statistical 
treatment of the field theory. But our most probable state method (hereafter 
called as Fermi-type theory) contains, in the physical contents, much more 
than the thermoequilibrium state, and is more rich in the physical interpre- 
tations. This point will be seen in detail in Scct. 7. An example to clarify the 
meaning of Fermi theory more concretely is given in Appendix B. 


4. — Imperfect meson gas. 


In Sect. 2 the structure of the S-matrix was discussed by the Ursel operat- 
ion, and in Sect. 3 the special cases were considered where the S-matrix is 
factorizable. When the S-matrix cannot be factorized due to the strong 
correlation between emitted mesons, one may treat it as an imperfect gas. 
But though, as will be seen in Sect. 7, the perfect gas model (Fermi-type theory) 
is sufficiently good, at least in the present stage, to explain the recent exper- 
imental data, we shall give here an imperfect gas model mathematically and 
abstractly. 

For this purpose a functional o[Z] is introduced as follows, where Z(k) 
is an arbitrary parametric function 


2a): (dp, dp, di; 
a) ofZ] = Soa] = IZ] 


- |Kp®pî|I|p1p:; ki... ka) |}S(Ap—Xki). 


As will be proved easily, this functional satisfies the following relations. 


1) gg = DI RR 

2) 5 P i azar. 
Sof Z 

3) soe = OCF) 


where o(k) is the cross section of production of the k-th meson. 


S2o[ Z| 


sz a(R) Oe Bs 


4) 
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where o(kk') is the cross section of production of the k-th and k’-th mesons. 
Or, setting o(Z) = exp[S(Z)] one obtains 


__ o(k) 

SAE ate oan 
CRY) 3S(Z) $912) 328(Z) of kk’) 
ae SZ(k) ee) Oia 


Introducing (2.4) into (4.1), 


(22)* (dp, dp, 
(CM VAS arf me (pips) TI Y 


Remembering that the 6 function is 604(v7) = (1/(27)*) [exp [— 7€-a]dé, c[Z] 
is expressed finally 


dp, dp, ; : 
of] = X al) =| PE elp.p.)| ass exp[—iApe + X (2, 8) 


n 


where 


(4.5) b{Z,é]= lS eae f, (Key... Ea) exp [i(k +... + a) €]-:Z(k) <Z(ka) « 


Ej 


Then, the most probable state of o[Z] is evaluated by the saddle point method 
and given 


dp,d : 
(46) of Z] =| pp (Pipa) exp [—Ap-B + YbZ,é=— id)], 


where 6 is determined by the conservation relation 


dk, dk, 
(4.7) Ap = 3} Sep VSS] 


E 


In such formulations, the branching ratio of production becomes 


o(k) x db[Z] 
Ttot a > dZ(k) yee 
oY = 1 kk,. CREA dk, dk 
=Zi=piffi k,) exp [B(k + si) | Sede 
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Consequently, if o(k), o(k, k’), ... ete. are obtained experimentally one may 
get a good information on the S-matrix structure, that is, the structure of 
the meson field. 


5. — Hydrodynamical behavior of the meson field. 


The fact that Fermi theory could not explain the distributions of angles 
and momenta in an air shower might show an incorrecteness of neglecting 
the correlations between the created mesons. Then, LANDAU gave a reasonable 
interpretation of these phenomena supposing that the once created mesons 
in the mode of FERMI are emitted like a jet of fluid, in place of gas eva- 
poration, because of the strong interaction and radiation pressure between 
mesons in the contracted volume. But his arguments were always classical 
(not quantum-mechanical) though very relativistic. So, from a field theore- 
tical consideration we shall derive some hydrodynamical features of the meson 
field which were the starting point of the Landau theory. 

For simplicity, we consider an ensemble of the charged free mesons whose 
momenta are k,, k,,...k,. In such state Wk, ... k,) the following quantities 
are examined. 


(5.1) o(0) = CP*, D*(w) D(a) PY) , 
(5.2) J,(2) = (vs, D* (a) dad DL) ; 
ua 
oS) Ta) Oe @* (x) i + DW) 1 


The meaning of these quanties is obvious, o(v) density, J,,(7) current and 
T ,(#) energy momentum tensor. By the Fourier transformation 


(5.4) D(x) = [ow exp [tka] dk , 
they are rewritten 
(5213) o(~) =| dkdk' exp [— dk — k')ar\o (kk) ; 


(5.2') J (a) = | dk dk’ exp [— i(k — k')a] (43) o(kk’) , 


(53) Ta = fas dk’ exp [— i(k — k’)a] (3°) (3 ) o(kk') , 
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where 


(5.5) o(kk') = <P, D*(k)D(k'YL) . 


Or, for the physical interpretation, it will be useful to transform them as 
following 


(5.1”) 0(#) = favor. x) ; 
(5.2") J,{0) = | App ,o(P, 2); 
(5.3") Tuo(®) = farpupseto x) ’ 
where 
ie se ik ; 
(a sr k_-k'=1, o(kk') = o(p,1), 
(5.6) 


| o(P, x) = fe l) exp [tla] dl . 


As ¥Y contains only the free mesons, the matrix elements o(kk’)=<k|o|k’> 
have no vanishing values only when k? + w?=k'? + u?= 0, where uw is the 
meson mass. Then, J,(x) and 7,,(4) are proved to satisfy the continuity 
equation 


0 
(5.7) da, Jc 0, 
(5.8) Sa, 
n Tuo) =0, 


and also it is seen immediately 
(5.9) DT. (0) = — u<P*, O* (x) Da) Py 
u 
And, J (x) being the current density, the quantities u,, defined by 


(5.10) U,(v) = 


correspond to the velocity of fluid which fulfils the 4-dimensional velocity 
relation — w,u,=1. 
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As is known in the fluidmechanics the average value of particle velocities 
k,/e represents the fluid velocity, and the mean square of the difference of 
fluid velocity and particle velocity gives rise to the pressure or viscosity. 
To understand these characters it is convenient to transit to the coordinate 
system where the fluid velocity is zero, that is, in this new coordinate the 
velocity vector becomes 


(5.11) (e) = (0, 0, 0, 1) . 
This Lorentz trasformation is represented by 4,,; 


| w,(2) me a (0) (€) = Cy 
k= A yy ky 
| Ty(%) = Ty (@)a,ry 


SIIT 


Now, we suppose here that in the coordinate system flowing with the meson 
fluid the distribution of momentum of composite mesons should be isotropic. 
For a free field this may be a rather natural assumption. Then, k+-k'/2=p 
is also isotropic in this system, and the energy momentum tensor 7’, is tran- 
sformed as following 


T'uy(®) Far Op! T p'y(®) Ay’ ’ 


(5.13) È? tun [APP APO, L)Ayry ’ 


o Ign Spa ee 
= D duilvi 5 [appre + Uras APPEOLPD) 
t= 
By introducing the two quantities p(w) and e(x) (*) 


; | dep*elr, x) = p(x), 
(5.14) 
Janvier. x) = ela). 


(*) In the case where the momentum distribution of mesons is not isotropic, the 
definitions of p(#) and e(#) should be taken as 


È _ Se 
p(x) = +> T(x) , 
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T,, takes finally a familiar form 


(5.15) 
| = play(6,, + 4,4) + ela)u,u, 


where the orthogonality condition of Lorentz transformation is used. 

As is readily seen from the defining equations, p(r) and e(x) represent the 
pressure density and energy density respectively. Consequently, this T,, 
is just the energy momentum tensor of a perfect liquid proposed first by 
LANDAU. 

Between the pressure and energy there exists a relation 


(5.16) p(@) < 3e(2), 


because of 


= — pO OV <0. 


In the limiting case p(x) = e(x)/3 which is called STEFAN-BOLTZMANN relation 
realized in the thermoequilibrium state, this is just the state equation of 
LANDAU. Therefore, one sees that an ensemble of free mesons behaves, to a 
certain degree, as like a fluid obeying the hydrodynamical relations of LANDAU. 


6. — Comments on Landau’s theory. 


LANDAU considered the multiple production process in three steps: 1) The 
formation of local thermoequilibrium in a contracted meson cloud because 
of the strong interaction and high density of mesons. 2) A hydrodynamical li- 
near expansion of the meson cloud, due to the radiation pressure, in the direction 
normal to the contracted disk. 3) The three dimentional emission of free 
particles owing to the thin density of nuclear materials. Such description was 
proved to be very effective to explain the angular distribution of air showers 
that the Fermi theory failed to do. 

As was discussed in the previous section the meson field is described as 
a fluid if the meson velocity distribution is isotropic with respect to the fluid’s 
velocity. In the Landau theory where the particles arrive at the state of thermo- 
equilibrium as in the Fermi theory this condition of isotropy is satisfied. How- 
ever it may appear that this isotropy should contradict the linear expansion 
of the meson fluid of LANDAU. In fact however, the equilibrium region being 
contracted like a disk the pressure has a large gradient in the normal direc- 
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tion, say 2, of the disk, and it is this gradient that forces the linear expansion 
of the Landau fluid. The pressure is given by (see Eq. (5.14), 


; (vpi if 
(6.1) p(v) = [teo x) = IUZAL exp [il-w]o(p, |) . 


On the other hand, Fermi’s equilibrium volume will expand isotropically with 
the thermal velocity 


(6.2) cane Jo (2) o(p, »/ [et 2) ap] 


Then, in order to realize Landaw’s linear expansion the spurting velocity 
\l| le must be much larger than the thermal velocity v,. By definition 
of p and J we have a relation for the momentum of emitted mesons 


(6.3) k=pt 


So, the above condition of linear expansion becomes 


PRIORE 
(6.4) day = pd + > CED, <ED 


namely, at the instant of creation the « component of meson momentum must 
have values very much larger than the other components. 


YA 


| 


Fig. 1. Fig. 2. 


Accordingly, in order to understand Landau’s hydrodynamical theory as an 
approximation of the field theory in the high energy region, the momentum distri- 
bution of mesons at the initial time of creation should be pictured as in Ried: 
Or, in the configuration space, this may correspond to the contraction of the me- 


son cloud as in Fig. 2, where the uncertainty relation {Ax)<Ak,) ~ 1 exists. If 
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the multiple production process is a shaking off of the contracted proper field 
(as in Fig. 2), the initial condition of momentum (as in Fig. 1) should be rather 
reasonable. In fact, as will be seen in the next section, only by this initial 
condition of momentum one can obtain almost all the same results as that 
of LANDAU. Namely, in his theory only the initial condition is of vital impor- 
tance, and the complicated hydrodynamical treatment may not be essential. 
In other words, if only his initial condition is taken one can obtain equally 
good results by a more familar and usual method. Then, the success of 
Laudau’s theory does not necessarily mean that the meson field is degenerate 
like a fluid in the extremely high energy region. 

According to Kosa (*), in the Landau theory the essential features of 
the extremely high energy events are determined principally by the low energy 
processes in these phenomena. If so, the ultra high energy events may be 
not so valuable, in spite of our expectation, to obtain useful information 
on the interactions of elementary particles in the high energy region. However, 
the bases of the hypothesis are less clear and one should be more careful to 
give them up, and it should be noted that the characters of the Fermi theory 
have not yet been clarified. The experimental basis of the conclusion that 
the ultra high energy characters are rather dependent on the low energy 
process is the branching ratio of 0° and x° production. Being impossible to 
explain it by the Fermi theory, KoBA had recourse to the strong final inter- 
action of the Landau theory. But there exists also another interpretation 
of this branching ratio basing only on a creation mechanism in high energy 
processes without any strong final interaction (fluid approximation) in low 
energy processes. In fact, recurring to the Fermi theory as an approximation 
of the field theory (Fermi type theory), we shall show in the next section that 
one can get, in a rather general way, almost all the same results as those 
of LANDAU. 


7. — A simpler model (Fermi type theory) 


FERMI derived his results basing on the picture that in a small contracted 
volume of meson cloud the thermoequilibrium states are attained immediately 
after the high energy collision of nucleons. Then, the velocity distribution 
of particles in equilibrium state is necessarily isotropic, and in order to ex- 
plain the sharply forward collimated angular distribution he had to resort 
to the restriction of phase volume by the conservation of angular momentum 
in collision. In this case the details of the interaction matrix elements are of no 
more importance. 

As described roughly in Sect. 3 a theory very similar to Fermi’s but much 
richer in physical contents was obtained by treating the field theory with 
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analogy to the statistical mechanics. Namely, in the case that many mesons 
are produced, ignoring the detailed knowledge of the probability one has only 
to consider the most probable state attained by the produced particles. 

In the simplest case where the meson fields are shaked off at random 
by a shock of collision the most probable distribution becomes 


a AN 
~ exp [Bex,] —1° 


(T.1) n(k) 


depending on the functions f(k) related to the form of the eigenfield of the nu- 
cleons. As was seen in the case of LOW these f(k) belong to the spectrum 
of the eigenfield of the nucleons. But we have no detailed knowledge of the 
eigenfield of the nucleons, especially in the extremely high energy regions. 
Rather, we should obtain some information on the eigenfield by the analysis 
of the high energy events. Then, we shall suppose only the following 
qualitative characters of the eigenfield. 


1) The extension of the eigenfield of a rest nucleon is of the order of 
the Compton wave length of the meson, that is 
h 


(7.2) Ax) = <Ay> = (Az) ~ 
pe 


‘Then, for a nucleon moving in the X-direction, (see Fig. 2) 


Ey\ 1 
Ax) ~ Ge 


1 (ia oe == 1)5 


(7.3) 


where E, means the energy of colliding nucleons in the center of mass system. 


2) In the collision of nucleons the mesons are produced as shaking off 
of the eigenfield. For the wave-packet of mesons there exist the uncertainty 


relations (see Fig. 1) 
(7.4) (Aa) <Aky) > 1. 


3) Though our distribution n(k) contains the eigenfield spectrum f(k), 
the multiplicity predicted coincides with that of FERMI 


E} 
(7.5) (Ma) =V (3) ’ 


y, costant. 
Such hypotheses appear to restrict too much the spectrum f(k), especially 


the hypothesis 3) might seem to be impossible. In fact however, these assum p- 
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tions do not impose severe restrictions, and as will be shown in Appendix A 
a simple step function 


— (Hy /M)u < ky < (EM), 


ik) ife reese 
(7.6) 
—u<k.<w, 


= otherwise 


satisfies these conditions. The physical picture of this function f(k) is evident. 

In the original Fermi theory f(k) comes from the matrix element and is con- 
sidered as a constant factor which plays a role of contraction of the equilibrium 
volume. In our case this f(k) is responsible for the angular distribution of 
emitted particles. FERMI supposes a realization of thermoequilibrium state, 
while we think the production process as a shaking off of the eigenfield in the 
form of wave packet. This difference appears also in the branching ratio 6° 
and 7° production. 

By the hypothesis 3) the meson kinetic energy per meson emitted, i.e., the 
temperature 1/6 is given by the conservation of energy. 


— 
ESì 
ni 

— 


IENE 2 Po 
ese ee 


where v is an unknown constant. But from the case of Shein’s star (20 emit- 
ted mesons at a laboratory energy of 1012 eV), one can obtain v= 2.4. One sees 
that 1/6 is sufficiently larger than the pion mass w. Then, the average value 
of lateral component of meson momentum can be approximated as 


i _ Sie |n(kydk fk, |f(k)dk 
(7.8) Lae I “Tak | 


because k, integral is now mostly dependent on f(k). Consequently, 
(7.9) (kh > > <Ak,> ~ w = 140 MeV/c , 


which just explains the experimental results of NISHIMURA for the lateral com- 
ponent of meson momentum. 

Next let us consider the angular distribution of emitted particles. In our 
method the averaged angle that the particles make with the air shower axis 
is given simply (in C.M. system) 


} LI pl 
(LO GSS O, ; 
os Chad <kad 


THEORY OF MESON MULTIPLE PRODUCTION ih 


On the other hand, due to our assumptions 1) and 2) ky > ky = k, (this is just 
the initial condition of LANDAU), k, is estimated as 


E [k.,m(k)dk fen(k)dk 2H, _2(E } 
an ieee dk femk)dk 25 _ 2 (Mo 
ot he {n(k) dk Na) Na (a) Aah 


Accordingly, the average angle in the C.M. system becomes 


3 vu (Ho\* 
0° MEN e ee eta $ 
3 OM Ga 


Transformed to the laboratory system this angle is expressed as 


co, — vu M\} vu 2M \î 
4 4M Eo 4M Era i 


E ide t (ia M gem. 
—_— ~~ OM , ~~ a Da 4 
are used. 


Now taking the empirical value y= 2.4 based on Shein’s star, the final 


where the relations 


expression for 0"" is 


COS 0.15-( A 
i , ESA 


For Fermi’s and Laudau’s theories, TCHELFSKY and ROSENTHAL had given the 
half angle 0, in which one half of the energy is emitted, by calculating the 
energy flow eu , (€ 18 the energy density, «, the velocity). Our results are com- 
pared, in Table I, with those given by them. As is seen, our results are 
essentially the same as that of LANDAU. 


TABLE I. 
Era 1013 eV 1014 eV 1015 eV 
Qt (Ours) ise NO Do Ome 4.8-10-§ 
04 (LANDAU) Tease WO poi Ome 3.8: 10-° 
0, (FERMI) 6.7-10-4 Dale ies Gail Ome 


Finally we consider the branching ratio for 0° and n° production. As the 
simplest example, we take the spectrum f(k) as that of LOW, i.e., f(k) = 
= (g?/u?)(u/e,). Then the multiplicity of 7° becomes, as was described in Sect. 3, 


2 
(7.16) ae De di 
i br exp [Bex] —l & 
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and in the same way for 6° mesons, 


(7.17) Cnos nf ; ce 
Hs Ne> = — | — 
ee ug] exp [Beo] —1 & ’ 


where gy, g. are coupling constants, ug, “, the masses of 0° and 7° mesons re- 
spectively. The common temperature is determined, as previously written, 
by the conservation of energy. This temperature being very high the integral 
part of the multiplicity equations are almost all the same for 6° and 7°, and 
the branching ratio is given 


Ng) (2) (2) (2) 
7.18 ~2|—)\|—)~ 0.16 |—). 
( ) (Nn) Mo Gx Ir 


The ratio of the coupling constant of production (g,/g_)? is not known well, 
but from the data in the low energy regions one may regard it as the order 
of one. Then, the branching ratio is 


(7.19) Mao 


which is not so large as that of FERMI. So, it will be possible to get the cor- 
rect value of the ratio only by the creation mechanism without the complicated 
final interaction in low energy processes. 

In the original Fermi theory, owing to the picture of thermoequilibrium, 
the volume in which 0 and 7 mesons are confiend must be the same. While, 
in our picture where the mesons are shaked off as packets the extension of 
the packets may be different for 0° and 7°. This may correspond, in the ori- 
ginal Fermi theory, to the different equilibrium volume for 0 and x. 


Gn lin 
Ere ) 
Mx En 
(7.20) 
ehe 
Mo £0 


If the ratio of the extension of the packets (Go/to)? (Gul tn)? is equal to 1, one 
obtains of course the Fermi result (n,/n_) =2. 


8. — Conclusion. 


1) Fermi’s and Laudau’s theories of multiple production seem to be 
foreign elements in the field theory, but they can be formulated in the frame 
of the present field theory. 
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2) Landau’s theory is at present regarded as the best, however, one can 


obtain equally good results by a statistical treatment of the field theory with 
a much simpler model. 


APPENDIX A 


As the simplest example of the spectrum function f(k) which satisfies the 
hypotheses 1), 2) and 3) in Sect. 7, we take the following function 


| — (H,/M)u<k, < Hy/M, 

(A.1) f(key Ra he) = gio if = RT SI 
— <k,< +p, 

=0 otherwise. 


For simplicity we shall consider the Boltzman distribution in place of 
the Planck’s one. Then, the multiplicity becomes 


(E/M)p + + Ml 


ME 2 fan. fax, fax. exp [— B(kî + ki + kî4 TOS 
ee —p =p 

È (By!M) 10 
(A.2) Oe iat 

—(Eo/M) 

pees gy? 1 exp eS (E,/M)pu] 
Bu 
and 
Kn) (* — exp[— (E,/M)Bu] exp [— (E/M)ba] i 

AGS E = TI ——_-- ™ 2 ze < È 
SIRIA Bu Bu U 


Putting pul,/M = y, (A.3) is rewritten 


[2 bay 
gg 


M? 1-exp[ yl 


y 
—— SI, = ex — | ~~ 
Pu, 7 ae oa 


(A.4) 


But in the region 0 < y <1 (A.4) has no solution, because in the high energy 
case M?/g?uH, < 4. So, the unique solution of (A.4) is 


pi ae 


mM ’ 
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1.04 


ni 1 (E,\* 

(A.5) — =- (=) ; 
Inserting (A.5) into (A.2) 

M\}(E,\} 

wma) (ae) 


Therefore, the packet (A.1) satisfies the hypothesis 3). It is self-evident that 
(A.1) satisfies the hypotheses 1) and 2). 


APPENDIX B 


In order to see the relation between the Fermi theory and the statistical 
field theory, we shall give an example starting from the Ps(Pv) interaction 
Hamiltonian, H(z), 

Op(@) 
cn. 


(B.1) H(a) =i . B(wyysy (0) 


LE 
Then, the cross section of n pion production is, as Fukuda-Takeda (5) (*), 


2 


(B.2) 0, 


__ 2nM (1 M\|VE(y:w)VE 
— vob \u a M 


where v = po/E is the velocity of nucleons in the C.M. system, , its energy. 
(Ep), (e;k;) are energy momentum of final nucleon and of pions respectively. 
Now, this cross section is transformed into the representation of meson occup- 
ation numbers by enclosing the whole system in a large volume V. Denoting 
the phase volume Vdk;/(2x)? by A,, it turns out to be 


È 22M (1 M 
(B.3) On = ii? (E Ir, Vier. Ve 
(Ny +Nq +... =N) 
where 
V E(uyUy) VE, |? 
D(Ny No ’ Pee) = | ( Li 9 d(w ny = 200) 
(m+ Ae — 1)! (A) 
(B.4) py ie tA 7 
ian V fee 


= exp [S(m, my ...)]8(w — > evn) , 


— 
* 
— 


The factorization with respect to meson variables is done as an approximation. 
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With the above preparations, the occupation numbers specifying the most 
probable final state are obtained by solving the variational problem 


OS (N1, Na «.-) 


B.6 Sn, 
ded DOT dn; 


=" 0 
with the restrictions 

BIST 

1 

|Den=w. 


For large occupation numbers, the problem is readily solved, resulting 


(B.7) 


dk, 
SL SALTA 


(B.8) n= 


where « and f are Lagrangian multipliers determined by the conservation 
relations. 

For extremely high energy collisions, if no restriction is imposed on the 
total number the «chemical potential» x may be ignored. Then, the distri- 
bution of occupation numbers in the most final state is Planck’s one, 


Q,Ak 


(B.9) n(k) = AAC 


where the « temperature » 1/8 is determined by the energy conservation 


E Q,ex dk 
(B.10) w <a È 


Further, Q, given by (B.5) may be rewritten as 
(B.11) Ore VIDI, 
where v, represents the velocity of the emitted pion, and 
9° 
(B.12) Q) = DE ’ 


is, as is easily seen, a quantity with the dimension of a volume and its mag- 
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nitude is comparable with Fermi’s volume of thermoequilibrium 


eed 


All of these results derived from the quantized wave field by introducing 
statistical considerations are very similar to the results of FERMI. (B.11) is 
our counter part to Fermi’s assumption of Lorentz contraction of the equi- 
librium volume 


(B.14) OO a 


SUPPLEMENTO AL VOLUME VII, SERIE X N. Ic 1958 
DEL NUOVO CIMENTO 1° Trimestre 


Statistical Theory of Particle Multiple Production 
in High Energy Nucleon Collisions. 


V. S. BARASENKov, B. M. BARBASEV and E. G. BUBELEV 


The Joint Institute for Nuclear Research 
The Laboratory of Theoretical Physics - Doubna 


(ricevuto il 25 Luglio 1957) 


Contents. — 1. Introduction. — 2. The calculations of phase volumes. — 
3. The model of compound particles. — 4. The results of calculations. — 
5. Discussions. 


1. — Introduction. 


It has been shown in a series of papers by the late professor S. Z. BELENKIJ (4) 
and his co-workers that in the energy region up to 3-5 GeV the probabilities 
of pion production in pion and nucleon collisions with nucleons, caclulated 
using Fermi’s statistical theory may be brought into agreement with experiment 
only by taking into account strong resonance pion-nucleon interaction in the 
state p= 3, T= 3. This idea is equivalent to the introduction of a pheno- 
menological particle N’ («isobar») with spin P= 3, isotopic spin TS 
and mass M,=1.32M,, where M, is the mass of a nucleon. It is reason- 
able to keep the phenomenological particle N’ also in the calculations made 
in the higher energy region of the colliding particles. At energies > 10 GeV, 
when intensive production of strange particles is possible besides the usually 


(1) S. Z. Beren'K1J and A. I. NIKISEV: Zu. Bksper. Teor. Fiz., 28, 744 (1955); 
A. I. Nixigev: Zu. Éksper. Teor. Fiz., 30, 601, 990 (1956); V. M. Maxstmenko and 
A. I. Nikisev: Zu. Éksper. Teor. Fiz., 31, 727 (1956); S. 7. BELEN'KIJ: Nucl. Phys., 


2, 259 (1956). 
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used conservation of energy-momentum, isotopic spin and its projection it 
is necessary to take into account the conservation of strangeness and baryon 
number. 

At present in the frame of statistical theory of multiple production one 
fails to take into account momentum conservation. An attempt in this di- 
rection made by FERMI (?) as it has been shown in the papers (*) is not correct. 
Therefore, when making the calculations we neglected momentum conservation 
but assumed the statistical weight of a set of particles produced in the 
particle collision with spins S,, S,,..., S, to be proportional to the product 
TI 28,4 1). Since the spins of strange particles are not exactly known at 
ent) for the estimation of the value of the spin the minimum of all possible 
spins has been chosen: S= 0 for K-mesons and S = 3 for all the other strange 
particles. 

In the papers (‘) we have shown that one may hope to bring into agreement 
the results of the calculations on statistical theory with experiment only by 
taking into account the internal structure of the colliding nucleons, what from 
the mathematical standpoint is equivalent to the introduction of the additional 
parameter into the theory. In our theory the value of this parameter is chosen 
from the requirement of bringing the calculated ratio of the number of pro- 
duced strange particles to the number of produced pions and nucleons into 
agreement with experiment. 

It is of great interest to consider this questions at the maximum energies of 
10 GeV achieved at present in an artificial way. 


2. — The calculations of phase volumes. 


In the mentioned paper the respective probabilities of pion, nucleon, anti- 
nucleon and strange particles production in the nucleon collisions with 
energy 7 and 10 GeV have been calculated using statistical theory. Phase 
volumes in the momentum space in all cases of two and three particle pro- 
duction have been calculated using the exact Block’s formulae (5) and in the 
cases in which a great number of particles were produced, using the 


E. FERMI: Phys. Rev., 81, 683 (1951). 
E. L. Fesnperc and D. S. CerNAVSKIJ: Dokl. Akad. Nauk SSSR, 81, 795 
; I. L. RozentAL' and D. S. CERNAVSEIJ: Usp. Fiz. Nauk, 52, 185 (1954). 
V. S. BARASENKOv, B. M. BarBaSEvV, E. G. BuBELEV and V. M. MAKSI- 
MENKO: Zu. Eksper. Teor. Fiz. (in press); Nucl. Phys. (in press). 

(9) M. M. Brock: Phys. Rev., 101, 796 (1956). 


STATISTICAL THEORY OF PARTICLE MULTIPLE PRODUCTION ETC. 119 


formula: 


(1) W mn(&) = done dem DS) Dy eats II M}- 
t=1 


i=1 
gqen+&m_-3 


Tn + 8(m—-1)) 


F(2n; 3; 3n + 3(m—1);— a), 


where m and n are respectively the number of non-relativistic and estremely 
relativistic particles 


E is the energy of the colliding particles in the center of mass system; 
M; are the masses of non-relativistic particles ; 
Fis a hypergeometrical function. 


Formula (1) is the consequence of the expressions for phase volumes ob- 
tained in the papers (9). 

For small values of the parameter « (practically « < 0.70.8) the hyper- 
geometrical function FP may be presented in the form of a series expansion in 
powers of « which is quickly convergent and convenient for the calculations: 


(2) F(2n; $; 3n + $(m—1); —@) = 


I(3n + 3(m—1)) 2 a (—1)'P(h + $) 7 
= = RE 2 On FG l 3n 3(m 1)) x >) 


where C%, ,,, are binomial coefficients. 
For the numerical calculations in the case of greater values of the para- 
meter « the function F is conveniently written in the integral form: 


(3) F(2n; 3; 3n + 3(m —1);—@) = 
l'(8n + fm 1) fe —oneinst 
— T(n + Hm —1) Ten) | (1 + of)! 


At great number of the produced particles (> 3) the formulas (1)-(3) ap- 
proximate well enough the exact expression of phase momentum volume and 
in cases when the part of the produced particles is deliberately « intermediate », 


(6) I. L. RozentaL' and V. M. MAKSIMENKO: Zu. Bksper. Teor. Fiz. (in press). 
(7) 8. Z. BeLEN'K1J, I. L. RozenTAL', A. I. Nixisev and V. M. MAKSIMENKO: 


Usp. Fiz. Nauk (in press). 
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i.e. cannot be referred to the category of extremely relativistic and non-rela- 
tivistic particles. In the cases when the precision of the formulae (1)-(3) 
appeared to be unsufficient the interpolation methods were applied. 

As comparative calculations have shown, the arising error is insignificant (*). 


3. — The model of compound particles. 


As well as in the calculations made earlier (4) sharp difference with exper- 
iment (8°) arises if we assume all the particles to be produced in the same 
space volume V_ or V,, (respectively with radius equal to the Compton wave- 
length of the pion or the K-meson). The choice of other values of the volume 
V leads to analogous results. The contradiction with experiment may be 
removed if we assume the different kinds of elementary particles to be 
produced in different space volumes. It can be understood in view of the 
formation of indivisible « compound particles » by the colliding particles where 
all the energy of the particles is concentrated in the interaction energy, but 
not in the proper energy of the particles (1°). The dimensions of the regions 
of interaction energy localization are specific for every concrete type of 
interaction. 

We have made the calculations for two different assumptions in respect 
to the character of elementary particle interaction: 


1) If we assume that the hyperons interact with pions as strongly as 
the nucleons, then these particles must be produced in the same volume V.. 
Corresponding antiparticles as well as isobars and antiisobars will be produced 
in the same volume. For K~ and K*-meson we assume a considerably less strong 
interaction with the pion field, so that these particles will be produced in the 
volume V,,. The interaction energy of K7 and K*-mesons with the baryon field 
which is produced by the pion field in the volume V_ may be neglected in 


(*) It is worth noting that the use for the calculations of W,,, of closed expressions 
obtained in the paper (7) is very difficult because of catastrophic disappearance of terms 
when subtracting in these expressions terms of nearly equal magnitude. This cireum- 
stance is clearly expressed also when calculating the numerical coefficients, determining 
the expressions W,,, in (7). 

(*) S. GOLDHABER: Proceedings of the Fifth Annual Rochester Conference (1955); 
L. vAN Rossum and L. T. KERTH: Bull. Amer. Phys. Soc., 1, 385 (1956); W. W. CHUPP 
et al.: Suppl. Nuovo Cimento, 4, 359 (1956); S. J. LinpEBAUM, and C. L. LUKE YUAN: 
Phys. Rev., 105, 1931 (1957); G. HARRIS, J. OREAR and S. TAYLOR: Phys. Rev., 101, 
1214 (1956). 

(*) G. D. RocHEstTER and C. C. ButLER: Rep. Progr. Phys., 16, 364 (1953). 

(°) D. I. BLoHINcCEV: Usp. Fiz. Nauk, 61, 137 (1957). 
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comparison with the interaction energy with the baryon field, localized in the 
volume Vx 

The structure of the compound particle thus created corresponds to the 
hypothesis of « global interaction » of baryons with the pion field (!!) sug- 
gested recently by SCHWINGER and GELL-MANN. It is closely connected with 
the nucleon model proposed by Academician I. E. TAMM (!°). In Tables II-IV 
this case of the choice of space volumes is marked by the subscript (DÒ 

2) If we assume that strange particles in contrast to the nucleons inter- 
act with the pion field comparatively weakly and the interaction energy of these 
particles with nucleons produced in the volume V_ by the pion field is negligibly 
small in comparison with the baryon field localized in the «core» of the compound 
particle V,,, then all the strange particles will be produced in the volume V,.; 
while the pions, nucleons, isobars and the corresponding antiparticles will be 
produced in the volume V_. 

Such a choice of the compound-particle structure corresponds to the hypo- 
thesis of nucleon structure suggested by NEGANOV (14) or PEASLEE (1°) and 
does not contradict the nucleon model suggested by I. E. TAMM. 

In Tables II-IV this case of choosing the compound-particle structure is 
marked by the subscript (2). 

Because of strong interaction in the regions Vy, and Vx statistical equi- 
librium can be established. The statistical theory of multiple production sug- 
gests that the time of the collision is greater than the time of the establishing 
of such an equilibrium. 

In the collision of two high velocity nucleons the strong interaction re- 
gion V, may be considered as forming as a result of the collision of nucleon 
cores and the region V_ as formed by the collision of external regions of nucleons 
and the extension of the region V,. 

In more strict theory together with the possible introduction of « the strong 
interaction volumes » V,, and V_ in the expressions for the statistical weights 
the interaction constants must be explicitly taken into account. 

The difference in the interaction constants for different kinds of particles 
in the considered variants of theory is shaped to some extent by the intro- 
duction of different volumes. 


(i.e., in the «core» of the compound « particle »). 


(4) V. Ju. Vorrusa, M. JA. Danysz and F. E. Law: Information on VII Annual 
Rochester Conference (Doubna, May 1957). 

(12) I. E. Tamm: Report at Theoretical Physics Congress held in Seattle (USA, Sep- 
tember 1956); Zu. Éksper. Teor. Fiz., 32, 178 (1957); V. M. Koprov and V. S. BARA- 
Sunxov: Zu Éksper. Teor. Fiz. (in press); E. G. BUBELEV: Zu Bksper. Teor. Fiz. (in 
press). 

(13) M. A. MARKOV: Dokl. Akad. Nauk SSSR, 101, 64 (1953). 
(14) B. S. NEGANOV: Zu. Bksper. Teor. Fiz. (in press). 

(15) D. C. PEASLEE: Phys. Rev., 105, 1934 (1957). 
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4. — The results of the calculations. 


When making the calculations we used the values of masses of elementary 
particles, which are given in Table I in terms of the nucleon masses: 


Apre 


The results of the calculations of respective probability of possible reactions 
in proton-proton and proton-neutron collisions are given in Tables II and III. 
In these tables W, is the probability of the reaction in percentages, 


(4) W.=YW(n), 


n=0 


the total probability of the reaction: 


pp 
(5) or }> (A... B-nz) NIE N 
pn 


Since the complete analysis of all the particles produced in the particle 
collision by the charge signs is rather difficult, in Table III the division 
by the charge signs is fulfilled only for strange particles. For the calculations 
of the statistical weights in the isotopic space the tables of the paper (1) were 
used; the values of the charge factor necessary for the division by the charge 
signs were taken from the papers (‘). In the cases when the probability of 
the reaction is very small the division by the charge signs was made only for 
antihyperons since these reactions are of great interest in connection with 
the study of the possibilities of discovery and experimental study of anti- 
hyperons with the big accelerator at the Joint Institute for Nuclear Research, 
at Doubna. 

In Table IV the average values of the number of pions, nucleons, anti- 
nucleons, K and K-particles, hyperons Y and antihyperons Y produced in 
(pp) and (pn) collisions at different energies of the incident nucleon £ and 
at two above discussed asusmptions on the compound-particle structure are 


(!9) V.5. BARASENKOV and B. M. BARBAZEV: Nuovo Cimento (submitted to publ.). 
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TABLE Li: 
| (pp); 7 GeV (pp);-10 GeV (pn); 7 GeV (pn); 10 GeV 
| _ pari i Pe: e c È 
Wy % Wa % Wi % Wa % Wi % Wi % Wi % W, % 
n NNnr 
0 0.084 0.097 | 0.0190 | 0.022 0.14 0.163 0.032 0.037 
I | SS 3.88 IR 1.26 3.43 3.99 1.19 LOS 
DAB ITS. 20.5 | 8.92 10.2 | 19.3 DITA 9.74 11.2 
| 3 36.0 41.4 | 26.2 30.0 34.4 39.7 25.5 29.3 
MAD 22 25.6 \29.1 BOLD DBD 25.5 29.0 33.4 
| 5 6.05 6.98 | 14.9 bri 5.78 6.70 14.4 16.6 
6 0.49 0.56 2.82 Bee 0.47 0.54 Dl 3.12 
| 0.58 0.65 0.43 0.50 
i | 
NAKnr 
| 0 0.68 | 0.035 0.28 0.0145 it} 0.059 0.47 0.025 
1 2.34 0.120 11,3383 0.068 2.00 0.103 1.18 0.061 
2 hea 0.090 1.79 0.091 1.90 0.098 1.95 0.099 
| 3 023 0.012 0.67 0.033 0.24 0.012 0.72 0.036 
4 0.022 0.0211 0.07 0.024. 0.022 0.0211 0.07 0.024 
| | 
NXKnr 
0 128 0.064 0.52 0.027 133 0.069 0.57 0.029 
1 4.04 0.207 2.44 0.124 4.36 0.225 2.66 0.136 
2 3.22 0.165 3.63 0.183 3.16 0.162 3.64 0.185 
3 0.34 0.017 1.28 0.064 0.33 0.017 1.25 0.063 
4 0.031 0.0715 0.14 0.027 0.030 0.0215 0.14 0.027 
2NKKnr 
0 0.022 0.025 0.020 0.023 0.028 0.032 0.026 0.030 
Il 0.069 0.079 0.074 0.084 0.063 0.073 0.069 0.080 
D 0.037 0.043 0.087 0.099 0.038 0.045 0.091 0.104 
3 0.0734 0.0239 0.032 0.077 0.0233 0.0238 0.034 0.076 
2A2Knr 
| & 
| 0 0.0220 0.0524 0.0224 0.0529 0.0733 0.0541 0.0741 0.0°50 
1 0.0336 0.0943 0.0212 0.0°15 0.0338 0.0°47 0.0213 0.0°16 
AX2Knr 
0.0240 0.0°49 0.0781 0.0°48 0.0743 0.0°53 0.0288 0.0°53 
1 0.0286 0.0°10 0.0238 0.0°45 0.0893 0.0511 0.0241 0.0°50 
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Tape II (continued). 


(pp); 7 GeV (pp); 10 GeV (pn); 7 GeV (pn); 10 GeV 
E : F = | 
Wee oa We el LIZ (War Soe Wa oes omnia e | Wi % 
n 2x2Knr 
0 0.0236 | 0.0544-| 0.0273 | 0.0588 | 0.0239 | 0.0548 | 0.0279 | 0.0596 
1 0.0355 | 0.067 0.0233 | 0.0539 | 0.0355 | 0.087 0.0233 | 0.0540 
NE2Knr 
0 0.0266 | 0.0326 | 0.012 0.0349 | 0.0284 | 0.0334 | 0.016 0.0363 
1 0.0265 | 0.0320 | 0.021 0.0982 | 0.0261 | 0.0317 | 0.020 0.0378 
2 0.0341 | 0.054 0.0264 | 0.0319 | 0.0341 | 0.055 0.0265 | 0.0319 
3NNnx 
0 0.072 0.083 0.52 0.059 0.092 0.106 0.66 0.76 
1 0.0218 | 0.0221 | 1.34 1.53 0.0219 | 0.0226 | 1.23 142 | 
Q | 1.03 17 1.09 12.680 
3 | 0.05 0.05 0.05 0.05 


given. The data for the energies E = 3 and H=5 GeV have been taken from 
the Tables of the papers (°). 

In Tables II-IV the exponent of the zero means the number of zeros after 
the comma, for instance: 0.0°57 = 0.00057. 


5. — Discussion. 


As it was mentioned above the taking into account of the compound 
particle structure is quite necessary for bringing into agreement theory with 
experiment (*°). However, in spite of the fact that the probabilities of the 
strange particle production by both variants of the choice fo the compound part- 
icle structure are different by one order it is difficult to determine which of the 
considered variants is more close to the experiment. An interesting conse- 
quence from the assumption on the compound particle structure is the maxi- 
mum of the strange particle ratio to the number of pions and nucleons, in 
the energy region of 7 GeV. 

The cause of it consists in the probability of the reaction (NN) > (NAKnz) 
and (NN) +(NXKnr) decreases when passing from E=5--7 GeV to E=10 GeV, 
the probability of the reaction with greater number of K-mesons, for instance 
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‘ (pp) 10 GeV (pn) 10 GeV 
Reaction = = Se È 
Wi % We % Wi% Wi % 
2NAA 0.11 0.0215 0.19 0.0726 
2NAAT 0.07 0.0210 0.05 0.036 
2NAZ(AZ) 0.15 0.0221 0.17 0.0223 
2NAZ(AZ)T 0.06 0.038 0.07 0.039 
PAND I» 0.19 0.0226 0.21 0.0229 
ON Dr 0.05 0.0°6 0.05 0.036 
2NEE 0.0293 0.0313 0.012 0.0716 
NAEY 0.0252 0.0524 0.0256 0.0527 
NE2A 0.016 0.0°74 0.027 0.0413 
3NAK 0.0327 0.0414 0.0335 0.0418 
2NNAK 0.0381 0.0441 0.07104 0.0453 
SNZK 0.0313 0.0567 0.0474 0.0°73 
2NNXK 0.0°40 0.0420 0.0°43 0.0422 
e. pp (np) 
Reaction = = 
Wi % Wa % Wy % We % 
NAK'nt 3.05 0.156 2.65 0.137 
NAK°nr 1.99 0.102 2.65 0.137 
NDS) IK The 2.14 0.110 LAT 0.060 3 
DTT 2.14 0.110 1.86 0.096 0 
Ke 1.79 0.092 1.56 0.080 6 
PO 1.10 0.056 TESO 0.080 6 
Diga ISO 0.069 1.86 0.096 0 
SI 0.33 0.017 TESTS 0.0603 
2NK° K'nr 0.040 0.046 0.030 0.035 
K° K° 0.034 0.040 0.035 0.040 
Kan td 0.034 0.040 0.035 0.040 
1k ie 0.018 0.021 0.040 0.035 
NE K*K nr 0.0221 0.0484 0.0383 0.0433 
=o K'R° 0.0743 0.0317 0.0237 0.0315 
= K'K* 0.0221 0.0484 0.0718 0.0473 
E° K° K° 0.0710 0.0438 0.0718 0.0473 
ie Eke’ 0.0719 0.0476 0.0737 0.0915 
E K° K® 0 0 0.0°83 0.0433 
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(pp) (np) 
Reaction ze 
hg Wi % Wi % W.% 
nz-NAK* 2.48 0.126 2.20 0.112 
Ke 1.66 0.084 2.20 0.112 
nz-ND*K* 1.86 0.094 1.08 0.091 
x*K° 1.86 0.094 1.64 | 0.128 
SVI" 1.58 0.080 1.40 | 0.111 
bo Ko 1.04 0.052 1.40 0.111 
ae 1.23 0.063 1.64 0.128 
= Ke 0.44 0.022 1.08 0.091 
nn: 2NK® K* 0.065 0.080 0.047 0.059 
K° Ko 0.053 0.066 0.057 0.071 
Tele: 0.053 0.067 0.057 0.071 
KK? 0.030 0.039 0.047 0.059 
nr:NE° K*Kt 0.0267 0.0825 0.0236 0.0814 
E° KtK? 0.0134 0.0850 0.0117 0.0844 
SK Ke 0.0267 0.0325 0.0258 0.0323 
| E° K° Ko 0.0242 0.0816 0.0258 0.0323 
| E K'R° 0.0283 0.0332 0.0117 0.0844 
= Ko Ko 0.034 0.041 0.0236 0.0314 
| (pp) (np) 
Reaction 
| Wi % Wi % Wea Oe Way 
| nr 2NAZ (AZ) 0.136 0.02184. 0.073 0.0398 
| x0 (Ax? 0.062 0.0882 0.087 0.02118 
SAD) 0.016 0.0321 0.073 0.0298 
| nm:2NDX- 0.122 0.02164. 0.084 0.02113 
| x0 0.075 0.02101 0.087 0.02117 
yr 0.042 0.0356 0.085 0.02114 
| 2A 
| nm-(224+2A)E" 0.019 0.0589 0.016 0.0576 
| RO 0.0220 0.0891 0.016 0.0576 
nr:2NE5- 0.0270 0.0495 0.0260 0.0482 
RO 0.0223 0.0432 0.0260 0.0482 
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TABLE IV. 
| TT I\ K K NE Ww 
ane i zi as : aa ue . 
| (1) | (2) (1) (2) (1) (2) (1) (2) | (1) (2) (1) | (2) 
is | x È bra i | 5 
| (pp)-collision 
3 | 1.92 | 2.00 = — 0.056 | 0.0226 | 0.0592 | 0.0596 | 0.056 | 0.0226 — — 
E30 |) Bags || PAT — = 0.13 | 0.0269 | 0.041 | 0.0347 | 0.13 | 0.0764 _ = 
7 | 2.88 | Sal | 0.0°74 0.0385 | 0.14 |0.0286 | 0.0213 | 0.0215 | 0.14 0.0271 — — | 
(0) | 3.35 | 3.60 | 0.029 | 0.033 | 0.12 | 0.0290 | 0.0221 | 0.0228 | 0.13 | 0.0763 0.0287 | 0.0912 
| | | 


(pn)-collision 


Si 


|199| — — | 0.065 | 0.0231 | 0.0411 | 0.0412 | 0.065) 0.0931) — | — 


3 | 1.88 | 

5 | 2.40 | 2.63 — — 0.14 | 0.0274 | 0.0846 0.0353 | 0.14 | 0.0269 =) 
7 | 2.84 | 3.08 | 0.0394 | 0.0211 | 0.15 | 0.0290 | 0.0713 | 0.0715 | 0.15 0.0275 — — 
0 


2530356001030 (0:03 5 0st3 0.0294 0.0222 0.0229 | 0.14 | 0.0266 | 0.010 | 0.0313 
| | | | 


| |— = — — 


though it increases with energy cannot, however, considerably change the 
total probability, because K-mesons are produced in the volume V, and 
the probability of their production is strongly suppressed in comparison with 
pions and nucleons. 

The decrease of strange particles number in the high energy region (H>10GeV) 
would be in contradiction with the cosmic ray experiments (°). 

However, at high energies, in variant 1 the reactions of the type: 


(NN) >2N2Ann); (NN) > (2N2Enx), ete. 


can give considerable contribution and the number of strange particles may 
increase. More detailed consideration in the case of variant 2 is also neces- 
sary. These questions are being studied now. 

The number of the produced antinucleons depends slightly on the variant 
chosen for the compound-particle structure and is approximately equal to one 
antinucleon per two-three tousand pions for the energy E = 7 GeV and one 
antinucleon per one-two hundred nucleons at the energy E = 10 GeV. This 
result can be put in agreement with experimental data by the order of 
magnitude: one antinucleon per 6:104:5-105 of the negative pions obtained 
when bombarding copper targets with protons energy H = 6.2 GeV (17:15), 


(17) O. CHAMBERLAIN, E. SEGRE, C. Wicaup and T. YpsiLantIs: Phys. Rev., 100, 
947 (1955); H. H. HECKMAN, F. M. Smirn and W. H. BARKAS: Bull. Amer. Phys. Soc., 
1, 186 (1956). 

(18) W. H. BARKAS et al.: Phys. Rev., 105, 1037 (1957). 
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if we take into account anormally small antinucleon range in nuclear matter 
in respect to the annihilation (1%), 


lan gi ge 


Due to this fact antinucleons produced on one of the nucleons in the nucleus 
of Cu practically annihilate on the following nucleon, and from the nucleus 
only antinucleons produced on the nucleus surface may emerge. 

Taking into account that the incident nucleon interacts in the Cu nucleus 
with the « tube » containing on the average3 — 4 nucleons and that 0,,/0,,4 ~ 0.25 
we obtain the probability of the incident nucleon reaching the last nucleon in 
the tube with small energy loss W~5-10-°. 

The basic contribution to the production of antinucleons on the nucleus 
of Cu is given by these nucleons, because the cross-section of the antinucleon 
production quickly decreases with decreasing energy (see Tables II-IV). Due 
to the above considerations the ratio of the number of the produced anti- 
nucleons to the number of pions decreases up to the magnitudes exper- 
mentally observed (19:29). 


In conclusion we express our gratitude to D.I. BLOHINCEY for his interest 
and discussions. We are grateful to L. G. ZASTAVENKO for fruitful discussions 
and V. M. MAKSIMENKO and I. L. ROZENTAL’ who presented us the results 
obtained in the papers (*7) before their publication. 


(19) It is worth taking into account that the nucleon in the nucleus of Cu possesses 
a Fermi energy of 25 MeV; due to this the colliding nucleons in the experiments (17-18) 
possess the same energy in the center of mass system as for a free nucleon collision 
at an energy E= 7 GeV. 

(2) An| analogous phenomenon is well-known in the atomic reactor theory. It is 
called « block-effect ». In our case the nuclei of Cu play the part of « blocks ». 
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ContENUTO. — 1. Introduzione. — 2. La capanna. — 3. Sistemazione generale 
del laboratorio. — 4. Comunicazioni e servizi. — 5. Il locale della ca- 
mera. — 6. Attuale impiego del laboratorio. 


1. — Introduzione. 


Presso l’Istituto di Scienze Fisiche della Universita di Milano si condu- 
cono, da più di un decennio, ricerche sulla radiazione cosmica. Per lo studio 
delle interazioni nucleari a grande energia (> 10 GeV) è stata in questi ultimi 
anni progettata e costruita una grande camera di Wilson avente dimensioni 
di ingombro di 156x156x98 em*, destinata, almeno nei prossimi anni, 
a ricerche a quota di montagna (*). 

Poichè un conveniente montaggio di tutta l’apparecchiatura avrebbe richie- 
sto un locale di non meno di 50 m? di superficie, non è stato possibile usu- 
fruire dell’esistente Laboratorio della Testa Grigia (*), essendo quest’ultimo 
già in precedenza prenotato da alcuni gruppi di ricercatori di altre Università. 

Si è allora cercato, in altre località di montagna, un edificio avente requi- 
siti tali da poter venire adattato a laboratorio. La Società Edison, da noi in- 
terpellata nella persona dell’Ing. CANDIANI, ci segnalò e successivamente ci fece 
visitare le seguenti località: la centrale idroelettrica dell’Alpe Devero (No- 
vara), il cantiere di costruzione degli impianti idroelettrici sull’alto Chiese 
(Trento) ed il cantiere della diga del Sabbione (Novara). 

La quota favorevole di quest’ultima località, l’esistenza in loco di servizi 


(1) E. FIORINI, R. GIACCONI e C. SUCCI: Nuovo Cimento (in corso di stampa) (1957). 
(2) G. BERNARDINI, G. Longo e E. PANCINI: Rie. Scient., 18, 91 (1948). 


9 - Supplemento al Nuovo Cimento. 
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per il trasporto di materiali pesanti (fino a 50 quintali), la possibilità di ac- 
cesso garantita tutto l’anno, l’esistenza di un collegamento telefonico diretto 
con Milano ed altre condizioni favorevoli di cui si dirà più avanti, hanno fatto 
scegliere la diga del Sabbione come sede del nostro Laboratorio. 

Per gentile interessamento del Dott. BoBBIo, Consigliere Delegato della 
Società Edison, l’Ing. Justi, Direttore Tecnico della Società Edison-Volta, 
autorizzava l’uso e l’impiego, a tempo indefinito, di una capanna del can- 
tiere posta a quota 2518 m s.m. a circa una cinquantina di metri dalla stazione 
terminale della funivia Morasco-Sabbione. 


2. — La capanna. 


La capanna adibita a Laboratorio, che aveva ospitato dirigenti e tecnici 
preposti alla costruzione della diga, ha pianta rettangolare con lati di 15.20 m 
e 5.20 m; l'altezza è di metri 2.50. Le coordinate geografiche sono all’incirca 
46° 25' 30” di latitudine Nord e 4° 08' 52” di longitudine Ovest (Monte Mario). 
L’asse maggiore della capanna giace nei quadranti Nord-Ovest e Sud-Est, e 
forma un angolo di circa 40° con la direzione Nord-Sud. 

Posta su di un piccolo rilievo roccioso sovrastante il lago e la diga del Sab- 
bione e, prospiciente il grande ghiacciaio dell’ Hochsandhorn la capanna (Fig. 1) 
è costruita su di un basamento in muratura di pietrame con malta di cemento, 
di altezza variabile da circa 50 cm a circa 3 m, ed ha un muro longitudinale 
a sostegno del pavimento. 


ce 
oe 


Fig. 1. — La capanna Laboratorio e la stazione della funivia. © 
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Le pareti perimetrali sono formate da pannelli di 100 x 250 em?, aventi 
uno spessore di circa 10 cm, con una intercapedine di 6 cm riempita di ver- 
miculite. La parete esterna dei pannelli, collegati fra loro a mezzo di montanti 
a doppio T, è di legno di larice, quella interna di abete. 

Il pavimento, sopraelevato rispetto alla roccia sottostante, è sostenuto 
a mezzo di travi e realizzato con pannelli da 100 x 300 em’, disposti longitu- 
dinalmente. 

La copertura, a falde inclinate di 30°, sostenuta con capriate di tavole 
di legno poste ad un metro di interasse, è realizzata con tavole di 2.5 cm di 
spessore, ricoperte con cartone catramato e lamiera ondulata da 6/10 di mm. 
Per aumentare la stabilità contro i venti che, durante l’inverno, soffiano im- 
petuosi da nord-est e sud-ovest, la capanna è puntellata e controventata con 
otto tiranti in corda di acciaio. 


3. — Sistemazione generale del laboratorio. 


Per permettere alla capanna di ospitare la grande camera di Wilson, si 
son resi necessari alcuni lavori di adattamento e variazioni alla disposizione 
preesistente dei locali. 

Data la imperfetta sutura dei pannelli, si è provveduto a migliorare l’isola- 
mento termico della capanna, per ridurre gli effetti delle variazioni esterne 
di temperatura, rivestendo le pareti interne con carta e con un foglio di poli- 
tene da 0.1 mm; sopra il politene è stato poi disposto un rivestimento di fae- 
site. Il pavimento è stato ricoperto con faesite e su questa è stato applicato 
del linoleum. 

Per sostenere alcune parti pesanti dell’apparecchiatura ed isolare la capanna 
dalle vibrazioni del grosso compressore della camera, si son dovuti costruire 
sulla roccia alcuni basamenti di calcestruzzo. Si son pure ancorate diretta- 
mente alla roccia, a mezzo di robuste putrelle di ferro, le due macchine foto- 
grafiche stereoscopiche della camera di Wilson, per impedire che le vibrazioni 
dovute alle espansioni della camera danneggiassero la bontà dei fotogrammi. 


La disposizione dei locali e la loro funzione sono mostrate in Fig. 2. 


0 100 200cm 


ne 


Fig. 2. — Pianta del laboratorio: A) Laboratorio; B) Locale per le camere; 0) Soggiorno; 


D) Camera da letto; #) Camera da letto; FY) Bagno; @) Corridoio d’ingresso. 
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Il locale A è il Laboratorio propriamente detto e contiene, oltre all’ap- 
parecchiatura della camera, i quadri di distribuzione dell’energia elettrica 
e l'impianto di termostatazione. Parte dello stesso locale è adibita ad officina 
meccanica. 

Nel locale B, destinato a contenere le camere, le lampade per Villumina- 
zione ed i condensatori sono disposti quattro travi di larice, di sezione qua- 
drata con 25 cm di lato, che sostengono, ad un’altezza variabile, un piano 
di putrelle di ferro, di 170 x 180 cm? di superficie, capace di portare fino 
a 100 quintali di materiale. 

In corrispondenza di tale piano di sostegno, per avere maggiore disponi- 
bilità di spazio, il soffitto è stato sopraelevato di circa 50 cm. 

Il locale C è adibito a soggiorno (Fig. 3) ed i locali D ed £ sono adibiti 
a camere da letto: la stanza D contiene un letto, la # due cuccette. 


Fig. 3. — Scorcio del locale adibito a soggiorno. 


F è il bagno, dotato di lavabo con acqua corrente calda e fredda, doccia 
eo WiG; 


4. — Comunicazioni e servizi. 


Al laboratorio si può accedere per mezzo della funivia Morasco-Sabbione: 
il percorso viene compiuto in circa 24 minuti primi; la base della funivia è col- 
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legata al cantiere di Valdo, e di qui a Domodossola (Novara) da una strada 
carrozzabile. 

D'inverno, quando il tratto Valdo-Morasco si chiude per la neve, la sta- 
zione inferiore della funivia può venire raggiunta attraverso le gallerie in cui 
sono sistemate le condotte forzate: il tempo che normalmente si impiega a 
superare il tragitto è di circa un’ora e mezzo. 

I materiali pesanti vengono trasportati al Laboratorio mediante una tele- 
ferica che, partendo dalla località Morasco, raggiunge il Sabbione in pros- 
simità della stazione della funivia. Durante la stagione invernale materiali 
anche pesanti, ma non troppo ingombranti, possono venire trasportati alla te- 
leferica attraverso le gallerie. 

Il laboratorio è dotato di energia elettrica che viene fornita da una ca- 
bina, distante dalla capanna una cinquantina di metri, per mezzo di un tra- 
sformatore trifase da 7500--8 000 volt primari e 230--240 volt secondari. La po- 
tenza attualmente disponibile, di circa 250 KVA, è aumentabile a 325 KVA 
impiegando un secondo trasformatore di riserva, collegabile in parallelo. 

Nel Laboratorio un trasformatore universale trifase da 20 kVA ed un tra- 
sformatore universale monofase da 15 kVA permettono di alimentare due 
altre reti interne di distribuzione a 160 volt. Una quarta rete interna alimen- 
tata da un raddrizzatore al selenio da 600 watt distribuisce corrente conti- 
nua a 24 volt. 

Una tubazione di 2 pollici di diametro, in gran parte sotterranea, fornisce 
acqua al Laboratorio prelevandola da una sorgente alimentata dal ghiacciaio 
del Blindenhorn, posta a circa 1 km dalla capanna. 

Il riscaldamento della capanna è ottenuto elettricamente a mezzo di 12 
stufe da 500 watt, alimentabili a coppie in serie o in parallelo e pilotate da 
termointerruttori bimetallici in modo da mantenere nell’interno una tem- 
peratura intorno ai 20°C. Per fronteggiare interruzioni prolungate nella di- 
stribuzione di energia elettrica la capanna è pure dotata di riscaldamento ed 
illuminazione a gas liquido. 


5. — Il locale della camera. 


Perchè le caratteristiche del locale che ospita la camera siano tali da ga- 
rantire un buon funzionamento è necessario stabilizzarne la temperatura entro 
qualche decimo di grado centigrado, anche quando le condizioni esterne sono 
particolarmente sfavorevoli. 

Per ottenere ciò, oltre a migliorare l’isolamento termico del pavimento 
e della parte esposta al vento, si sono impiegati tre termoconvettori aventi 
ciascuno una superficie utile di 20 m? disposti come è indicato in Fig. 4 ed ali- 
mentati con acqua circolante in ciclo chiuso a 20 SCE 

L'acqua proviene da un serbatoio di 450 dm? di capacità, in cui un ter- 
mometro tipo Vertex, sensibile al decimo di grado pilota due riscaldatori blin- 
dati aventi una potenza regolabile da 500 a 1000 watt. 

Una pompa della portata di 6 litri al minuto primo rimescola l’acqua nel 
serbatoio. Una seconda pompa della portata di 20 litri al minuto primo, pi- 
lotata da un secondo Vertex di eguale sensibilità, posto nel locale della camera, 
fa circolare l’acqua nei termoconvettori. Con il dispositivo impiegato la tempe- 
ratura è risultata stabile entro 3/10 di grado centigrado. 
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Per ostacolare la formazione di gradienti di temperatura nel locale della 
camera si sono disposti tre altri ventilatori. Ad un metro e mezzo di altezza 
dal pavimento ed in vicinanza della camera il gradiente orizzontale di tempe- 
ratura è risultato inferiore a 0.2 °C/m ed il gradiente verticale inferiore a 
0.5 °C/m. 


6. IONE 200m 


| 
Fig. 4. — Disposizione degli apparecchi nel la' oratorio: 1) Attrezzatura d officina; 


2) Centrale elettrica; 3) Compressore; 4) Serbatoi per aria; 5) Termostato per acqua; 

6) Apparecchi elettronici; 7) Ponte di sollevamento; 8) Rotaie; 9) Camera di Wilson 

grande; 10) Camera di Wilson piccola; 11) Lampade a flash; 12) Dispositivo di sicu- 

rezza pei flash; 13) Alimentatore elettronico; 14) Batteria di condensatori; 15) Lam- 

pada a incandescenza; 16) Macchine fotografiche; 17) Termo-convettori; 18) Sup- 
porti per piano di sostegno. 


Per ridurre eventuali riflessioni di luce, dannose per la fotografia, le pareti 
del locale contenente la camera sono state ricoperte con un drappo nero. 

Per facilitare le operazioni di montaggio e smontaggio, la camera è stata 
posta su di un carrello della portata di 10 quintali, scorrevole su rotaie fissate 
su una serie di otto pilastrini di calcestruzzo. Con semplice trazione la ca- 
mera viene trasportata, attraverso un’apposita porta quadrata, nel locale 
adibito a laboratorio sotto un ponte di sollevamento costituito da un paranco 
montato su travi di legno. 


6. — Attuale impiego del laboratorio. 


Attualmente al laboratorio si avvicendano alcuni ricercatori e tecnici delle 
Università di Milano e di Pisa (*), costituiti in un gruppo per una ricerca il 


(3) M. ConveRSI, E. FIORINI, C. GirorI, C. Russia, C. Succt e R. TORELLI: Progetto 


di una esperienza per la ricerca di particelle di massa 550 m,; Conferenza internazionale 
sulle particelle elementari, Padova, 1957 
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cui scopo è di accertare l’esistenza nella radiazione cosmica delle particelle 

di massa 550 m,, segnalate recentemente (*) da ALIHANIAN et. al. 
L’apparecchiatura è costituita da un telescopio di due grandi camere di 

Wilson (Fig. 5), associate in un selezionatore elettronico di particelle, ed ha 


(*) Comunicazione privata del prof. M. CONVERSI. 
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lo scopo di permettere la determinazione della massa, del decadimento e della 
vita di queste particelle. 

Le Fig. 6 e 7 mostrano alcuni aspetti del laboratorio durante la fase di 
montaggio degli apparecchi. 
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La costruzione del Laboratorio è stata possibile grazie al generoso e con- 
tinuo appoggio ricevuto dal Gruppo degli Amici dell’Istituto di Scienze Fi- 
siche, in particolare dalle Società Edison ed Edison-Volta, A.N.I.D.E.L., Pi- 
relli, Linoleum, Agipgas alle quali vogliamo esprimere la nostra più viva ri- 
conoscenza. 

È nostro gradito dovere ringraziare particolarmente il Dott. C. BOBBIO, 
Consigliere Delegato della Società Edison, per l'interessamento e la solleci- 
tudine con cui ha seguito il nostro lavoro; il Dott. Ing. L. JUSTI, Direttore 
Tecnico della Società Edison-Volta, per le sue numerose concessioni; ed il 
Dott. Ing. M. CALCIATI, Direttore di Zona della Società Edison-Volta, per il 
suo cordiale e generoso aiuto. 

Al Prof. G. PoLvanI, Direttore dell’Istituto di Scienze Fisiche dell’ Uni- 
versità di Milano, ed al Prof. P. CALDIROLA, Direttore della Sezione di Milano 
dell’I.N.F.N., esprimiamo la nostra più viva riconoscenza per il continuo ap- 
poggio ed incoraggiamento datoci durante l’allestimento del Laboratorio. 


SUPPLEMENTO AL VOLUME VII, SERIE X Nee le L958: 
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A Microwave Electron Spin Resonance Spectrometer (*). 


M. B. PALMA-VITTORELLI and M. U. PALMA 


Istituto di Fisica dell Universita - Palermo 
(ricevuto il 19 Agosto 1957) 
CONTENTS. — 1. Introduction. — 2. The magnet. — 3. The magnetic field. 


1. Design considerations. 2. The power stage. 3. The error-signal 
channel. 4. Composite loop behaviour. 5. Modulation and measu- 


rement of the magnetic field. — 4. The microwave system. 1. Design 
considerations. 2. The apparatus. — 5. Detection and recording of the 


magnetic absorption signal. 


1. — Introduction. 


The importance of the electron spin resonance technique has been growing 
as higher sensitivities and precisions have been reached. The equipment here 
described reaches a sensitivity which is considerably higher than the one of 
similar commonly described equipments. The high stabilization of both the 
magnetic field and the microwave frequency yields a precision which is es- 
sentially limited by the width of the absorption lines (in the case which is of 
interest for us, i.e. solid state). The result is achieved by the use of relatively 
simple circuitry. The use of a rotating magnet enables ready and accurate 
measurements of anysotropy and saturation to be performed. 

As at is well known, one is interested not only in the position of the ab- 
sorption maxima, but also in the whole plots of the real and imaginary com- 
ponents (y’ and y") of the complex magnetic susceptibility. Far from satu- 
ration, the complex susceptibility is an analytical function of the frequency, 


(*) This work forms part of a program of research undertaken with the financial 
support of the Consiglio Nazionale delle Ricerche, through the Comitato Nazionale per 
la Fisica. for which the A.A. wish to express their sincere gratitude. The work has been 
performed at the suggestion of Prof. M. SANTANGELO, whose large and valuable advice, 
criticism and help are gratefully aknowledged by the AA. 
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and y' and y" are related by the Kramers-Krònig relations, which may be 
found by means of a complex contour integration (*). 

The two plots (vs. the quasi-static magnetic field and vs. the microwave 
frequency) are not equivalent. In fact the simple relation which gives the 
maximum of the absorption curve, i.e. hr =gBH (where h= Planck’s con- 
stant, 6 = Bohr magneton, g = const = spectroscopic splitting factor) does 
not always hold, in the sense that, if for instance one ope rates at constant 
frequency, the ratio v/Hna May depend on the frequency at which one ope- 
rates (for instance, when a pre-splitting is present). Moreover, the line shape 
and width may be noticeably affected by the spin-lattice relaxation time, which 
is in turn H-dependent. It is a usual practice, in order to avoid remarkable 
experimental difficulties to get plots of 7’ and y' vs. H, although plots vs. v 
would give better information on the different quantities involved (relaxation 
times, exchange interaction etc.) and on their magnetic field- and tempe- 
‘ature-dependence. In our equipment too, the measurements are performed 
at constant frequency, but more than one frequency is available in order to 
compare experimental data; further, the available power in the K-band (about 
600 mW) renders frequency multiplication feasible. 

The electron spin resonance spectrometer which will be described in the 
following has the following characteristics: magnetic field homogeneity better 
than 3-10-4% in the sample volume (solid state samples) and of the order of 
10-2% in the volume of the whole cavity; magnetic field stability of the order 
of 2-10-3% short term, and of the order of 5-107? °/, long term (and a suitable 
precision in the measurement, using proton resonance); magnetic field sweep 
linearity, of the order of 0.1%; frequency stability of the order of LOmeg 
ultimate sensitivity, at room temperature, with a signal-to-noise ratio of about 
4:1, (using optimum modulation amplitude) about}101%-AH magnetic centers 
(where? AH’ is the line-width* in Gauss). 


Fig. 2. — Photograph of the cryostats for low temperature measurements. 


(*) The presence of poles in the area encircled by the integration contour invali- 
dates the relationships and corresponds to saturation in magnetic absorption (see, for 
instance, ref. (!)). Much in the same way, in circuit theory, the Bode relationship is 
found, which is of fundamental importance in feedback loop design and which relates 
cain characteristics and phase-shifts in four-terminal networks. In this case, the pre- 
sence of poles corresponds to the presence of non-minimun-phase-networks. 
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The above performances refer to a typical case of 3 cm wavelength and 
room temperature. The sensitivity may be further increased by an estimated 
factor of several hundreds, using the available system for low temperatures 
(down to liquid helium) and the A-band set up. 

The equipment is partially shown in Fig. 1; the low temperature system 
is shown in Fig. 2. 


2. The magnet. 


The electromagnet was designed to fulfil at the same time the requirements 
of electron spin resonance experiments and of a good general purpose laboratory 
magnet (*). 

It has a horizontal double yoke, rotating around a vertical axis (i.e. around 
an axis orthogonal to the direction of the magnetic field in the gap) and the 
rotation angle is easily measurable. The poles are of a diameter of 25 cm. 
The gap is continuously variable within large limits, and a mechanical system 
allows the symmetry around the center to be accurately preserved. There 
are a number of cylindrical precision machined pole pieces, with ground faces, 
some having Rose shims (2) suitable for air gaps up to about 11 em and others 
having independent shims (*). In every position corresponding to useful air 
gaps, the poles are in good mechanical balance so that the coincidence of their 
axes within strict tolerances is secured. Anyway, an alignement is required, for 
each setting, using three spacers and possibly tiny shims: the control is made 
mechanically first (using a micrometer) and then using the proton probe. In 
order to decrease the residual inhomogeneity, a saturation of the magnet in 
the reverse direction, followed by a slow cycling of the field (*) has been found 
useful. In this way the values specified in Sect. 1 have been obtained. 

In the coil system design, two requirements have been satisfied: high effi- 
ciency defined as the ratio (induction field in the air gap x gap)/ampere-turns; 
and low input power, in order to allow the use of the Laboratory storage bat- 
tery, as a power supply, and to simplify the problems concerning control, 
stabilization and cooling. The first feature has been obtained using the double 
yoke and a very good iron quality; this renders the required number of 
ampere-turns reasonably small. In addition, the required power has been 
kept low by the use of a large copper mass, in accord with the fact that this 
power is proportional to r?(N7)?/M (where M is the copper mass of the coils, 


(*) The magnet has been constructed by the Ansaldo San Giorgio S.p.A. using 
Terni iron of the type T11S: we are here very pleased to thank the Terni S.p.A., for the 
generous gift of the iron to this laboratory; the Ansaldo San Giorgio S.p.A. for accurate 
construction, and Dr. Ing. A. Lorri, of the Ansaldo San Giorgio, for very valuable 
advice on practical details and the design of the cooling system. 

(1) A M. Portis: Phys. Rev., 91, 1071 (1953). 

(*) M. E.. RosE:.Phys.-Rev., 58, 715 (1938). 

(*) J. D. BJORKEN and F. Birrer: Rev. Sci. Inst., 27, 1005 (1956). 
(4) K. D. Bowers, R. A. Kamper and R. B. D. Knicur: Journ. Sci. Inst., 34, 
49 (1957). 
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r the average radius of the turns and Ni the ampere-turns). Accordingly, the 
coil system consists of six independent sections, each of 100 turns, 5 ohms, 
three on each pole. The sections are separated from the yoke and each other 
by hollow pancakes in which the cooling water circulates. Moreover, suitable 
discharges are provided, in order to protect the magnet against harmful vol- 
tage surges. 

Roughly, up to 10 kilogauss, using a gap of about 10 em, the required 
current is about 1 ampere per kilogauss, when the six sections are series- 
connected. 

_ The possibility to rotate the magnet, together with the use of the cylin- 
drical resonant cavities which will be described in Sect. 42, allows precise 
and ready measurements of anisotropy. The sample is located at the center 
of the cavity, on a holder which may rotate around a diameter of the circular 
section. The radio-frequency magnetic field H,, on the sample is directed 
along the cylindrical axis of the cavity. The geometry of the system is shown 
in Fig. 3a. The directions of H, of the rotation axis of the sample, and of the 
rotation axis of the magnet, are rendered orthogonal to each other by the use 
of a Gauss self-collimating system which will be described in the following. 


Rotation axis of the magnetic field 


A Z=Cilindrical cavity axis = Hf A z= Cilindrical cavity axis=Hpf 
| =Rotation axis of the magnetic field | = 

| 

| 


xi 
a 
o) 
= o aka —— pe 
Y= Holder axis % Vi V= Holder axis 
ws 
Vi 
A 
vA 
7 2 X=H 
Fig. 3. — Illustration of the method for anysotropy measurements (see text). 


Measurements in the a’y'-plane of the crystal may be made simply by 
rotating the magnet around the ¢-axis (see Fig. 3a). Measurements in the 
z'v'-plane may be made by rotating only the crystal around the uh axis and 
keeping H along the « direction. Finally, measurements n the y'2 ~plane may 
be made by rotating first the sample holder by 7/2 around the y-axis (thus 
obtaining the orientation shown in Fig. 3b), and then by rotating the magnet. 
The three sets of measurements are sufficient for the determination of the 
anisotropy in magnetic resonance absorption, often independently of the initial 
knowledge of the directions of the principal magnetic axes (ae 

In order to orientate the cavity together with the sample holder, as shown 
in Fig. 3a, a Gauss self-collimating telescope is used, which fits into a small 
hole. This hole has been precision-machined in the yoke, in a direction which 


(5) M. B. PALMA-VITTORELLI, M. U. Patma, D. PALUMBO and M. SANTANGELO: 
Nuovo Cimento, 2, 811 (1955). 
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passes through the center of the gap and is orthogonal to both the pole and 
the rotation axes. The sample holder has a head, having a flat face ortho- 
gonal to the rotation axis. On this face a small mirror is secured, which is 
used to obtain the self-collimation condition, thus rendering the holder axis 
coincident with the y-axis. A diametral line is engraved on the mirror: 
if a suitable microscope is used in the hole instead of the self-collimation 
telescope, the rotation of the holder may be measured (*). 

This method eliminates the lack of accuracy in anisotropy measurements 
due to the setting and re-setting of the sample in three known orthogonal 
directions. Furthermore, it allows in any case the measurements to be made, 
in one single run, with a consequent saving of time and refrigerants, (if used). 

One extra feature of the rotating magnet is the possibility of ready and 
precise measurements of saturation parameters, operating at a constant micro- 
wave level (*). Indeed, using a cavity in which H,; (on the sample) lies in the 
xy-plane of Fig. 3a, one may measure the departure from the law which gives 
the absorbed power (proportional to H?, sin? H,H.;), as the magnet is being 
tated. The saturation parameter may be then easily computed. 


3. — The magnetic field. 


3:1. Design considerations. — In order to stabilize and control the magnetic 
field of an electromagnet which is energized by a motor-generator set, the 
motor-generator itself may be used as the power stage of the control loop, by 
feeding the amplified error-signal into the generator field-control winding. 
Such a system, which was originally developed by SoMMERS, WEISS and 
HALPERN (7), has been widely applied (*°4), in spite of considerable difficulties 
arising from the non-linearity of the generator and large phase-shifts. 

If a fixed voltage power supply or a large storage battery is available, as 
in our case, a very simple stabilization system may be obtained by the use 
of a variable impedance in series with the magnet. Such a variable impedance, 
controlled by a servo-loop, may correct the possible fluctuations or shifts of 
the magnetic field. For instance, one may use as a variable impedance element 
a number of tubes in parallel, the grids of which are driven by the amplified 
error-signal. This is the principle of operation largely used in many low-power 
stabilized power supplies of common use. 

If the dominant lag of the loop, up to the frequency corresponding to a unitary 
open-loop gain, is given by the magnet, the loop itself is inherently stable, 

(*) The whole optical system has been constructed by the Officine Galileo S.p.a.; 
the help of Prof. G. Giorrr is gratefully ocknowledged. 

(6) G. FEHER: Bell. Syst. Tech. Journ., 26, 449 (1957). 

(*) H. S. Sommers, P. R. Weiss and W. HaLPERN: Rev. Sci. Inst., 20, 244 (1949); 
22, 612 (1951). 

(8) J. R. WoLFF and M. S. FREEDMAN: Rev. Sci. Inst., 22, 736 (1951). 

(9) M. W. P. STRANDBERG, M. Tinkuam, I. H. Sort and C. F. Davis: Rev. Sci. 
Inst., 27, 596 (1956). 

(9) F. T. Hepecock and F. Hunt: Rev. Sci. Inst., 27, 970 (1956). 


(1) H. W. BopE: Network Analysis and Feedback Amplifiers Design (New York, 1945). 


(1°) J. G. THomASsoN: Linear Feedback Analysis (London, 1955). 
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since the slope of its controlled high-frequency cut-offe haracteristic has atmost 
a value of — 6 db/octave (1'). In such a case the phase-margin would be 90°. 
It follows that we may use at most in the error signal channel, either a single 
lag having a time constant comparable with the one of the magnet, or a number 
of lags of smaller time constants up to the reduction of the phase-margin to, 
Sava TUE (BED 
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Fig. 4. — Block-diagram of the composite loop for stabilization and control of the ma- 


onetic field, and detail of the ‘athode-follower power stage (for tubes and parasitic 
suppressors, see text). 


In order to further narrow the bandwidth of the error-signal channel 
(i.e. to narrow the bandwidth and to obtain a related reduction of drifts and 
noise of the high-gain amplifier) and in order to get the extra advantages 
which will be specified later, we have used a composite loop. The block- 
diagram of the composite loop is shown in Fig. 4a. The part of the circuit 
within the dotted line (power stage) is shown in detail in Fig. 4b, where a 
single tube stays for a bank of up to 96 parallel EL34 tubes. The gain of this 
stage, supposing that the magnet is mounted as a plate load, is indicated by A. 
The cathode-follower action is symbolized by the introduction of the feedback 
amount f. The load consists of the magnet. The error-signal channel is 
actuated by: the feedback /°, the error-signal amplifier, which has a gain A”, 
and a matching cathode-follower. 

Let us think of the composite loop in terms of a single loop, the feedback 
amount of which we. will label by Bum. Let Au be its gain without feedback, 
and let 
Aes 


5 5 
i Prg A Sia 


y 


(1) G= 


be the equivalent closed-loop overall gain. By equivalent open-loop gain we 
will mean the product Acme: 
One may easily see, by inspection of Fig. 4a, that 


TAR 
- C= —~ 4 
(2) 1—AA'(p'+ (B/A')) 
so that 
(3) Aan Bean = AA'(B'+ BIA‘) © 
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The above equation will be written in its explicit form in Sect. 3°4. There 
we will see that, for a given high-frequency cut-off characteristic of the equi- 
valent open-loop gain, we may let A’ drop faster, with the result that the 
bandwith of the error-signal amplifier may be narrowed or the gain may be 
increased. 

Let us now look for the behaviour of the different stages. 


32. The power stage. — The power stage, as shown in Fig. 4b, contains 
up to 96 EL34 tubes in parallel, having the magnet as a load, in a cathode- 
follower mounting. The power supply is the laboratory storage battery which 
gives up to 250 V. The heating current is supplied from a stable d.c. source. 
A very well stabilized, variable voltage supply is used in order to bring the 
tubes in the proper working range, and to give the possibility of a manual 
control of the field. Furthermore, it allows an easy switching of the magnet, 
without harmful voltage surges. 

The cathode-follower mounting offers the following extra advantages: since 
its voltage gain is about unity, larger noises and drifts are allowed, at constant 
overall gain, at the input (i.e. at the output of the high gain error-signal 
amplifier). Further, it makes the phase-shifts in the power stage negligible. 
It is to be noted anyway that the cathode-follower action reduces the 
voltage-to-voltage phase-shifts, while the magnetic field-to-voltage phase-shifts 
are determined only by the time constant of the magnet. 

An elementary analysis shows that in the entire frequency range in which 
we are interested (parasitic suppressors disregarded), it is: 


(4) 4  ER4+0) + OL + jLol 
Ja = ph (R+0) + wl? 


3 


where 4 is the amplification factor and o the plate resistance of the EL34’s 
which are regarded as a single tube, and , L are respectively the resistance 
and self-inductance of the magnet. Eq. (4) shows that the gain A is a function 
of the number of tubes used, of the resistance of the magnet (depending on 
the connection of the different sections), and of the air gap in the magnet. 
For typical values, related to 50 triode-mounted tubes, series-connection of 
the magnet sections, and a gap of less than 10 cm, o = 12 Q, uw = 6, R = 30 Q, 
L=100H; A is almost frequency-independent and may be taken as equal to 
about 5, up to © 1000. At higher values of © the parasitic suppressor’s 
time constant becomes of importance and determines the high frequency cut-off 
characteristic, which therefore has a slope of — 6 db/octave at most. 
Using the mounting exemplified in Fig. 4b, we have: 


(5) 1 


but lower absolute values are obtained if some of the sections are mounted 
as plate load. 


3°3. The error-signal channel. — Two different types of error-signal channels 
have been developed, using respectively a magnetic field and a current diseri- 
minator. 
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In the first, the principle developed by Kip and ARNOLD (12%) and 
then applied also by STRANDBERG and co-workers (°) is used: a signal 
is given by a coil which rotates in the magnetic field at a speed of 
25 r. p.s. This signal is compared in a Poggendorf bridge (potentiometer), 
with the one given !by a reference generator driven by the same axis; 
the latter signal may,be conveniently adjusted in amplitude and phase, for 
proper operation of the difference-voltage detector. The latter is a cathode- 
follower, in order to preserve the linearity of the system, and gives the dif- 
ference between the flip-coil voltage and the fraction of the reference gene- 
rator voltage corresponding to a given setting of the potentiometer. The latter 
is a helical potentiometer, which may be set manually or swept by linking 
its shaft to the drive-shaft extension of the Speedomax recorder (see Sect. 5). 
A simple switching system allows an easy change of the sweep range and 
amplitude, by putting suitable resistances in series and in parallel with the 
potentiometer. For a given potentiometer setting, a field fluctuation or shift 
modulates the 25 Hz carrier generated in the rotating coil and gives rise to 
an error-signal which corrects the field. Inversely, a given change of the po- 
tentiometer setting (i.e. of the used fraction of the reference voltage), gives 
rise, through the loop, to the wanted shift in the field value. The error-signal 
is amplified in a low frequency feedback amplifier and phase-detected in a Brown 
Converter. The RC filter which follows the phase-detector may have a time 
constant up to the order of the magnet’s and its bandwith is in any case small 
compared with the bandwidth of the amplifier. So doing, the RC filter is the 
dominant lag in the error-signal channel. (Should the side-bands undergo 
appreciable phase-shifts in the a.c. amplifier, they would be opposite in sign, 
and the phase-detector action could not cancel them, but would give in this 
case the average of their absolute values (!*)). This system offers the advan- 
tage of giving a good linearity when sweeping the field from very low values 
up to a number of kilogauss. For higher field values, the linearity is no longer 
preserved (°). Moreover, in order to have negligible phase-shifts in the 
amplifier, one has either to have a large bandwidth, or to limit the gain. This, 
together with the fact that the brushes of the flipping coil are noisy, gives in 
any range a limit of the order of half a gauss in the field stabilization even 
with the use of the composite loop. The flipping coil provides in any case a 
system for an approximate measurement of the field (7*%’*). 

When a higher stability is wanted, a d.c. error-signal channel is used- 
A high-rating 10 resistor in series with the magnet gives a voltage propor. 
tional to the current, and having the same phase of it. The same helical potentio- 
meter bridge as above is now used to compare this voltage with a fraction of 
a reference battery voltage. The resulting error-signal is then amplified in 
a d.c. amplifier. The high-frequency cut-off characteristic of the channel is 
again determined by some form of RC filter like the one following the phase- 
detector in the first case. The d.c. amplifier bandwidth again has to be larger 
than the filter’s bandwidth, but now, compared with the flipping coil case, a 
given high-frequency cut-off slope, in db/octave, and a given amplifier gain 
correspond to a much smaller bandwidth. 


(1248) R. D. ArnoLD and A. F. Kip: Phys. Rev., 75, 1199 (1949). 
(3) M. B. PaLma-VirtoRELLI, M. U. Patma and D. PaLumBo: Nuovo Cimento, 
6, 1211 (1957). 
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In both cases, the resulting error-signal is fed into the power stage through 
a matching cathode-follower stage, which uses two push-pull tubes in order 
to have a larger output voltage swing. 

The expressions of f’ are the same in both cases, and the numerical values 
too, as long as magnetic saturation is negligible: 


=~ 1 — NONE 
(6) p' “è py 9 3 ) 
L == 2? Ty 

(where 85) ~ — 1/30 is the £' value corresponding to w=0 and ty, is the 
magnet time constant which, for the typical values given above, is about 3 s). 
Eq. (6) may be easily obtained, since fp’ is equal to the ratio between the non- 
amplified error-signal voltage and the magnet voltage variation to be cor- 
rected. 

The A’ expressions and values are likewise the same for both cases, namely, 
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where A, is the d.c. (or carrier-frequency) gain, and t, the time constant of 
the dominant lag. 


S4. Composite loop behaviour. — Using eqs. (4) to (7), one gets 


AA‘Bo 


(8) Peg 3 : 9 TRA. 
(1 + w?t5,)(1 + 0275) 


(1 + w?ti,)(1 +23). . 
di sù JO(Tr+ Ty) 


1 = @2TxTr — aT 4 
Pod 0 


The above equation, which yields both gain and phase-shift of the equivalent 
open loop as a function of the frequency, proves that the equivalent loop is 
«absolutely » stable, since the phase-shift cannot reach 180°. Therefore, the 
loop gain may be raised as long as other time constants do not invalidate 
eq. (8) (which does not happen in our case up to loop gains of about 300 using 
the d.c. error-signal channel). Finally, it is worthwhile to note that the asymp- 
totic value of the gain as given by eq. (8), is A: the real high-frequency cut-off 
characteristic is determined by the parasitic suppressors (see Sect. 3°2). 

The foregoing analysis exemplifies a typical situation which occurs when 
one uses the spectrometer in X-band. Anyway, the working conditions of 
the loop may be largely varied (and the above analysis may be easily extended 
to the different conditions), by varying, for instance, the number of the EL34 
tubes; or the value of f; or feeding independently some of the sections of the 
magnet; or by drawing part of the current which energizes the magnet through 
a high wattage resistor which parallels the EL34’s. The latter case would be 
a rather unusual one with our equipment, because of the large number of tubes 
available in the power stage. Anyway, since the present work has been com- 
pleted, an apparatus using such a resistor across a small number of parallel- 
tubes, has been described (!*). 


(14) R. J. ABRAHAM, D. W. OvenaLL and D. H. Wuirren: Journ. Sci. Inst.. 34. 
269 (1957). 
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3:5. Modulation and measurement of the magnetic field. — It is very con- 
venient to modulate the magnetic absorption signal via a modulation of the 
magnetic field which is small compared with the linewidth. Also, it may be 
convenient sometimes to use a field modulation as large as the width of the 
whole line or spectrum, in order to get a convenient oscillogram. 

In our case the modulation may go up to 30 gauss peak-to-peak and is ob- 
tained using a pair of flat coils arranged in the gap in such a way that the 
modulation occurs almost only within the gap rather than in the massive iron 
core. The coils are transformer-coupled to a power amplifier which uses a 
bank of eight EL34 tubes in parallel, and is driven at 59 Hz by a very stable 
oscillator. The same oscillator provides the ‘reference signal used for the phase- 
sensitive detection of the absorption signal (see Sect. 5). A small auxiliary 
coil and a V.T.V.M. provide monitoring and measurement of the magnetic field 
modulation. 

An approximate measurement system of the magnetic field is given by 
the rotating coil and a Poggendorf bridge, similar to the one used in the a.c. 
error-signal channel (see Sect. 3°3). 

In addition, a precision measurement is given by a proton-resonance 
signal provided by a small probe which is part of a circuit similar to the one 
described by KNOEBEL and HAHN (1°). In our case if may be used either 
with the original frequency-modulation system, or by taking advantage of the 
magnetic field modulation. The proton-resonance frequency is measured by 
comparison with a secondary frequency standard. Since the proton-resonance 
oscillator is very stable, the system provides a very simple and accurate mo- 
nitoring of the magnetic field stability and homogeneity, simply by the use 
of a very small field modulation. The effect due to « wiggles » in the resonance 
signal is, for our work, negligible. 


4. — The microwave system. 


41. Design considerations. — As we have seen in Sect. 1, the absorbed 
power (far from saturation) is proportional to Hî;. It follows that the sensi- 
tivity in the detection of magnetic absorption is increased putting the sample 
in a high-Q cavity. The absorption may be detected measuring the variations 
of the r.f. power (or amplitude) reflected or t ‘ansmitted by the cavity. The 
use of a transmission or a reflection cavity yields almost identical sensitivities 
as it is shown, for instance, in the interesting paper by G. FEHER (°) on the 
ultimate sensitivity. We have chosen the use of a reflection cavity, since for 
low temperature work it offers the advantages of less space requirements, 
less thermal dispersion and less gas-bubble noise. 

When a reflection cavity is used, the magic-T technique is convenient. Let 
us connect for instance the cavity and a variable load respectively to the 
side-arms 1 and 2 of the magic-T. Let the r.f. power come from arm 4 
(H-arm) and a matched detector be mounted on arm 3 (H-arm). When mag- 
netic resonance occurs, a variation of 7" gives rise to a variation of the reflected 
amplitude; a variation of y' gives rise to a variation of the resonance frequency 


(15) H. W. KNOEBEL and E. L. HAHN: Rev. Sci. Inst., 22, 904 (1951). 
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of the cavity i.e. to a variation of the phase of the reflected wave, together 
with a second-order variation of its amplitude (if the resonant-frequency shift 
is small). If H, is the reflected amplitude from the cavity, if #, is the reflected 
amplitude from the variable load and if g is the relative phase at the magic-T, 
the detector is reached by a power proportional to 


Ren Oe end Bi 
(9) Hi = Se Ss — E,E, cos , 


whereas the reflected power toward the r.f. source is proportional to 


te wD eee 
(10) Ee=>+ 5+ HE, cose. 


If the effect of the magnetic resonance is to change H, by an amount dH! 
and g by an amount dg (where dH, and dp are small, as we will assume if not 
otherwise specified), the power on the detector is now (to the first order) 
proportional to: 


(11) E,; = E + (#,— E, cosg)dE, + E,E, sing dp 
and the reflected power toward the r.f. source is proportional to 
(12) EY = HE? + (#,+ E, cosq) dH, — HE, sing do . 


Eqs. (11) and (12) show that the total change in the reflected power is 
proportional to 2H,d#,, i.e. it equals the change in the power reflected by 
the cavity, which equals the change in the absorbed power. But eq. (11) 
shows further that the variation of the power at the detector may be even 
much larger than the absorbed power (*). The series expansion of Hj, as given 
by eq. (11) shows that, provided that | E — Hj| < Ej, the response of either 
a linear or a square-law detector is always proportional to the power variation 
in the detector arm. It is such a variation that we want to maximize. 

Let us now give a brief account of the most important situations which 
may occur using the above magie-T mounting. 


a) If the load in arm 2 is matched, then E,=0, and the variation of 
the power to the detector is proportional to H,dH,, i.e. it is one half of the 
reflected power from the cavity. The linear or square-law detector response 


(*) The system behaves as a phase-detector, since the signal to be detected is pre- 
viously mixed with the coherent signal H,. In this case the situation in the detector 
arm is illustrated by the following example. Consider two voltage sources of e.m.f. 
F and f respectively, and negligible internal resistances. If separately shorted on a 
resistance /, the two sources supply, respectively, powers equal to /?/R and f2/R, 
whereas if connected in series (or in opposition), they supply an extra power 2(fF/R) 
or — 2(ff/R). The simultaneous presence of the two sources in the circuit forces 
each of them to supply an extra power + (fF/R). 
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is then proportional to x’ to the first order. The same result holds strictly if 
a square-law detector is used and if the r.f. oscillator frequency and the re- 
sonance frequency of the cavity are made to coincide. In this case one maylock 
the r.f. oscillator to the resonance frequency of the cavity, using an A.F.C. 
system. This minimizes the effects of the cavity vibrations induced by the 
modulating field, which, being H -dependent and coherent with the modulating 
field, may give a shift in the base-line of the recorded spectrum (°). The 
optimum sensitivity is reached if the power reflected from the cavity is max- 
imized. This is obtained by the use of a coupling iris which renders the 
V.S.W.R. equal to 2 + 4/3 (**). It seems worthwhile to note, anyway, that 
the maximum is very broad and it does not necessarily coincide with an 
optimum signal-to-noise ratio, when a crystal detector is used, due to the 
anomalous dependence of the crystal noise on the r.f. power level. 


b) If an amplitude unbalance is left, namely if in eq. (11) is A, AH, 
and gp=0 or p= the power variation at the detector is proportional to 
(H, — E.) dH, and the detector response is proportional to 7" (!9!). In order to 
maximize the signal, one has to maximize both dH, and (E, + E.). The 
maximum of the reflected amplitude occurs for E,= 0, but for this value an 
inversion on the signal sign occurs (*): it is thus convenient to have H, small. 
though larger than dE, and £, large. But if a crystal detector is used, a com- 
promise in the bucking has to be made (°). Eq. (11) shows that the same results 
hold even if some phase unbalance is left, provided that the dg term is can- 
celled by keeping the cavity at resonance. 


c) If a small phase-unbalance is left, such that cos » 1 but p> dy, 
the system is sensitive only to phase variations and its response is propor- 
tional (11) to y’. Considerations similar to the ones of the above case, may 


be made concerning sensitivity. 


d) If complete balance is used (aside from the bridge stability problem) 
the response is in general not simple but it depends on both x’ and x" (1°). 


From the above considerations, it follows that if one is interested in 7’, 
the optimum sensitivity is obtained in case b); furthermore cases b) and c) 
require the same microwave mounting, which is a very important feature if 
one is interested in both 7’ and 7’. Conversely, also case a) may be sometimes 
convenient, due to the possible minimization of the effects of the cavity 
vibrations. 


42. The apparatus. — The block-diagram of our microwave apparatus 
(for X- and K-bands) is shown in Fig. 5: measurements corresponding to 
cases a), b), c) and d), may be made without any change in the mounting. 
Usually, most measurements are performed in conditions b) and ¢). 

The (medium power) klystron is immersed in a sealed, water-cooled oil 
bath, which provides proper shielding, cooling and temperature stability, and 


(16) N. BLOEMBERGEN, È. M. PurceLL and R. V. Pounp: Phys. Rev., 73, 679 (1948). 
(7) C. A. WHITMER, R. T. WEIDNER, J. §. Hsrane and P. R. WEISS: Phys. Rev., 
74, 1478 (1948). 
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which preserves the possibility of manual adjustment of the klystron CO, 
The highly stabilized power supply is suitable for high-voltage klystrons. The 
x-band klystron is followed by a ferrite isolator. The circuitry within the dotted 
Hewlett-Packard closed line in Fig. 5, which includes the magic-T No. 1, is used 
for the Pound I.F. automatic frequency control(?) when measurements in 
conditions d) and c) are wanted. In these cases, it is more convenient to lock 
the frequency on a cavity other than the one which contains the sample. in 
fact, in so doing, magic-T No. 2 may be unbalanced as one wishes, using the 
variable probe which is in one of the side-arms. 


stab magnetic - field ei sample RI Li 
| cavity {BB to amplifier, demodulator 
= E wattmeter & recorder 
non stab. klystron inser bolometer} 
& frequency shifter 
power ferrite measurement matched 
pera eS a ] 
detector | 
| and meter | 
| 
| | variable 
| i 11 [phase adjustable| | probe 
Î AME shifter cavity | 
| | 
| ) 
L-|-=f=------— : 
O 
30 Mc. oscillator 
& buffer 
amplifiers 
JA 
strip 
lock-in 
detector 
Fig. 5. — Block-diagram of the microwave instrumentation. 
eo oD 


When condition a) has to be actuated, the above mentioned microwave 
circuitry (magic-T No. 1) and the variable probe in magic-T No. 2 are no 
longer used. Magic-T No. 3 serves as a power divider, and the klystron fre- 
quency is now locked on the sample cavity. In this case, the intermediate- 
frequency A.F.C, in the Zaffarano-Galloway mounting (!8) has been found more 
convenient. 

A sperry self-balancing barretter bridge is used for measurements of the power 
incident on (and reflected from) the sample ‘avity. The measurement is per- 


(19) F. P. ZAFFARANO and W. C. GALLOWAY: Tech. Rep. No. 31, R.L.E., M.I.T., 1947, 
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formed using the directional coupler mounted on the sample cavity arm. A 
quantity proportional to the incident power may be measured by substituting 
the cavity by a short, or setting the cavity far from resonance. The measu- 
rement of a quantity proportional to the power reflected by the cavity at 
resonance allows the optimum sensitivity conditions to be reached; moreover, 
this measurement (or the one made on the directional coupler of magic-T No. 1) 
is convenient in order to lock the A.F.C. to the maximum of the cavity re- 
sonance (°). 

The frequency being locked on a cavity, the temperature drifts limit the 
frequency stability. Anyway, the stability which may be reached (of the 
order of one part in 10°) very largely fulfils our requirements. The frequency may 
be measured using a good secondary frequency standard. For most work in 
X-band, a good cavity wavemeter (having a precision better than 107!) is 
enough. However, a system comprehensive of a secondary frequency standard 
oscillator, frequeney multiplier, mixer and receiver will be available in the 
near future. 

Since the sensitivity is proportional to the unloaded @ of the cavity, high-( 
cavities are used, which are similar to the one designed and used by the authors 
for measurements performed at 
M.I.T. (7°), and described by STRAND- 

BERG et al. (*). They are precision 

machined cylindrical cavities, using 

the TE, mode. The sample holder, 

ag shown in Fig. 6 lies orthogonal to 

the cylindrical axis, and passes 7M,-mode lifting 
through the center of the cavity. groove 
Fig. 6 shows also a groove which is 

used in order to lift the degeneracy Sample holder 
with the TM,,, mode (2°). Since, at 
the center of the cavity, H" has a 
maximum and is directed along the 
axis of the cylinder, the methode 
descibed in Sect. 2 may be used. One 

extra advantage offered by this type of cavity, is given by the fact that the field 
modulation is in this case less effective in producing the mentioned vibrations. 

The absorption signal is given by the detector on magic-T No. 3. A erystal 
detector is usually used, but a bolometer may be used, particularly when 
absolute absorbed power measurements are required. 


Wave guide 


Coupling iris 
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Fio. 6. — High-Q resonant cavity. 


5. — Detection and recording of the magnetic absorption signal. 


A modulate absorption signal is obtained using a magnetic field modulation 
which is accomplished as described in Sect. 3°5. The modulated absorption 
signal is injected through a well shielded step-up transformer into a conven- 


(19) M. B. PALMA-VITTORELLI, M. U. Parma, D. Patumpo and F. SGARLATA: Nuovo 


Cimento, 3, 718 (1956). 
(29) C. G. MONTGOMERY (editor): Technique of Microwave Measurements, R.L.5.. 


No. 11 (New York, 1948). 
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tional frequency-selective amplifier (2!) having a few Hz bandwith and 
135 db gain, similar to the one described by STRANDBERG et al. (*). The trans- 
former makes the noise figure of the amplifier approach unity because at the 
same time it increases the effective impedance of the signal detector and the 
input signal amplitude. Its primary impedance is small, so that the short- 
circuit characteristics of the crystal detector hold. 

The amplified signal undergoes phase-detection in a Brown converter driven 
by an amplified and phase-adjusted signal given by the same oscillator which 
drives the field-modulation power amplifier. 

When a modulation is used which is small compared with the linewidth, 
the absorption signal is proportional to the derivative of the absorption curve. 
The search for the H-values for which the absorption maxima occur is made 
by a manual setting of the helical potentiometer which drives the magnetic 
field (see Sect. 3°3). The maximum condition is monitored by a central-zero 
meter at the output from the phase-detector. A plot of the whole spectrum 
is obtained by feeding the output from the phase-detector into a Speedomax 
recorder. The drive-shaft extension of the latter drives, through a gear, the 
shaft of the helical potentiometer which in turn drives the magnetic field. 
In so doing, linearity is preserved. 

In order to have a good fidelity when recording a spectrum, the modulation 
amplitude has to be small. This causes a decrease of the sensitivity, by a factor 
roughly equal to the ratio between the modulation amplitude and the line- 
width. If one is interested in the detection of the absorption only, the modu- 
lation amplitude may be raised up to the order of the linewidth. The sensi- 
tivity may also be raised using a larger time constant in the filter which follows 
the phase-detector, i.e. reducing the bandwith. Anyway, if a good fidelity is 
wanted, the recording time has to be in this case proportionally increased. 
A typical compromise is a recording time of the order of half an hour and a 
related overall bandwith which depends on the spectrum to be recorded. 

The ultimate sensitivity has been given in Sect. 1. 
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